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Preface 


The analysis of structures is a common core course requirement in U.S. universities (and those elsewhere) for 
students majoring in civil, mechanical, aeronautical, naval, architectural, and other engineering disciplines. The stu- 
dents are trained in structural analysis through a series of courses beginning with strength of materials and the analy- 
sis of simple indeterminate structures. In these elementary courses, the students are exposed to the fundamental 
structural analysis concepts in a simplified form. Comprehension of these principles becomes essential because these 
basic courses lay the foundation for other advanced structural analysis courses. The usefulness of this subject cannot 
be overemphasized because structural mechanics principles are routinely used in different engineering disciplines. 
We can even speculate that some of the concepts were used consciously or subconsciously for millennia by the mas- 
ter builders — the Romans, Egyptians, Eurasians, Chinese, Indians, and many others — of magnificent edifices, cathe- 
drals, temples, bridges, ships, and other structural forms. The analytical formulation of the principles, however, is 
much younger. It is popularly traced back to the cantilever experiment (ref. 1) of Galileo (1564—1642) depicted in 
this figure. 



Galileo's cantilever experiment. 


Even though some of his calculations were underdeveloped, Galileo’s genius is well reflected in the solution of 
the problem, especially because Newton (1642-1727), bom in the year of Galileo’s death, was yet to formulate the 
equilibrium laws and develop the calculus that we use in the analysis of structures. Industrial revolutions, successive 
wars, and their machinery requirements assisted and accelerated the growth of stmctural mechanics because of its 
necessity and usefulness in such cases. Several textbooks have been written on the subject, beginning with a compre- 
hensive treatment by Timoshenko in his Strength of Materials (ref. 2), first published in 1930. Therefore, the logic 
for yet another report on this apparently matured subject should be given. 
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We address this issue by considering the analysis of a symmetrical, four-legged table as an example (see Illustra- 
tive Example 15). Navier (1785-1836) attempted to determine the four reactions along the four legs of the table, 
but he could write only three equilibrium equations in terms of the four unknown reactions. He could not solve the 
(3 x 4) rectangular system of equilibrium equations because there were four unknowns but only three equations were 
available. Navier thus identified the indeterminate nature of structural mechanics problems. Their solution required 
additional conditions that were needed to augment the equilibrium equations. Deformation compatibility represented 
the additional conditions. If Navier or another researcher had completed the formulation of the compatibility condi- 
tions during the past century and a half, there would have been no need to write this report. As it turned out, the com- 
patibility formulation that was available in the literature remained either incomplete or ad hoc in nature both for 
structural analysis and the theory of elasticity. We have researched and understood these conditions, which for 
Navier’s table problem provide the additional (1 x4) rectangular compatibility equation. Coupling Navier’s (3 x 4) 
equilibrium equations and our (1 x 4) compatibility condition provides a system of four equations, which in matrix 
notation can be symbolized as [£]{/?} = {P}. The four equations with a (4 x4) square matrix [5] can be solved for the 
four reactions { R } given an external load { P } . Thus, the direct force determination method, which can complement 
currently available indirect methods, has now opened up. This formulation, which couples the equilibrium equations 
and the compatibility conditions to directly solve for the internal forces without an intermediate step of solving for a 
displacement or a redundant, is referred to in the literature as the Integrated Force Method (IFM) for structures. 

In elasticity, Saint Venant (1797-1886) formulated the compatibility conditions in the field of an elastic 
continuum, but he overlooked the conditions on the boundary, which have since been formulated by Patnaik. 1 Aug- 
mentation of Saint Venant’s field conditions with our boundary compatibility conditions completed the classical 
formulation, which in the literature is referred to as the Completed Beltrami Michell Formulation (CBMF). The IFM 
and CBMF are equivalent force and stress formulations, for structures and elasticity, respectively. 

Structural analysis, despite its immaturity with regards to compatibility, progressed, but only through indirect 
methods. These are the stiffness method, which uses displacement as the primary unknown, and the redundant force 
method, which treats redundant forces as primary unknowns. Both are indirect formulations because the internal 
forces (or the reactions in the case of Navier’s table problem) are not the unknowns of either method, but they are 
back-calculated — from displacements in the stiffness method and from redundant forces in the redundant force 
method. 


Stiffness Method 

Navier observed that the table had three displacement unknowns, these were the transverse displacement, w, and 
two rotations, 9* and 9 y (see fig. 1.1); and three equilibrium equations were available to him. The three equilibrium 
equations, when expressed in terms of the three unknown displacements, resulted in the set of three stiffness equa- 
tions. Solution of the three stiffness equations yielded the three displacements, from which reactions could be back- 
calculated. Generalization of this procedure, which is credited to Navier, became the popular displacement, or 
stiffness, method. 


Redundant Force Method 

In this method, one leg of the table was “cut” to obtain a three-legged auxiliary table. This three-legged auxiliary 
table had three reactions, which were calculated from the three equilibrium equations for the external load P. For the 


hn three-dimensional elasticity, the boundary compatibility condition can be written in terms of stress (a.t) on the boundary of an 
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auxiliary table, the displacement A p at the “cut” was back-calculated from the three known reactions. The solution 
process was repeated for a load R , referred to as the redundant force, in place of the reaction for the leg that was 
“cut,” and the displacement A R at the “cut” was obtained in terms of the redundant force R. Since the physical table 
had no real cut, the “gap” was closed (A p + A R = 0), which yielded the value of the redundant, or one reaction for the 
table problem. Thus, the table problem was reduced to an equivalent three-variable problem with two loads, consist- 
ing of a given external load P and a known reaction equal to R , which was treated as a load. This three-variable 
problem could be solved with the three equilibrium equations. Generalization of this procedure became the redun- 
dant force method, which was popular at the dawn of computer-automated analysis. Currently, for all practical pur- 
poses, the redundant force method has disappeared because it was cumbersome and had limited scope. 



Equilibrium equations and compatibility 
conditions in elasticity. 


The fundamental equilibrium and compatibility concepts of structural mechanics can be depicted in a pie dia- 
gram. The immaturity in the compatibility condition is represented by the shaded quarter. This portion was recently 
completed by Patnaik. It is true that indirect analysis methods, which included the popular displacement method, can 
be developed utilizing the information contained in just three quarters of the pie diagram, bypassing the shaded quar- 
ter. The direct Integrated Force Method, however, was developed utilizing all the information in the pie diagram. 

The direct formulation was envisioned by Michell (1863-1940), and it is described by Love (1863-1946, ref. 3) in 
the following quotation. 

“It is possible by taking account of these relations [the compatibility conditions] to obtain a complete 
system of equations which must be satisfied by stress components, and thus the way is open for a direct 
determination of stress without the intermediate steps of forming and solving differential equations to deter- 
mine the components of displacements.” 

Philosophically, it is not difficult to conceive of a deficiency in a solution that could be obtained without the 
explicit use of certain conditions, such as the compatibility conditions for structural problems. The following quota- 
tion by Todhunter (1820-1884, ref. 4) describes one such situation. 

“Important Addition and Correction. The solution of the problems suggested in the last two 
Articles were given — as has already been stated — on the authority of a paper by the late Astronomer 
Royal, published in a report of the British Association. I now observe, however — when the printing of 
the articles and engraving of the Figures is already completed — that they cannot be accepted as true 
solutions, inasmuch as they do not satisfy the general equations (164) of § 303 [ note that the equations 
in question are the compatibility conditions ]. It is perhaps as well that they should be preserved as a 
warning to the students against the insidious and comparatively rare error of choosing a solution which 
satisfies completely all the boundary conditions, without satisfying the fundamental condition of strain 
[ note that the condition in question is the compatibility condition], and which is therefore of course not 
a solution at all.” 


vii 
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Compatibility conditions are required for the analysis of the indeterminate problems of structural mechanics. 
Because the compatibility conditions were not fully comprehended, indeterminate analysis was attempted without the 
explicit use of the compatibility condition. Now that we understand the compatibility conditions, it is natural to use 
these conditions to analyze indeterminate problems. The use of compatibility systematizes indeterminate analysis, 
especially for problems with temperature variations and the settling of supports, because the initial strains are 
directly accounted for through compatibility. 

This report, which presents the IFM to solve indeterminate structural mechanics problems, does not duplicate 
any existing textbook. Because it complements indeterminate analysis, it should be read by undergraduate students in 
strength of materials and elementary structural analysis courses. This report should also be valuable to researchers 
who wish to work on the new method of forces or to complement their understanding of the compatibility condition 
of structural mechanics. 

Industry required solutions to structural mechanics problems. This demand was fulfilled through the two indirect 
methods. However, completing the structural analysis formulations to develop a direct method, such as the IFM, was 
not the central worry of industrialists. Industry’s demand for solutions (not methods), complemented by the availa- 
bility of computers, promoted the rapid growth of structural analysis via the indirect displacement method. This 
method is very popular in industry, and for the foreseeable future, it will remain the method of choice. The displace- 
ment method, however, may have entered an era of diminishing returns because of the intensive, worldwide research 
of this method for the past half century. Therefore, structural analysts can afford to take time off from the displace- 
ment method to research the little-explored method of forces, which was temporarily abandoned when engineers 
began to use computers to automate structural analysis. 

In summary, there is the single direct stress determination method, which is the Integrated Force Method. In 
contrast, there are several indirect methods wherein stress is back-calculated, such as the stiffness method, the redun- 
dant force method, and others, as listed in table I (also see app. A). The direct IFM should be learned because it can 
bring value-added benefits to the analysis, design, and testing of structures. 

The IFM for advanced analysis is already available in the literature (see the IFM references, refs. 5 to 32). 
Through this report we would like to introduce the IFM to undergraduate engineering students. Although written 
especially for this young group, this report can be useful to others interested in learning the new method. For the 
benefit of advanced readers, an elementary treatment of the IFM concept for finite element analysis is also included. 


TABLE I.— METHODS OF STRUCTURAL MECHANICS AND ASSOCIATED VARIATIONAL FUNCTIONALS 


Method 

number 

Method 

Primary variables 

Variational 

functional 


Elasticity 

Structures 

Elasticity 

Structures 


1 

Completed Beltrami-Michell 
Formulation (CBMF) 

Integrated Force 
Method (IFM) 

Stresses 

Forces 

IFM variational 
functional 

2 

Airy formulation 

Redundant force method 

Stress function 

Redundants 

Complementary 

energy 

3 

Navier Formulation (NF) 

Stiffness method (DM) 

Displacements 

Displacements or 
deflections 

Potential energy 

4 

Hybrid method (HF) 

Reissner method (RM) 

Stresses and 
displacements 

Forces and deflections 

Reissner 

functional 

5 

Total formulation (TF) 

Washizu method (WM) 

Stresses, strains, and 
displacements 

Forces, deformations, 
and deflections 

Washizu 

functional 


The subject matter of the report is presented in four chapters — Introduction, Basic Theory of Indeterminate 
Analysis, Solution of Indeterminate Problems, and Integrated Force Method and Dual Integrated Force Method for 
Finite Element Analysis. Chapter 1 introduces the three analysis methods: the Integrated Force Method, the displace- 
ment method, and the redundant force method. Chapter 2 develops the analysis equations — equilibrium equations, 
compatibility conditions, deformation displacement relations, and force deformation relations. Chapter 3 solves a set 
of 15 examples, some with thermal loads and settling of supports. The integrated and dual force methods for finite 
element analysis are introduced and illustrated in chapter 4. The first four appendices (A to D) provide classification 
of variables, superposition techniques, strength of material formulas, and sign conventions. Appendix E lists the 
important symbols used in this report. 


viii 
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Scientists 


Summary 


The theory of structures that evolved over the centuries (see the following figure) camouflaged a deficiency in 
the compatibility formulation. It was a deficiency that would block the growth of the primal method of force (or 
IFM) for analysis of indeterminate structures, causing it to split into the stiffness method and the redundant force 
method, as depicted in the figure on the next page. Our research on the compatibility theory has now restored the 
natural course of growth, leading to the Integrated Force Method for structures and the Completed Beltrami-Michell 
Formulation in elasticity. 
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Determinate structures 

(equilibrium) 


Indeterminant structures 
(equilibrium and compatibility) 



Compatibility barrier prevented extension of force method for indeterminate structures. 

IFM allows free movement between variables: from stress to displacement and vice versa. It can be specialized to 
derive the stiffness and other methods given in table I. In addition, the IFM variational functional can be modified to 
obtain the functionals of the other methods. The derived methods cannot outperform the primal IFM, which has also 
been demonstrated numerically. Textbooks on the strength of materials, theory of structures, elasticity, and related 
subjects need to be revised so that they incorporate the new information on the theory of compatibility. This report, 
which utilizes the new information on compatibility conditions to solve indeterminate problems, is the precursor of 
a strength of materials textbook. 
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Chapter 1 

Introduction 


Solving indeterminate strength of materials problems via conventional methods can be cumbersome and difficult in 
comparison to analyzing determinate structures. Conventional analysis methods include redundant analysis, wherein 
members are “cut” and the “gaps” are subsequently closed; superposition techniques; and solutions through the applica- 
tion of energy theorems. Some of these concepts are illustrated in appendix B. Even though such procedures can solve 
simple indeterminate problems, their generalization is not straightforward, especially for complex structures. Further- 
more, such procedures, though applicable for linear analysis, cannot be easily extended for nonlinear or dynamic analy- 
sis. After comprehending the conventional methods, readers will rightfully realize that these redundant-based ad hoc 
techniques can be problem dependent. Their generalization, when possible, can be cumbersome at best. 

Analyses of indeterminate problems, in addition to the equilibrium equations (EE), require the compatibility condi- 
tions (CC). The natural equilibrium concept common to the analysis of both determinate and indeterminate problems is 
straightforward; and its development, understanding, and use for analyzing structural mechanics problems is complete. 
However, until recently, the same could not have been said about the compatibility conditions. These CC, which are 
unique to the analysis of indeterminate structures, have been neither adequately researched nor utilized in analyses. The 
immature state of development of the CC can be considered to be the primary impediment to the conventional analysis 
of indeterminate structures. We have researched and now understand these illusive compatibility conditions. Our under- 
standing of the CC has systematized the analysis of indeterminate structures. Conceptually, such an analysis can be rep- 
resented by the following equation: 


Equilibrium equation 
Compatibility condition 


{Force} = 


[ Mechanical load 
[initial deformation 


( 1 . 1 ) 


A balance of the internal force and external mechanical load is achieved through the equilibrium equation, which 
forms the upper part of equation (1.1). The compliance of force and initial deformation is achieved through the compati- 
bility condition, representing the lower portion of equation (1.1). This equation, which bestows appropriate emphasis on 
equilibrium and compatibility, provides both necessary and sufficient conditions for determining force in an indetermin- 
ate structure. The displacement, if required, can be easily back-calculated from the force. 

In advanced finite element structural analysis, the new force method that has been developed is based essentially on 
equation (1.1), which couples the EE and CC. In the literature, this method is referred to as the Integrated Force Method 
(IFM). Research publications on IFM are listed in the IFM references under the headings Integrated Force Method — 
Basic Theory (refs. 5 to 10), Integrated Force Method— Design Optimization (refs. 1 1 to 18), Integrated Force 
Method— Elasticity (refs. 19 to 22), and Integrated Force Method— Finite Element Analysis (refs. 23 to 32). Analysis of 
indeterminate problems, the subject matter of this report, though simpler, also requires the same EE and CC concepts. 

To maintain consistency, we refer to this procedure, which utilizes EE and CC simultaneously to solve structural 
mechanics examples, as the IFM solution to indeterminate problems. The equilibrium equations can be transmuted to 
generate the popular displacement method of structural analysis. Likewise, the redundant force method can be devel- 
oped through a cumbersome transformation of a set of selected continuity conditions (members are “cut” and the “gaps” 
are closed) that bears some semblance to the CC representing the lower portion of equation (1.1). In other words, the 
IFM bestows simultaneous emphasis on the EE and CC, whereas the displacement method and redundant method 
emphasize the equilibrium and the compatibility, respectively. Since the IFM utilizes both EE and CC, this method can 
be developed systematically and can produce reliable solutions even for difficult structural analysis problems. 
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Background to Analysis Concepts 


Some of the fundamental structural analysis concepts are illustrated through the example of a table problem. A 
four-legged table subjected to a load P with eccentricities e x and e y is shown in figure 1.1(a). The problem is to deter- 
mine the four reactions (R v R 2 , /? 3 , R^) along the legs of the table. It is assumed that this symmetrical table, which is 
made of wood, is resting on a level floor made of a rigid material, such as stone. The distances between the legs along 
the x- and y-directions are la and lb, respectively. The cross-sectional area and the modulus of elasticity of the legs are 
A and E, respectively. The problem is referred to as “Navier’s Table Problem” because Navier (1785-1 836) was the 
first to attempt its solution. The problem is solved as Illustrative Example 15 in this report. Here, the basic analysis con- 
cepts are examined through the table problem. 

Three equilibrium equations can be written for the problem. These are obtained through the summation of the 
forces along the ^-direction, E.V = 0, and the summation of the moments about the x- and y-axe, s, "L X M = 0 and ZyV/ = 0, 
respectively. The three equilibrium equations can be written in matrix notation as 







V 

r 2 



1 

1 
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—a 

-b 
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r 3 
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e y P 
e x P 





*4. 


l X J 


The three equilibrium equations are expressed in terms of four unknown reactions. The four reactions cannot be deter- 
mined from the three equilibrium equations. Thus, the problem is indeterminate. The indeterminate nature of the table 
problem was recognized first by Navier. Solution of the indeterminate table problem required augmenting the three 
equilibrium equations with one additional condition. The compatibility condition represented the additional condition. 

The table problem can be solved following three different methods. These are (1) the IFM, which is the subject 
matter of this report; (2) the displacement method, which is currently popular; and (3) the redundant force method, 
which for all practical purposes has disappeared from current practice. 


Integrated Force Method 

In the IFM, the compatibility condition is formulated in terms of the four reactions as 

(^) /?1_/?2+/?3_/?4) = 0 (L3) 

The three equilibrium equations (eq. (1.2)) and the compatibility condition (eq. (1.3)) represent four equations that 
can be solved for the four reactions. The IFM, which calculates the force unknowns directly without any reference to 
displacement, represents the direct force determination method. The basic IFM concept is to couple equilibrium and 
compatibility (as shown in eq. (1.1)), to calculate the forces. The IFM could not be developed earlier because the gen- 
eration of the compatibility condition, as given in equation (1.3), had not been known to Navier (1785-1836) (see the 
Integrated Force Method — Basic Theory references, refs. 5 to 10), nor to other structural analysts, until Patnaik’s 
formulation. In IFM, displacements, if required, can be back-calculated from reactions, see Illustrative Example 15. 


Displacement Method 

Navier, who recognized the indeterminate nature of the table problem but did not formulate the compatibility 
condition, developed a displacement solution to the problem. He observed that each equilibrium equation is associated 
with a displacement. For example, the table problem with three equilibrium equations has exactly three unknown 
displacements (see figs. 1.1(b), (c), and (d)). These are (1) the transverse displacement w along the ".-axis, correspond- 
ing to the transverse EE (the first equation in equation (1.2)); (2) the rotation 0 y about the x-axis, corresponding to the 
first moment EE (the second equation in equation (1.2); and (3) the rotation 0 y about the y-axis, corresponding to the 
second moment EE (the third equation in equation (1.2)). The displacements (w, 0 y , 0 ) represent the dual variables of 
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(a) Indeterminate table problem. 



(b) Average displacement along z-axis. (c) Average tilt about x-axis. (d) Average tilt about y-axis. 



(e) Redundant force method. 


Figure 1.1.- Navier's table problem. 

the EE given by equation (1.2) (see app. A). Navier expressed the three EE given by equation (1.2) in terms of the three 
displacements by eliminating four reactions (R v R 2 , R 3 , R 4 ) in favor of three displacements (w, 0^, 9 V ). Navier’s equa- 
tions can be written in symbolic form as 
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Here, (i.j = 1, 2, 3) are the stiffness coefficients that can be calculated from the material and the geometrical proper- 
ties of the table. 

Solution of the stiffness equation (1.4) yields the displacement values (w, 0^, 0 v ). From the displacements, the reac- 
tions are recovered through back-substitutions. The displacement method is an indirect method because even when only 
the reactions are required, it is necessary to proceed through an intermediate step of forming and solving the stiffness 
equations. Even though Navier’s genius did not formulate the CC given by equation (1.3), he did give us the displace- 
ment method that is currently popular. 


Redundant Force Method 

In the redundant force method, one of the table legs, for example the first leg of the table, is “cut” as shown in 
figure 1.1(e). The resulting table is referred to as the auxiliary determinate structure. The auxiliary structure with three 
unknown reactions (R n , f? 3 , R^) is solved using the three equilibrium equations given by equation (1.2), without the first 

( P P P\ 

Ro , R% , R 4 J obtained for the external load P is repeated next for a redundant 
force in place of the first leg of the table (with P = 0) to obtain R^ t R %- , R ^ . Deformation A a at the cut is obtained 
for the auxiliary structure for the external load P (or r£ , R 3 , R 4 ) and for the unknown redundant force 
(or R ^ , R ^- , Rf- j . The deformation A a is a function of the unknown redundant force , which represents the dual 
variable of the compatibility condition given by equation (1.3) (see app. A). 

Deformation A^at the cut for the first leg only is also obtained for . The actual structure has no “cut.” This 
condition is restored by closing the gap through the following continuity condition of the redundant force method. 

A°+A i? =0 (1.5) 

Equation (1.5), which is a linear function of the redundant is solved to obtain the redundant force The reactions 
(Ry,R 2 , Rj, R 4 ) and displacements (w, Q x , Q y ) are back-calculated by using the value of !%_. A solution using the redun- 
dant force method is given in appendix B for a beam problem. 

The redundant force method was developed before the end of the 19th century by Clebsh (1833-1872), Maxwell 
(1831-1879), Castigliano (1847-1884), Mohr (1853-1918), and others. The IFM is compared with the displacement 
method and the redundant force method in references 29 and 30, respectively. Comparison of the three different meth- 
ods is not the objective of this report. The objective is to describe the IFM to engineering students, who are encouraged 
to learn all three methods first, then compare and contrast the three methods. The purpose of this report is to show that 
indeterminate IFM analysis can be developed in a systematic manner. The straightforward procedure does not use 
redundant analysis. The theory will be developed from basic principles, and it will be self-contained. No prior know- 
ledge of compatibility conditions or IFM theory is required. Discussion of the theory is limited to simple linear analysis 
problems suitable for strength of materials and elementary structural analysis courses. The IFM theory, however, is 
problem independent, and it can be easily generalized for solving difficult problems in linear, nonlinear, and dynamic 
regimes. For the sake of completeness, an elementary treatment of the IFM for finite element analysis is given in chapter 4. 

This report is written to complement the existing textbooks on strength of materials and elementary structural 
analysis. Indeterminate IFM analysis for elementary problems will be presented in the subsequent two chapters: 

Basic Theory of Intermediate Analysis and Solution of Intermediate Problems. The IFM for finite element analysis is 
introduced in chapter 4. Five appendixes have been added. Appendix A reviews the classification of variables and 
methods of structural mechanics, appendix B discusses the superposition principles for indeterminate problems, 
appendix C summarizes some standard strength of materials formulas, appendix D summarizes the sign conventions 
used in the report, and appendix E lists the important symbols used in the report. 
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Chapter 2 

Basic Theory of Indeterminate Analysis 

Equilibrium and compatibility are the two fundamental concepts of the structural mechanics theory. The EE, how- 
ever, are written in terms of forces, which can be the axial force, shear force, bending moment, or the torque. But the 
CC, in their original state, are expressed in terms of deformations, which can be elongations P a , curvatures p h , or rela- 
tive angles of twist P f . It is necessary to express the CC in terms of forces because it can then be coupled to the EE, 
which are already available in terms of forces. For this purpose, two more equation sets besides the EE and CC are 
required. These are the deformation displacement relations which are required to derive the compatibility conditions; 
and the force deformation relations that are used to eliminate forces in favor of deformations in the compatibility condi- 
tions. The four sets of structural mechanics equations are the 

(1) Equilibrium equations (EE) 

(2) Deformation displacement relations (DDR) 

(3) Compatibility conditions (CC) 

(4) Force deformation relations (FDR) 

The four basic relations in terms of the three variables (forces, deformations, and displacements) are de veloped 
and illustrated in this chapter. Solving indeterminate problems using these relations is the subject matter of chap- 
ter 3. The IFM solution yields forces and displacements. Stresses and strains that can be induced in a structure can 
be back-calculated from forces by using standard formulas that are summarized in appendix C. 


Equilibrium Equations 

Force balance is the central concept behind the equilibrium equations. The four types of forces— axial force F, 
shear force V, bending moment M, and torque T — at a point B for a beam of length l, depth d, and width b which is 
oriented along the abscissa of a Cartesian coordinate system are shown in figure 2.1 . The axial, or normal, force F is 
along the x-axis. The shear force V is along the y-axis. The bending moment M in the x-y plane is directed along the 
z-axis. Torque T in the y-z plane is directed along the x-axis. For simplicity and clarity, we will confine our discus- 
sions to two-dimensional problems only. Even two-dimensional problems utilize the z-direction to define moment. 
Standard sign convention is followed — that is, the forces are considered positive when they are directed along posi- 
tive axes (as described in app. D). 

The equilibrium equations for the four types of forces common to both determinate and indeterminate analysis 
can be generated as follows (see fig. 2.1): 


(1) The sum of all axial forces along the x-axis is zero: 

(2) The sum of all shear forces along the y-axis is zero: 

(3) The sum of all moments along the z-axis is zero: 

(4) The sum of all torques along the x-axis is zero: 


IF =0 

(2.1a) 

jt-axis 


o 

II 

(2.1b) 

j-axis 


= 0 

(2.1c) 

2-axis 


M 

II 

o 

(2. Id) 

x-axis 
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Enlarged 
cross section 


Figure 2.1 . — Forces in a beam. 



(b) Reactions as unknowns. 



(c) Beam divided into two elements, (T) and (2), and 
three nodes, Q], [2], and [3]. 



r a *1 I" a 

Shear force and moment Two moments 

(d) Selection of force unknowns. 


Figure 2.2. — Analysis of a propped cantilever beam. 
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Simultaneous application of all four EE given by equation (2.1) is seldom required in the solution of simple struc- 
tural mechanics problems. For example, truss analyses use equilibrium along coordinate axes; beam analyses require 
moment and shear equilibrium equations, and shaft torsion analyses use the equilibrium of torques. 

Formulation of the EE is illustrated by considering a propped cantilever beam that is fixed at A and simply sup- 
ported at B, as shown in figure 2.2(a). The beam, which has a span of (6 = 2d), is subjected to a transverse external 
load P at the center of the span. Equilibrium equations can be developed by considering either reactions or internal 
forces as the unknowns. Both approaches, which for the purpose of analysis are equivalent, are illustrated for the 
problem. The sign convention summarized in appendix D is followed. 

Approach 1: Reactions as unknowns . — This method considers the three reactions as the unknowns of the prob- 
lem. These are the reaction R A and moment M A at support A and the reaction R B at support B, as shown in the free- 
body diagram in figure 2.2(b). Summation of the forces along the y-axis and the moments along the 4 -axis at point B 
yields the two EE: 


XV = 0 r a +r b -p = 0 

V-axis 


IM = 0 

x-axis 


M a -{R a )( + — = 0 


In matrix notation, the two EE can be written as 


( 2 . 2 a) 


-1 

t 


or 

In equation (2.2c), [B] is a rectangular equilibrium matrix of dimension (m x n), where m = 2 is the number 
of rows or the number of EE, and n- 3 is the number of columns or the number of unknown forces; { F} is the 
unknown force vector of dimension n- 3; and {P} is the load vector of dimension m = 2. The components of the 
load vector {P} in the EE must be aligned along the positive directions. The load (-P) is along the positive y-axis 
and the moment (PCI 2) is directed along the positive z-axis (see app. D). Since the two EE given by equation (2.2b) 
are expressed in terms of three unknowns ( R A , R B , M A ), the EE alone cannot be solved to determine the three reac- 
tions. The problem is indeterminate. The degree of indeterminacy designated by r can be calculated as 

r = number of unknown forces ( n ) - number of equilibrium equations (m) 

- or (2.3) 

r = n-m = degree of indeterminacy 

The structure is one-degree indeterminate since r - 1 . One additional equation is necessary to solve the problem. 

This is the compatibility condition of the problem, which is discussed later. 

Approach 2: Internal forces as unknowns. — In this approach, the beam is divided into elements and nodes. 

The cantilever beam is divided into two elements, which are the circled 1 and 2, and three nodes (1,2, 3), which are 
shown in squares in figure 2.2(c). Typically, boundary points, load application points, change of member orienta- 
tions, and change of material constitute nodes. An element is a member connecting two nodes. For a beam element, 
the force distribution can be determined if two internal forces anywhere in its span are known. The pair of force 
unknowns can be a moment M and a shear force V or two moments (M p , M q ) as depicted in figure 2.2(d). The oppo- 
site direction is used for M p and M q because the moment equilibrium is satisfied, or M p - M q . Although both 


-1 0 
0 -1 


Ra 

r b 


M 


A) 


-P 

PC 


[B]{F} = {P} 


( 2 . 2 b) 


( 2 . 2 c) 
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force systems are equivalent, we prefer the two-moment system because both moments have the same unit of 
measure — force times distance. The two-moment system is further expanded in figure 2.3(a), where associated 
shear forces induced at nodes p and q are determined by taking moments at locations q and p, respectively, as 


V„ 



a 




M n 


(2. 


The moments (M p , M q ) and associated shear forces (V , V g ) for the element depicted in figure 2.3(a) satisfy both 
transverse and rotational equilibrium conditions. In this approach, the EE are written along the free displacement 
directions for all nodes in the structure. We avoid writing the EE along the restrained directions, such as at 
support A, or node 1 (along the ^-displacement and z-rotational directions), and at support B, or node 3 (along the 
y-displacement direction) in figure 2.3(c), because this process excludes the reactions as additional unknown force 
variables, thereby reducing the number of equations. Reactions at the restrained nodes are back-calculated from the 
forces. For this example, both cases — internal forces, as well as internal forces and reactions, as unknowns — are 
illustrated. 

Case 1— Internal forces as unknowns: Each beam element has two moments, and the total structure 
has four unknown moments (Afj, M 0 , M 3 , Mf). Rotational EE at the simple support B yield M 4 = 0, and the 



(a) Beam element with two unknown moments (M p , M q ). 


P 



Forces at node 1 . Forces at node 2. Forces at node 3. 

(c) Equilibrium at nodes. 

Figure 2.3. — Internal force system. 
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unknowns reduce to three moments: Afj, A/ 2 , and M 3 , as shown in the free-body diagram of the structure in figure 
2.3(b). Equilibrium equations are not written for the clamped node 1, which is fully restrained. Node 2 can displace 
and rotate along the y-axis and s-axis, respectively. Two EE that can be written for node 2 along the two displace- 
ments are as follows (see fig. 2.3(c)): 


Along the y -direction 


x^° 

y 


-¥ 


(M 2 -Ml) 

a 


P + - 


M, 


■ = 0 


a 


Along the z -direction 


^ A/ — 0 — > — M 2 + A/ 3 — 0 

z 


For node 3, the rotational equilibrium (M 4 = 0) was used earlier to reduce the moment unknowns from four (M j, M,, 
Af 3 , M 4 ) to three (M v M 2 , M 3 ). The transverse equilibrium at node 3 was avoided to exclude the reaction R B as an 
additional unknown. Node 3, in other words, provides no additional equilibrium equation. The two EE of the canti- 
lever beam can be written in matrix notation. 


1/a 

0 


-1/a 

1 



M, 

-1/a 

-1 

m 2 


m 3 _ 



(2.5) 


The EE given by equation (2.5) is one-degree indeterminate because two equations are expressed in terms of 
three unknown moments. One additional CC is required for the determination of the moments. Once the moments are 
known, the reactions at the clamped support and the simple support can be back-calculated by writing the EE along 
the restrained directions, such as along the transverse and rotational directions at support A and along the y-direction 
at support B, as follows: 


^ V at node 1 

y 

^ M at node 1 

z 

^ V at node 3 

y 


Ba J u *- m 'K o 

a 

-M ^ + A/j = 0 

m 3 


or Ra = 


or 


m 2 -M i 


R b — = 0 or Rr = 


M a =M x 

m 3 


B 


(2.6a) 


The reaction vector {%} in matrix notation can be written as 


Ra 


-1/ a 

11 a 

0 


AT] 

m a 

■ = 

1 

0 

0 

< 

M 2 ( 

r b. 


0 

0 

1 la 


m 3 


or in matrix notation 

{R} = [B s ] {F} 

where \R%\ is the equilibrium matrix required to back-calculate reactions from the internal forces [F\. 


(2.6b) 
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Case 2 ■ — Internal forces and reactions as unknowns: In this approach, the EE are written at all nodes — that is, 
restrained as well as free nodes. There are two EE for each of the three nodes (1,2, 3), or there are six EE in terms of 
seven unknown forces, consisting of four internal moments and three reactions (My M My M 4 : M A , R A , R B ). These 
six equilibrium equations follow: 


TV = 0 at node 1 

—1/ a 

Ha 

0 

0 

0 

-1 

0 ' 

' M x 


0 ■ 

TM = 0 at node 1 

-1 

0 

0 

0 

1 

0 

0 

m 2 


0 

TV = 0 at node 2 

1/a 

-Ha 

-Ha 

0 

0 

0 

0 

m 3 


-p 

TM = 0 at node 2 

0 

1 

-1 

0 

0 

0 

0 

m 4 

■ 

0 

TV = 0 at node 3 

0 

0 

Ha 

0 

0 

0 

-1 

m a 


0 

TM = 0 at node 3 

0 

0 

0 

1 

0 

0 

0 

Ra 

Rb. 


0 


(2.7) 


The EE at the nodes given by equation (2.7) represent a concatenation of equation (2.5), the boundary moment con- 
dition M 4 = 0, and equation (2.6a). 

Different choices of unknown forces neither increase nor decrease the degree of indeterminacy of a structure. 
For example, the propped cantilever beam is one-degree indeterminate, or r = 1 , for all three choices of force 
unknowns: (1) reactions (R A , R B , M A ), n = 3, m = 2, and r-n-m= 1; (2) internal moments (M ] , My Mf), n = 3, 
m = 2, and r = 1; and (3) moments and reactions (My M 2 , . . . , R B ), n = l,m = 6, and r - 1 . It is preferable to work 
with fewer equations (eq. (2.2b) or (2.5)) rather than a large number of equations (eq. (2.7)), especially for manual 
calculations. Thus, we can understand and appreciate the basic principles of indeterminate analysis without intensive 
numerical calculations. 


Compatibility Conditions 

For indeterminate problems, compliance of deformations (/’(|3 1 , |3 0 , . . . , (3 n ) = 0) is the central compatibility 
concept. The CC can be derived in two steps: 

Step 1 — Derive the deformation displacement relations (DDR). 

Step 2 — Eliminate the displacements from the deformation displacement relations to obtain the compatibility 
conditions. 

Step 1 — Derive the defonnation displacement relations: The DDR is an important structural mechanics relation. 
The DDR is the central ingredient behind both the equilibrium equations and the compatibility conditions. The EE 
that was generated earlier from the force balance principle can be alternatively derived from the DDR by using varia- 
tional calculus. Likewise, the reverse is true; that is, the DDR can be derived from the equilibrium equations. The CC 
can be generated from the DDR by direct elimination of the displacements (see app. A). 


Deformation Displacement Relations 

We will derive the DDR from the EE that have been formulated earlier. The deformations are associated with 
each type of force variable. Extension is the deformation for normal or axial force. Likewise, for the bending moment 
the curvature is the deformation. For shear force it is the shear deformation, and for torque it is the relative twist 
angle. In the derivation of the DDR, it is not essential to know the exact nature of the deformations. These will be 
further explained during the discussion of the force deformation relations. Here, it is sufficient to understand that half 
the product (force times deformation) represents the internal strain energy IE that is stored in a structure. The internal 
energy can be written as 
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( 2 . 8 ) 


IE = ±{F} T $} = ±(F 1 pj +F 2 p 2 + , . . . , + F„P„) 

where the deformations (P } , p„ . . . , p n ) correspond to the n internal forces (F p F 2 , ■ ■ ■ , F n ), respectively. Likewise, 
the work done W by the external loads can be written as 

W = ^{P} T {X} = \{P x X l + P 2 X 2 +, . . . , + P m X m ) (2.9) 

where the displacements (X V X 0 , . . . , X m ) correspond to the m external loads (P,, P 2 , . . . , P m ), respectively. In 
equations (2.8) and (2.9), deformation {p} and force {F} are n-component vectors, whereas displacement {X} and 
load {P} are m-component vectors (with degree of indeterminacy r = n- m). 

According to the work-energy conservation theorem, the internal energy (IE) stored in the structure is equal to 
the work done by the external load ( W ), or 


IE =W 

or (2.10) 

In equation (2.10), the load vector {P} can be eliminated in favor of internal forces [F] by using the EE 
(ffi]{P} = {P }) to obtain the following relation: 



{F} r {P} = {F}W(X} 

(2.1 la) 

or 

{F} r ({P}-[B] r {X}) = 0 

(2.11b) 


Because the n forces can be arbitrary and ( F } is not a null vector, its coefficient should be zero, which yields the 
DDR as 


{p} = [Z?f{X} (2.12) 

The DDR are easily defined through the equilibrium matrix [B], which is essential for analysis. Generation of the 
DDR does not require additional effort once the equilibrium matrix [ B ] is known. 

In the derivation of the DDR given by equation (2.12), it is assumed that load { P) and displacement {A} are in 
the same direction. The directions of the displacement components {A} in the DDR given by equation (2.12) are 
along the corresponding directions of load components {P}. Thus, orientation of the load vector { P} in the EE 
([B]{F} = {P}) along positive axes will yield displacements in the positive axes, see appendix D. 

Step 2 — Eliminate the displacements from the deformation displacement relations to obtain the compatibility 
conditions: In the DDR given by equation (2.12), n deformations (P) are expressed in terms of m displacements 
{A}. Elimination of m displacements from the n DDR yields r = n-m equality constraints as 

[C]{P} = {0} (2.13a) 

Equation (2.13a) represents the r compatibility conditions of the indeterminate structure with n force and m displace- 
ment unknowns. The compatibility matrix [ C] with r rows and n columns has full row rank r. 

The deformation { P } in the CC given by equation (2. 1 3a) represents total deformation consisting of an elastic 
component {p} e and initial component { P } () as 


{PMPr + {p}° 


(2.13b) 
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The CC in terms of elastic deformation can be written as 


[C]{P} = [C]{P} e + [C]{P}°-0 (2.13c) 

or [ c HP} e = { 8 ^} 

where {8/?} = -[C]{P}° (2.13d) 

The compatibility condition, when expressed in terms of total deformation, represents a homogeneous equation, 
such as equation (2.13a). The CC becomes a nonhomogeneous equation when it is written in terms of elastic defor- 
mations, such as in equation (2.13d). 


Initial Deformations and Support Settling 

Initial deformations {|3} 0 are included in the right side of the compatibility conditions through the effective ini- 
tial deformation vector {8 R] defined in equation (2.13d). Such deformations, when due to thermal effects, represent 
temperature strains, which can be written as the product of the coefficient of thermal expansion a and the tempera- 
ture change AT as e r = a AT. Initial deformation due to thermal strains ({(3}° = {el') can be easily calculated and 
included in the compatibility conditions. These are not further elaborated on here but are discussed in the solution of 
the examples in chapter 3. 

Initial deformations {p}° due to support settling can be calculated from energy considerations. Let us assume 
that a support settles by X amount and that the corresponding reaction induced at the support is R. The reaction that 
can be back-calculated from forces (F), by using equation (2.6b), can be rewritten for a single reaction R as 


R = {‘B r } T {F} 


(2.14a) 


where {misrepresents a row of the equilibrium matrix [fi-J in equation (2.6b). 

The work done by the induced reaction R and the prescribed displacement X can be written as 

W = --XR (2.14b) 

2 

The work term is negative because the reaction R and displacement X are in opposite directions. 

The strain energy IE can be written as 


where { {5}° represents initial deformation due to an X amount of support settling. 
The work-energy conservation theorem (W = IE) can be written as 


(2.14c) 


(2.14d) 
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The reaction R is eliminated between equation (2.14a) and (2.14d) to obtain 


{P}° +{®r} 2 A ' 2+- • • ->+{®r} p 

or {P}° =[»»■]{*} 

nxl nxp pxl 

Because {F} is arbitrary and it is not a null vector, 

{P}° = -{<B r }X (2. 1 4f) 

Equation (2. 14f) yields the initial deformation vector for the settling of a single support by X amount. The equation 
can be generalized when a p number of supports settle simultaneously by the amount X i , X 2 , • • • . X p as 

{P} 0 = {®r } ! ^1 + {®r} 2 *2+> • • ;«+{®r} p %p 
or {P}°=[®r]{*} 

nxl nxp p x 1 

m°=-[Zz] T {x} (2-15) 

where the p columns of matrix [® r ] correspond to p rows of the equilibrium matrix [fi_J in equation (2.6b) for the p 
number of prescribed displacements |x 1 due to the simultaneous settling of the p supports. Equation (2. 1 5) is simi- 
lar to the DDR ( { P> = [fi] 7 ! X) ) except that the EE are written for the restrained nodes. 

Null Property of the Equilibrium Equation and Compatibility Condition Matrices 

The following null property {[B][C] T - [0] or [C\\R\ T = [0]) of the equilibrium and compatibility matrices can be 
verified from equations (2.12a) and (2.13a). Elimination of the deformation between equations (2.12) and (2.13a) 
yields 

[C][Z?] r {X} = [0] (2.16) 

Because displacement {X} is arbitrary and it is not a null vector, its coefficient matrix must vanish, or 

[C][sf=[0] or [F][cf=[0] (2.17) 

For correctness, the null property of the equilibrium and the compatibility matrices should be verified after the gen- 
eration of the matrices. 
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Illustration for Compatibility Conditions 


The example of the propped cantilever beam used earlier to illustrate the generation of the EE is used again to 
illustrate the calculation of the CC. The two simple steps mentioned earlier, (1) generate the DDR and (2) eliminate 
displacements from the DDR to obtain the CC, are followed. The CC are generated for all three cases for which EE 
were developed earlier in equations (2.2b), (2.5), and (2.7). 

Case 1: Reactions (R A , R B , M A ) as unknowns (see eq. (2.2b )). — 

Step 1 — Derive the deformation displacement relations: These relations ({[}} = \B] T \X]) are obtained from 
equation (2.2b) as 


p!=-Al +« 2 

P 2 =-*l 

P 3 =-A 2 (2.18) 

where deformations (p 1? p 2 , P 3 ) are associated with the reactions (R A , R B , M A ), respectively. The displacements (Xj, 
X 2 ) are also referred to as the dual variables of the first and second EE given by equation (2.2b) (see app. A). 

Step 2 — Eliminate the displacements from the deformation displacement relations to obtain the compatibility 
conditions: Two displacements (A, and X 2 ) are eliminated between the three DDR given by equation (2.18) to obtain 
the single CC in deformations as 


Pi— p 2 +fp 3 =o 


or 


[c]{P} = [i 


4 


= o 


(2.19) 


Thus, 


[c] = [i -i e] 


The null property ([B][C] r ) for the EE and CC matrices can be verified as 





' 1 ' 


-1 

-1 

o' 


0 




-1 

— 


i 

0 

-1 



0 




e 


L J 


( 2 . 20 ) 


Case 2: Internal forces (M v M 0 , Mf) as unknowns (see eq. (2.5)). — For this choice of force variables (M v M 2 , 
M 3 ), the following CC is obtained by writing three DDR from the two EE given by equation (2.5) and then eliminat- 
ing the two displacements: 
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2Pi +P 2 + p3 


or 


[2 1 


fPil 


!] 


P 2 

P 3 . 


= 0 


(2.21a) 


Thus, 


[C] = [ 2 1 1 ] 


The reader can derive and verify the CC. The null property of the EE and CC matrices, \B}\C\ T = [0], can be 
verified as 





"2" 



1/ a 

-1/a 

-1/a" 

1 



"o' 

0 

1 

-1 

1 


0 







(2.21b) 


Notice that the elements of matrix [C] given by equation (2.21a) are dimensionless. This is because the 
deformations (pj, P-,, P 3 ), which are rotations corresponding to the moments (M j, M 2 , AT,), have the same unit of 
measure. In contrast, the units of measure of different elements of the [C] matrix given by equation (2.19) differ 
because deformations pj and p 7 , which correspond to reactions, are extensions measured in length units; whereas 
the deformation P 3 , which is due to a moment, is a rotation that is dimensionless. 

Case 3: Internal forces and reactions as unknowns (see eq. (2.7 )). — Determination of CC for this selection 
offerees also starts with the formulation of the seven DDR ({p} = \B] T {X\) expressed in terms of six displace- 
ments (see eq. (2.7)): 


Pi = — x 2 +~ P 4 ~ x 6 

a a 

p 2 =— —+X 4 P 5 = X 2 

a a 

P 3 =- — +— P6 =-^l 

a a 

P 7 =-*5 ( 2 - 22 ) 

Elimination of the six displacements (Xp X 0 , . . . , X 6 ) from the seven deformation displacement relations 
given by equation (2.22) yields the single CC in terms of the seven deformations (pj, p o , . . . , p 7 ) as follows: 
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(2.23a) 


Thus, 


[21102 -1/a 1/a] 



P 7 

a 

a 


V 


P 2 


P 3 

1 /a} 

P4 


P 5 


Pe 


P 7 . 


= 0 


[C] = [2 1 1 0 2 -1/a 1/a] 


(2.23b) 


The null property of the equilibrium and compatibility matrices can be verified from the EE and the CC, which 
are given by equations (2.7) and (2.23b), respectively, as 


[B] [Cf 

6x7 7x1 


'O' 

0 

0 

0 

0 

0 

6x1 


(2.24) 


The problem has only one compatibility condition, irrespective of the choice of force unknowns, such as reactions, 
internal forces, and their combinations. However, the EE and CC matrices, [ B ] and [C], respectively, differ depend- 
ing on the choice of force unknowns. 


Force Deformation Relation 

The equilibrium equations are expressed in terms of force variables {F} as [B]{FJ = { P}. Likewise, the compat- 
ibility conditions are written in terms of deformations { (3 } as [d{|3} = {0}. Since indeterminate analysis requires the 
coupling of the EE and CC, it is necessary to establish relations between forces and deformations. The force defor- 
mation relation (FDR) can be used to express the CC in terms of forces, which can then be coupled to the EE, which 
are already available in terms of forces. The FDR of strength of materials is equivalent to the familiar Hooke’s law of 
elasticity, which relates stress a to strain e through the Young’s modulus E of the material (o = Et). The FDR can be 
obtained from Hooke’s law by relating stress to force and deformation to strain. For a normal, or axial, force F acting 
in a bar with a cross-sectional area A and length l, as shown in figure 2.4, the FDR can be obtained as follows: 

F 

o = — and e = — (2.25) 

At 
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Hooke’s law can be written in terms of force F and deformation (3° as follows: 


or 


e = l! = £ = -L 

£ £ E AE 

(2 - 26 *> 


where the flexibility coefficient is defined as 


g = 


t 

AE 


(2.26b) 


The flexibility coefficient (g = i/AE) represents the deformation in the bar for a unit value of the force (F = 1). 

Force deformation relation from energy considerations . — The FDR given by equation (2.26a) can be alterna- 
tively derived from energy considerations. The first derivative of strain energy U with respect to force F is equal to 
the deformation (3° corresponding to that force F, 

l I 

or (2.27) 

oF 

Before equation (2.27) can be used, the strain, or internal, energy has to be defined. For the case of the normal force 
acting in a uniform bar of cross-sectional area A (as shown in fig. 2.4), the strain energy can be defined as 

i 

U=\ — Adx (2.28) 

l 2 
0 

When stress is eliminated in favor of force as o = FI A and strain in favor of force as e = F/AE, the strain energy for a 
uniform bar of area A can be obtained only in terms of force F as 


U = 


f p2 j p2e 

dx — 

J 2 AE 2 AE 

0 


or 



(2.29) 


Notice that equation (2.26a), which is derived from Hooke’s law, and equation (2.29), which is obtained from energy 
considerations, are identical. Experimentation is the basis of Hooke’s law. The alternative, energy-based derivation, 
however, is analytical in nature even though the definition of strain energy requires material constants such as the 
Young’s modulus E. In the analytical indeterminate analysis, we prefer the energy-based derivation for FDR because 
in this approach all four analysis equations (EE, CC, DDR, and FDR) can be derived analytically. 
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Strain Energy Expression for Force Deformation Relations 


Strain energy is a powerful scalar quantity that can be used to derive various analysis equations and formulations 
of structural mechanics (see app. A). In this elementary treatment, a definition of the energy scalar that is sufficient to 
derive the force deformation relations for beams will be presented. For the discussion here, the strain energy scalar 
can be defined as follows: 


U- 


2 

f ° 

+ o 

J 2 E 


ocA T 


dv 


(2.30) 


where 


E 

Young’s modulus 

AT 

change in temperature 

V 

volume of the beam 

a 

coefficient of thermal expansion 

a 

stress in the beam 


Equation (2.30) is specialized next to obtain the explicit strain energy expression for axial force, bending moment, 
and shear force, as well as for pure torsion. 

Case 1: Axial force in a bar . — For the bar shown in figure 2.4, the strain energy scalar given by equation (2.30) 
can be specialized with the following definitions: 


Elemental volume, dv = Adx 


F 

0 = — 
A 



(2.31) 


where U a is the strain energy stored in the bar subjected to an axial force F and a AT change in temperature. 
For a uniform bar of length l and cross-sectional area A, the strain energy U a is simplified to 


F f 

U a = +FaATC 

2AE 


(2.32) 


The axial elongation in the bar is obtained as the first derivative of the scalar U a with respect to the axial force F 
as 


y ^‘{ii) F * a&TI <233) 

In the absence of temperature (AT = 0), the force deformation relation given by equation (2.33) simplifies to equation 
(2.26a). 


NASA/TP— 2004-207430 


18 



Area A 


□ 

Cross section 
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Figure 2.4. — Elongation, p a , due to 
axial force F in the bar. 


Case 2: Bending of a beam. — The bending of a beam (see fig. 2.5(a)) involves the interaction of a bending 
moment and a shear force V (as shown in fig. 2.5(b)). However, in the energy expression it is customary to include 
only the strain energy due to bending moment. Because the strain energy due to shear force is small, it is neglected 
for simple structural mechanics applications without an appreciable adverse effect. By considering the beam shown 
in figure 2.5(a) as an example, we can obtain the strain energy in flexure from equation (2.30) with the following 
specialization. 

Stress is calculated from the standard flexure formula, see appendix C. Stress o at a location x along the beam 
length and at a distance y from the neutral axis shown in figure 2.5(c) can be written as 


S-~ — 
y ~ i 


(2.34) 


where 7 is the moment of inertia of the beam cross section and 'Mt is the bending moment. The script TM. is used for 
moment function at x as Ttt or %(x). The incremental volume dv for the beam with a uniform width b shown in figure 
2.5(a) can be written as 


dv = bdy dx (2.35) 

where dy and dx are incremental lengths along the depth and length of the beam, respectively. 

When we substitute the stress from equation (2.34) and the incremental volume from equation (2.35) into the 
energy expression given by equation (2.30), the following equation is obtained for the strain energy \J b in flexure: 


U b 



Y 


d/2 

J EaATbydy 


\dx 


0 {-d/2 


(2.36) 
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(a) Loads in the beam. 



(b) Moment 7K and shear V. 



(c) Bending stress distribution, a = TKyll. 
Figure 2.5— Simply supported beam. 


where d, b, and l represent the depth, width, and length of the beam, respectively. Bending moment for a beam shown 
in figure 2.3(a), which produces compression in the top surface and tension in the bottom surface, gives rise to con- 
cave, deformed shape. However, positive temperature at the top surface (and negative temperature at the bottom 
surface) produce convex deformation in the beam. To account for the two opposite types of deformations, a negative 
sign is introduced in the strain energy expression in equation (2.36). 

The moment of inertia is defined as 


dll 

/= J by 2 dy 

-dll 


(2.37) 


Moment ( M T ) due to thermal effect can be defined as 


dll 

Mf - J E(aAT)bydy (2.38) 

-dll 
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The strain energy due to bending can be simplified by substituting equations (2.37) and (2.38) into equation (2.36) as 
follows: 


U b 


J 


u 

' M 2 ' 

I 2 

l EI ) 



(2.39) 


Deformation p* due to a bending moment M is obtained as the first derivative of the flexural strain energy given 
by equation (2.39), with respect to that bending moment M, as 


b _ dU b _ r f M t <m \ 

P ~ dM _ J 1 £7 dM El dM j 
0 


(2.40) 


Likewise, deformation p fc '' due to shear force V (which produces bending, thereby contributing to the strain energy) 
is obtained as the partial derivative of the flexural strain energy with respect to the shear force V as 


b - s _ du b _ r aag M T dm ] 

p j \ei1)v ei avj 

0 


(2.41) 


Case 3: Torsion of a uniform circular shaft . — We can obtain a strain energy expression for a uniform circular 
shaft under torsion by specializing the strain energy formula given by equation (2.30) as follows: 


(1) The Young’s modulus E should be replaced by the shear modulus G. 

(2) The moment of inertia I should be replaced by the polar moment of inertia J. 

(3) The normal stress o should be replaced by the shear stress x. 

(4) The thermal coefficient a is set to zero because, for an isotropic material, the temperature effect does 
not produce thermal shear strain. 


The shear stress formula for the circular shaft shown in figure 2.6 can be written as 


x _ T 
r J 


(2.42) 


where x represents the shear stress at a distance r from the neutral axis and *T is the torque. Torsion as a function of x 
is represented by (T= T(jcj). Strain energy due to torsion can be written as 


t y 

X H 


c 



n 




J 


L X- 

- — a *-j 





Shear stress distribution (t = Tr/J) shown on enlarged cross section X-X. 
Figure 2.6. — Uniform circular shaft subjected to torsion T. 
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u s 



t 0 

r R o, 

1 

{—■ 

f 271 r dr 

2 

J JG 

J J 


0 

lo J 


dx 


(2.43) 


because 


J = 2nr dr 


I 

-iJ 


1 r T“ 

U s = ~ f dx 

JG 

o 


(2.44) 


Deformation due to a torque T can be obtained as the derivative of the strain energy due to torsion with respect to 
that torque T as 


dU s 


= f- 

J JG 


T 3T 


ar J jg dr 
0 


dx 


(2.45a) 


For a uniform shaft under constant torsion, the deformation in equation (2.45a) can be specialized for T — T as 


or 


T( 

B* = — 
JG 

P s =g s T 


(2.45b) 

(2.45c) 


where g . is the flexibility coefficient due to torsion and it is defined as 


g s = 


JG 


(2.45d) 


Illustrative Examples 

The deformations induced by axial and torsional effects can be calculated by direct substitutions of axial force F, 
torsional moment T, and temperatures in equations (2.33) and (2.45b), respectively. For flexure, the deformation cal- 
culations are also straightforward, but the coupling of bending moment and shear force has to be considered. Flexural 
deformation calculations are illustrated for the beam shown in figure 2.7. As mentioned earlier, beam response requires 
two internal unknown forces that can be either (1) two bending moments (M^, M 2 ), as shown in figure 2.7(a), or (2) a 
bending moment ( M) and a shear force (V), as shown in figure 2.7(b). We will derive the deformations for both cases. 
For simplicity, the temperature effect will not be included here, but it is included in the examples given in chapter 3. 

Case 1: Two bending moments (M v M 2 ). — For the beam shown in figure 2.7(a), M { and M 2 are considered to be 
the two unknown moments at the ends (A and B) of a beam of length l, respectively. The reactions can be calculated 
from the equilibrium conditions as 


Ra - -Ri 


Mo -M 


(2.46) 


The moment at any point x in the beam axis can be written as 


Mo -Mi „ x s ,, x 

W(x) = M x + 2 .. — ±x = M l Q.--) + M 2 - 


(2.47a) 
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(a) Moments M- \ and M 2 are the force variables. R A and Rg are 
the equilibrating reactions that represent shear forces. 



(b) Moment M and shear V are the force variables. M A and R A are 
the equilibrating reactions. 

Figure 2.7. — Deformation in a beam. 


The derivative of the moment 7K with respect to My is 

3^(x) x 

3Mj 1 r 

From equation (2.40), we can write the deformation due to the moment My at location A as 




nb J±YMl + Mi 
PA l £/A 3 6 


The deformation p* represents the rotation at location A due to the action of both bending moments (M ] 
wise, the deformation at location B can be calculated as 


b dU b £ f My | M 2 


n = 


dM 2 EI{ 6 3 
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(2.47b) 


(2.48a) 


, A/ 0 ). Like- 

(2.48b) 



The two flexural deformations (pj( , pg) can be written in matrix notation as 


[p£l 


"2 

r 

[«, 1 

IpbJ 

| ~ 6£7 

1 

2 ’ 

\m 2 j 


(2.48c) 


{P } b =\G] b {M} 


(2.48d) 



'2 1 
1 2 


(2.48e) 


where the (2 x 2) coefficient matrix [Gf is referred to as the flexibility matrix for the beam when M, and M 2 are con- 
sidered to be the moment unknowns. In the FDR given by equation (2.48c), the deformations (p^ , Pb ) follow the sign 
convention for associated moments (M j, M 2 ), which have opposite directions. 

Case 2: Bending moment M and shear force V. — For the beam shown in figure 2.7(b), M and V are considered to 
be the unknown moment and shear force at location B of the beam of length L The reactions (R A , M A ) at location A 
shown in figure 2.7(b) can be calculated from the equilibrium conditions as 

R a = V and M A = M + V( (2.49) 

The moment at any point x and its derivative with respect to the shear force V can be written as 

M(x) = M+ V(£ - x) 


dM{x) 

__ 


(f-jc) 


The deformation due to the shear force V at location B can be obtained from equation (2.41) as 


(2.50) 


or 



(2.51) 


The deformation p^ due to the shear force V represents the displacement at B along the shear force V. Likewise, the 
deformation due to the bending moment at location B can be calculated as 


or 






vr 


El 


+ M( 


The deformation p^ represents the rotation at B. 


(2.52a) 
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In matrix notation, the two deformations can be written as 


= C (2 72 \ M ~ 

[p£j"U/J//2 e 3 /3 L y _ 

or {P} = [G] to {^} (2.52b) 


where 


[( G] bs 


j_r i 1 1 / 2 

El[e 2 / 2 i 3 /3_ 


(2.52c) 


The (2 x 2) coefficient matrix [G] /w is the flexibility matrix of the beam corresponding to the choice of force variables: a 
moment and a shear force (M, V). 


Flexibility Matrix 

The relationship between the force variables and the deformation variables is established through the flexibility 
matrix [G] as 

{P} = [G]{F} (2.53) 

The flexibility matrix [G] is a symmetrical matrix of dimension ( n x n), n being the number of force unknowns. For 
axial force F and torque T, the flexibility matrices become (lxl) matrices defined by equations (2.26b) and (2.45d), 
respectively. For beam flexure, which involves two force variables, the flexibility is a (2 x 2) matrix. For the choice 
of two bending moments the flexibility matrix is defined in equation (2.48e). An alternative flexibility 

matrix for a bending moment and a shear force is defined in equation (2.52c). 


Summary 

Analysis of an indeterminate problem requires the coupling of the equilibrium equations and the compatibility con- 
ditions. The equilibrium equations can be obtained as a vectorial summation of internal forces and external loads. The 
compatibility conditions can be calculated in two steps: First, the deformation displacement relations are obtained. Next, 
displacements are eliminated from the deformation displacement relations to generate the compatibility conditions. Such 
compatibility conditions expressed in terms of deformations are rewritten in force variables by using the force defor- 
mation relations. The force deformation relation can be obtained from strain energy considerations. Together, the equi- 
librium written in forces and the compatibility, also expressed in terms of forces, represent a sufficient number of 
equations for determining the internal forces. Displacements, if required, can be back-calculated from forces. 
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Chapter 3 

Solution of Indeterminate Problems 


In this chapter, the theory presented in the previous chapter is used to solve a number of indeterminate problems, 
some with thermal loads and support settling. The basic steps that we will use follow: 

Step 0 — Solution strategy. 

At the initial problem-formulation stage, the sign convention is specified, the force unknowns are identified, and the 
displacement components are specified. The number of equilibrium equations, compatibility conditions, and degree of 
indeterminacy are determined. 

Step 1 — Formulate the equilibrium equations. 

Step 2 — Derive the deformation displacement relations. 

Step 3 — Generate the compatibility conditions. 

Step 4 — Formulate the force deformation relations. 

Step 5 — Express the compatibility conditions in terms of forces. 

Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve for 
the forces. 

Step 7 — Back-calculate the displacements, if required, from the forces. 

The solution of indeterminate problems requires the inversion of a square matrix, which, though trivial with a com- 
puter, can become cumbersome for manual calculations. Since the objective here is to master the basic concepts through 
manual solution, treatment will be confined to simple problems with small matrices. When the calculations become 
involved, a computer can be used for the solution. The IFM for computer solutions is introduced in chapter 4. 


Illustrative Example 1: Thermomechanical Solution for a Fixed Column 

A column of length 31, shown in figure 3.1(a), is restrained at both ends. It is made of steel with a Young’s modulus 
E of 30 000 ksi and a coefficient of thermal expansion a of 6X10 -6 per °F. The cross-sectional area of its central span 
(2A) is twice that of its boundary spans, which have an area of A. Solve the problem for forces and displacements for 
l = 10 in., A = in. 2 , and the following three load cases: 

Load case 1: A mechanical load (P , = 10 kips and P 7 = 20 kips) applied at the one-third and two-thirds span 
locations, as shown in figure 3.1(b) 

Load case 2: A uniform temperature variation (AT = 2000 per °F) along the central one-third span, as shown in 
figure 3.1(c) 

Load case 3: A uniform temperature variation (AT = 2000 per °F) along the entire length of the column, as shown 
in figure 3.1(d) 

Step 0— Solution strategy: The coordinate system is marked in figure 3.1(a), with the x-axis along the length of the 
column. The model for the column, consisting of three axial bars and two nodes, is shown in figure 3. 1 (e). The problem 
has three force unknowns (Fj, F„ F 3 ), one for each of the three bars, and two unknown displacements (X v X 2 ), one for 
each node, along the column x-axis. The column is one-degree indeterminate (r = n-m= 1) because there are three 
unknown forces (n = 3) but two unknown displacements (m = 2). The column requires two EE and one CC for the deter- 
mination of the three forces. To solve the problem, we follow the seven steps. 
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(f) Equilibrium conditions. 


Figure 3.1 . — Thermomechanical analysis of a three-span column. 
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Load case 1: Solution for mechanical loads . — 

Step 7— Formulate the equilibrium equations: The two EE of the problem are obtained by summing forces along 
the displacement directions (Xj, X 9 ), see figures 3.1(e) and (f). 

EE along Xj , F\ - F 2 - P\ = 0 (3.1.1a) 

EE along X 2 , ?2 - h ~ p 2 = 0 (3 - Llb) 

The EE in matrix notation ([B]{F} = (Pj) can be rewritten as 


-1 1 
0 -1 



(3.1.1c) 


Because the two EE represent an insufficient number of equations to determine three unknown forces, one additional 
CC is required. 

Step 2 — Derive the deformation displacement relations: The DDR, which are the main ingredients of the CC, are 
obtained using the EE matrix as ({(3} = [fi] r {X}). For example, the kth DDR can be obtained as the dot product of the 
let h column of the EE matrix [B] and the displacement vector { X} . The DDR for the problem has the following explicit 
form: 




p2 — ^"i — %2 


P 3 =*2 


(3.1.2) 


Step 3 — Generate the compatibility conditions: The single CC ([CliPl = 0) for the problem is obtained by 
eliminating two displacements from the three DDR as 


Pi +P2 +P3 - 0 


or 

Pi' 

[1 1 1] p 2 =0 

p 3 . 

(3.1.3a) 

or 

[C] = [1 1 1] 



The CC given by equation (3.1.3a) constrains the total elongation of the three bars to zero, which for the problem could 
have been asserted by observation. Correctness of the CC can be verified from its null property ([£][C] 7 ' = [0]). For the 
problem. 


-1 

0 


(3.1.3b) 
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The deformation {(3} in the CC given by equation (3.1.3a) represents the total deformation. The total deformation is 
the sum of an elastic component (fl'O and a thermal component ((30 as ((3 = (3 e + (3 r ). For mechanical loads only, (3 = \5 e 
because p' = 0. The CC given by equation (3.1 .3a) is written in terms of deformations, whereas the EE given by equa- 
tion (3.1 .lc) are expressed in forces. To couple the EE and CC, we need to write the latter condition in terms offerees; 
for this purpose, the force deformation relations are required. 

/ F£\ 

Step 4 — Formulate the force deformation relations: The FDR for the three bars for axial force I p = I, given by 
equation (2.33), are obtained as 


Bar 1, i = 10, A = 1; or 



Bar 2, l = 10, A = 2; or 



Bar 3J = \0,A = l; or 



(3.1.4) 


Step 5 — Express the compatibility conditions in terms of forces: The CC given by equation (3.1.3a) can be 
expressed in terms of forces by using equation (3. 1 .4): 


f[> 1/2 1] 


h 

A 




(3.1.5) 


In the homogeneous CC given by equation (3.1.5), the coefficient (10/E) can be set to unity without any consequence. 

Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve for 
the forces: The EE given by equation (3.1.1c) are coupled to the CC given by equation (3.1.5) to obtain the IFM equa- 
tions for the problem: 


■-1 

1 

O ' 

a 


-P\ 

0 

-1 

1 

h 

• = - 

- p 2 

1 

1/2 

1 

A 


0 


Solving equation (3.1.6) for P { = 10 kips and P 0 = 20 kips yields the forces as 


(3.1.6) 


a 

case 1 

14' 

pi 

= * 

4 

A 


-16 


(3.1.7) 


Step 7—Back-calculate the displacements, if required, from the forces: The two displacement components can be 
calculated from the DDR and FDR given by equation (3.1.2) and (3.1.4), respectively. Any two of the three DDR can 
be used for calculating displacements from forces. We use the first and third deformations in equation (3.1.2) to calcu- 
late the two displacements as follows: 
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(3.1.8) 


x l =-Pj =-1^L = -4.67x 10 -3 in. 

E 

X 2 = p 3 = = -5.33 x 10 -3 in. 

Because of the compatibility compliance, the satisfaction of the remaining DDR (P 2 =X l -X 2 ) can be verified. 

P, = ^_ = 0.667 xl(T 3 in. 

E 

and Xj -X 2 =0.667 xl0~ 3 in. (3.1.9a) 

Thus, P 2 — -^1 — (3.1.9b) 


The displacements (Xj, X,) are along the negative x-axis, or along the loads (Pj, P 2 ), as expected. 

Loadcase 2: Central span thermal toad.— Solution for thermal effect follows the procedure for mechanical loads 
with nontrivial {8P} in CC. Both thermal and mechanical deformations have to be included in the definition of total 
deformation. Deformations including thermal effects (case 2, when the temperature increases in the midspan by 
2000 °F) are as follows: 

Thermal deformations, 

P(=0 p 2 -aAT l = 0.12 p3=0 (3.1.10) 

0 

0.12 • 
t 0 



Total deformations (P = P e + PO, 


R case2_f£n 10*1 
Pl UpJl E 

p case2 = lH + 0n 

E 

ocase2 _ 10 *3 

P3 E 


(3.1.11) 
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The effective initial deformation vector can be calculated as 


{5*} = -[C]{p'} 


= -[l 1 1] 


0 

0.12 

0 


= {- 0 . 12 } 


(3.1.12) 


The CC given by equation (3.1.5) can be rewritten incorporating the thermal effect {8R} as 



1/2 



= {- 0 . 12 } 


(3.1.13) 


The left side of equation (3.1.13) is identical to equation (3.1.5), but its right side is replaced by the equivalent thermal 
load (8R = -[C]{P r }). Earlier in the homogenous CC in equation (3.1.5), the coefficient (10 IE) was set to unity without 
any consequence. However, in the presence of thermal loads, the coefficient has to be retained since the CC in equation 

(3.1.13) is not homogeneous. The solution for thermal load, case 2, is obtained by incorporating the following changes 
into equation (3.1.6): 

P j =P 2 = 0 since there are no mechanical loads. 




■-1 

1 

0" 

>1 


0 
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-1 
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■ = ■ 
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case 2 

144' 

or 

* 

F2 

- = — ■ 

144 



A 


144 


kips 


(3.1.14) 


(3.1.15) 


Step 7 — Back-calculate the displacements, if required, from the deformation displacement relations: Displace- 
ments are back-calculated from any two of the three DDR as 

10F, 

Xx =-Pj = =0.048 in. 

E 

X 2 = P 3 = — ^ = -0.048 in. (3.1.16) 

E 

In the displacement calculation, thermal strain is bypassed because deformation in the first and third bars, P, and P 3 , are 
not explicitly affected when the central span temperature is increased. 
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Let us calculate deformation p o by using the FDR (eq. (3.1.11)) as well as the DDR (eq. (3.1.2)). 

From the DDR equation (eq. (3.1.2)), 

P 2 =X 1 -A r 2 or p 2 =0.096 (3.1.17a) 

From the FDR equation (eq. (3 . 1 . 1 1)), 

ft-, =^- + 0.12 = -0.024 + 0.12 = 0.096 (3.1.17b) 

E 

Deformation (P 2 ) calculated from both the DDR and FDR agreed, as expected. The central span expands because of 
an increase in the temperature, resulting in a positive displacement (X j) at node 1 and a negative displacement (X 2 ) at 
node 2. Because of symmetry, |Xj| = |X 2 |. 

Load case 3: Uniform thermal load . — Uniform thermal increase (A T = 2000 °F), case 3, can be solved by follow- 
ing the procedure shown for thermal load, case 2. Deformations due to thermal effects, case 3, are as follows. 

P case 3 = Bfl + (oA7T) ] = -^ + 0.12 

p case3 = ^2. + 012 
- E 


P ease 3 
3 



The CC and the effective initial deformation vector {§R} are calculated as follows: 


10 

E 


[1 1/2 1 ] 


h 


= -0.36 


The IFM equation (3.1.15) is modified by incorporating 8 R for case 3 to obtain 


'-1 1 Ol f Fj 1 0 ' 

0 -1 1 \f 2 • = ■ 0 - 

1 1/2 iJ[f 3 J [-1080 


(3.1.18) 


(3.1.19) 


or 





(3.1.20) 


NASA/TP— 2004-207430 


33 



Step 7 — Back-calculate the displacements, if required, from the deformation displacement relations: Displace- 
ments are calculated from any two of the three DDR as: 


*1 



= 0.024 


X 2 = P 3 = -0.024 


(3.1.21) 


The uniform thermal expansion induces a positive displacement (Xj) at node 1 and a negative displacement (X,) at 
node 2. 

Simultaneous solution . — The thermomechanical response for combined thermal and mechanical loads can be 
obtained by superposing mechanical and thermal solutions. An efficient alternative, however, is to combine both 
mechanical and thermal loads together in the calculation of forces. In this example, for clarity, mechanical and ther- 
mal loads are treated separately. The solution procedure illustrated for this simple example is quite general and 
applies equally well to other problems, irrespective of their complexity. 

Proration rule . — For mechanical loads only, when the force vector is multiplied by a constant factor, the displace- 
ments are prorated by the same factor, and vice versa. For thermal loads, however, this constant proration rule need not 
apply. Take, for example, the forces given by the thermal loads for cases 2 and 3, which differ by a factor of 3: 


h 

case 2 

K 

case 3 

'432' 


* = < 

h 

> = — < 

432 

kJ 


*3. 


432 


However, for the two thermal load cases, 



the displacements do not differ by a factor of 3: 

j 0.144 1 1 | fX] 1 case3 _ | 0.024 

[-0.144 jj * j|X 2 J “ [-0.024 


(3.1.22) 


(3.1.23) 


Illustrative Example 2: Propped Cantilevered Beam Under a Uniform Load 

A propped cantilevered beam of length l is subjected to a uniformly distributed load of intensity q per unit length 
as shown in figure 3.2(a). The beam, which is made of steel with a Young’s modulus E of 30 000 ksi, has a depth d 
of 2 in. and a width b of 1 in. Solve the problem for forces and displacements. 

Step 0 — Solution strategy: For this problem, conventional coordinate axes (x, y) and the origin at A are defined 
in figure 3.2(a). In addition, the abscissa (x )with its origin at B is defined because this choice reduces some calcula- 
tions. We solve the problem by considering three reactions as unknowns (n = 3). These are the moment M and shear 
force V at the fixed support and the shear force R at the hinged support (see fig. 3.2(b)). For the beam, two overall 
equilibrium equations — one rotational EE and one transverse EE — can be written (m = 2). The dual, or displacement, 
variables associated with the rotational and transverse EE are 0 and v, respectively. The problem is one-degree 
indeterminate because r = n-m- 1 . The displacement function w(x) is calculated because of the distributed nature 
of the load q. Since the beam is uniform, the solution is obtained in terms of El, where 1 is the beam’s moment of 
inertia. 
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Figure 3.2. — Propped beam under uniform load. 


Step 1— Formulate the equilibrium equations: Summation of the moment and shear force provides two EE in 
terms of three unknown forces (V, R , M). Rotational or moment equilibrium at the hinge point B yields the first EE. 


M-W + — = 0 
2 

The force equilibrium along the transverse direction is the second EE: 

V + R-qi = 0 

The two EE in terms of the three reactions (V, R, M) can be written in matrix notation as 


(3.2.1) 


(3.2.2) 


" e o -f 

V 


ql 1 ' 


R 

■ = ■ 

2 

-1 -1 0 

M 




(3.2.3) 


The equilibrium equations are one-degree indeterminate because three unknown reactions ( V,R,M ) are expressed in 
terms of two equations. One compatibility condition is required to solve the problem. 

Step 2 — Derive the deformation displacement relations: The DDR ({(3} =[B\ T {X}) have the following form: 



(3.2.4) 


In the DDR, the deformations ((3 V , (3 r , p m ) correspond to the reactions (V, M, R), respectively. The displacements 
(0, v) are the dual variables of the moment and transverse equilibrium equations, respectively. 
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Step 3 — Generate the compatibility conditions: The single compatibility condition is obtained by eliminating the 
two displacements (0, v) from the three DDR: 


[1 -1 


4 


Pv 

p r 

p m 


= { 0 } 


The correctness of the CC can be verified from its null property ([J3][C] 7 =[0]): 


i 0 
-1 -1 



0 

0 


(3.2.5) 


(3.2.6) 


Step 4— Formulate the force deformation relations: We can express the compatibility conditions ([C]{ (3} = 0) in 
terms of reactions by using the FDR, which can be obtained as the derivatives of the strain energy U as given by 
equation (2.40). 


Pv 

Pr 

P m 


du I 

r( 7ft dfft 

3V ~J 

o El dv 

3(7 | 

r - 7ft dm 

3 R ~ J 

0 El 3 R 

dll _ 

7ft dm 

3 M ~ ^ 

1 0 El dM 


(3.2.7) 


where the strain energy U for the beam can be written in terms of moment 7ft as (see eq. (2.39)) 


U = 


’ 7ft 
~2E1 


dx 


(3.2.8) 


For strain energy calculations, either point A or B (see fig. 3.2(b)) can be selected as the origin without any adverse 
consequence. When point A is selected as the origin, the moment (Tftix)) has to be written in terms of two unknowns 
(M and V). However, when B is selected as the origin with an axis x from B to A as shown in figure 3.2(b), the 
moment ( 7ft(x )) can be written in terms of a single reaction unknown (R). To reduce the number of calculations, 
we select B as the origin with axis x . The reader, however, could select the conventional axis x and arrive at the 
same results. 


7ft(x) = RI-^L- (3.2.9) 

The derivatives are 3 TflldV = 0, dM/dR = x , and 3 Tftl 3 M = 0. Thus, the deformation (3 V = p m = 0. The 
deformation p r can be calculated as 


Pr = 




( x)dx 


(3.2.10) 
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Upon integration, the deformation is obtained as 


Pr = 


El 


Rt 


sL 

8 


4 


(3.2.11) 


Step 5 — Express the compatibility conditions in terms of forces: In terms of force variables, the compatibility 
condition can be written as 


r 

-i 

V 

0 

0 

R 

3 El 

- 

M 



or 


[0 1 



3 qi 
8 


(3.2.12) 


Notice that the CC, which is an uncoupled equation, is nonhomogeneous because of the distributed nature of load q. 

Step 6— Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve for 
the forces: The compatibility condition can be coupled to the equilibrium equations to obtain three equations in terms 
of three unknown force variables: 


' i 

0 

-1 ' 

V' 


qi 2 ' 

2 

-1 

-1 

0 

R 


-qi 

0 

1 

0 

M 


3 qi 
8 


Solution of the equation yields the three reactions as 


V' 


5q£ 

8 

R 


3 qi 
8 

M 


qi 2 



8 


(3.2.13) 


(3.2.14) 


The moment function MCx ) given by equation (3.2.9) can be rewritten in terms of load q as 


M{x) = 


3 qix 



Step 7 — Back-calculate the displacements, if required, from the forces: Because the load is distributed along the 
span of the beam, the transverse displacement is a function of x ; that is, w(x). The moment curvature relation 
(k(x)) has to be integrated to determine the displacement function w(x). 
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Elementary Derivation of the Moment Curvature Relationship 

The moment curvature relationship derived in standard strength of materials textbooks is reviewed for com- 
pleteness. For a beam oriented along the x-axis and with the transverse displacement w(x), the curvature defined in 
elementary analytical geometry textbooks can be written as 


K = 


d~w 

dx 2 


1 + 


f dw ' N ~ 
dx 


-i3/2 


(3.2.15) 


where p = 1 /|k| is the radius of curvature. In simple beam analysis, it is assumed that the square of the rotation 
(i dwldx ) 2 is much smaller than unity; thus, the denominator is set to unity, and the curvature K simplifies to 


K = - 


d~w 
dx 2 


(3.2.16) 


Curvature is the deformation associated with the bending moment 7ft. Their product represents strain energy U. The 
strain energy due to flexure can be written as (see eq. (2.39)) 


or 



(3.2.17a) 


(3.2.17b) 


Since the moment function {7ft) is arbitrary and is not a null function, 


7ft d 2 w 

This simplistic derivation is sufficient because the process provides the moment curvature relationship, 
lem, the moment curvature relation can be written as 


(3.2.18) 
For the prob- 


K(» 


d 2 w(x) 
dx 2 


7ft 

f-1 

^3 qlx 

qx 2 \ 

YY' 

U/J 

l 8 

2 J 


Integrating the moment curvature relation yields the displacement function: 


(3.2.19) 


w(x) 


1 

Yl 


1 q£x 3 

16 



+ CjX + C2 


A 


) 


(3.2.20) 


The integration constants in the displacement function can be determined from the kinematic displacement 
boundary conditions, which are essential for the stability of the structure. The number of kinematic conditions N kbc 
can be calculated as the difference between the total number of displacement boundary conditions N tbc and the 
number of compatibility conditions r: 


^ kbc ^tbc r 


(3.2.21) 
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For this problem, N tbc = 3, r = 1, and N kh( . = 2. The three displacement boundary conditions are 


fiC-1 


w(x) = 0 

at 

x=0 


BC — 2 


o 

il 

IK 

T 

at 

X = l 


BC—3 


dw{x) _ 
dx 

at 

x = £ 

(3.2.22) 


Any two of the three displacement boundary conditions are sufficient for determining the two constants of integration 
(cj, c 2 ) in equation (3.2.20). The boundary conditions (SC-1) and (SC-2) are used to determine the constants, as 
follows: 


c 2 =0 


c\ =- 


si 1 

48 


The displacement function has the following explicit form: 




El 


( qlx 3 

16 


qx A 

24 


( 3 ^ 


48 


The reader can verify the compliance of the slope boundary condition (SC-3) at the fixed end. 
The maximum displacement can be determined from principles of calculus as 


0 . 0054 ^ 

El 


at x = 0.4215f 


(3.2.23) 


(3.2.24) 


(3.2.25) 


where w is considered to be positive along the y-axis. The load q, however, is along the negative axis, which is also 
the direction of displacement; or displacement is along the negative y-axis. 

Maximum slope, which occurs at the propped end at x = 0, is 


0 mflY 


48 El 


(3.2.26) 


The slope — = - — is positive along the positive z-axis (which also represents the negative z - axis). Likewise, 
\dx dx ) 

the maximum value of the bending moment is 


M 


max 


9 qi 2 
128 



(3.2.27) 
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Illustrative Example 3: Two-Span Beam Under a Uniform Load 


A two-span beam of length 2 1 is subjected to a uniformly distributed load of intensity q per unit length, as shown 
in figure 3.3(a). The beam, which is made of steel, has a Young’s modulus E of 30 000 ksi, a depth d of 2 in. and a 
width b of 1 in. Solve the problem for forces and displacements. 

Step 0 — Solution strategy : For the problem, the coordinate system (x, y) is defined with its origin at A (see fig. 
3.3(a)). Also, to reduce calculation, another abscissa (x) is defined with its origin at C, see figure 3.3(b). The prob- 
lem is solved by considering three reactions (R v /? 3 ) as the unknowns, or n = 3 (see fig. 3.3(b)). Since the beam 
is uniform, the solution is obtained in terms of El , where I is the moment of inertia. The three unknown reactions can 
be reduced to two because of symmetry (f? 3 = R^), which reduces to n = 2. Because the condition of symmetry is 
used, both rotational and transverse equilibrium yield the same single, independent EE, or m = 1. The single dis- 
placement unknown considered is v along the y-direction. The problem is one-degree indeterminate ( r = n-m= 1). 
The seven-step procedure is followed to solve the problem. 



cross section 


(a) Beam. 



(b) Unknown reactions. 



(c) Moment, 3£(x). 

Figure 3.3. — Two-span beam under uniform load. 
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Step 1 — Formulate the equilibrium equations: Since the condition of symmetry has already been used, only the 
transverse equilibrium is considered: 

2R l +R 2 ~2qi = 0 (3.3.1) 

The reader can verify that the moment equilibrium condition does not produce any additional independent EE. In 
matrix notation, the single EE can be written as 


[-2 -!]{*'} = (-2 q*} (33.2) 

The equilibrium equation is one-degree indeterminate because two unknown reactions (R v R 2 ) are expressed by a 
single EE. One compatibility condition is required to solve the problem. 

Step 2 — Derive the deformation displacement relations: The DDR ({(}} =\R\ T \X}) has the following form: 

Pl= — 2v 

p 2 =-v (3.3.3) 

In the DDR, the deformations (Pj, P 2 ) correspond to the reactions (R } , R 2 ), respectively. The displacement v along 
the y-direction is the dual variable associated with the transverse EE. 

Step 3 — Generate the compatibility conditions: The single compatibility condition is obtained by eliminating the 
single displacement v from the two DDR: 


The correctness of the CC can be verified from its null property ([B][C] r =0). 


[-2 - 1 ] 



= [ 0 ] 


(3.3.4) 


(3.3.5) 


Step 4 — Formulate the force deformation relations: The force deformation relations can be obtained as the 
derivatives of the strain energy U as 


Pi 


dU 

dR] 


P 2 


dU 

<)R 2 



dm 

dR 



(3.3.6) 


The strain energy function U for the beam has the following familiar form: 


C7 = 2 J 



dx 


(3.3.7) 


Because of the symmetry condition, the coefficient 2 and the limits of the integral (0 - 1) for a single span are used in 
equations (3.3.6) and (3.3.7). The strain energy stored in the structure is twice that of each span. For the calculation 
of the strain energy, we can select either point A or C as the origin without any consequence. The moment at any 
point x from origin A can be written as (see fig. 3.3(b)) 
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(3.3.8) 


. qx 2 

W(x) = R x x — 

The derivatives are = x and dM/dR 2 = 0. Thus, the deformations (p 2 = 0) and (p,) from equation (3.3.6) can 

be written as 


Pi 



or 


Pt = 


El 


*. 


8 


4^ 


(3.3.9) 


Step 5 — Express the compatibility conditions in terms of forces: The compatibility condition is obtained in terms 
of reactions by eliminating deformations between equations (3.3.4) and (3.39) as 


2 

i i ' 

C 

i 

pn l 

'ql A \ 

Yi 

L 3 J 

1*2 J 1 

4EI 


(3.3.10) 


Notice that the CC, which is an uncoupled equation, is nonhomogeneous because of the distributed nature of load q. 

Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the 1FM equations, and solve for 
the forces: The compatibility conditions can be added to the EE to obtain two equations in terms of two unknown 
reactions as follows: 


~ -2 

-f 

*t 



1/3 

0 

*2 


00 I'S, 


Solution of the equation yields the two reactions as 


(3.3.1 1) 


From symmetry. 


*1 


3qi 

8 


1*2 


R 3 =R 


5 ql 

. 4 . 

_ 3 qt 
1 ~ 8 


The moment function M(x) for the first span AB can be written as 


7K(x) = 


3 qlx qx 1 
~ 8 2 


The moment at the central support B is obtained for x-l as 


(3.3.12a) 


(3.3.12b) 


(3.3.13) 


M b =- 


qi 1 


8 


(3.3.14) 
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Step 7 — Back-calculate the displacements , if required, from the forces: Because the load is distributed along the 
span of the beam, the moment curvature relations (derived in Illustrative Example 2) can be integrated to obtain the 
displacement function w( x). For the first span AB, 


K(X): 


d~w(x) 7%. 


dx 1 El 

The displacement function w(x ) for the first span is obtained by integrating the moment curvatures relations as 

a 

w(x ) i | jqzx t 

2 


(3.3.15) 


dx 


Integration yields the displacement function as 


d 2 w(x) 1 
“ El 


^ 3q£x qx~^ 


8 


W{x): 


El 


f q()C qx 4 

16 24 


+ CjX +Co 


For the problem, the three displacement boundary conditions are 


(3.3.16) 


(3.3.17) 


BC-l => 

s 

II 

o 

at 

x = 0 


BC-2=> 

s 

II 

o 

at 

x = e 


BC- 3=> 

w(x) = 0 

at 

x = 2 i 

(3.3.18) 


For span AB, the boundary conditions BC-1 and BC-2 are used to determine the two constants as follows: 

c 2 = 0 


c l = — 


qj 1 

48 


The displacement function for span AB has the following explicit form: 


v(x)-- 


( qtx 2 ’ qx 4 


El 


( 3 x^ 


16 24 48 


(3.3.19) 


(3.3.20) 


Slope, 


dw(x) _ 1 ( 
dx El 


3 qlx qx 

6 


16 


3 ^ 
48 J 


Slope at x = 0, 


dw 

dx 


x = 0 


qf 3 

48 El 


Slope at x = i. 


dw 

dx 


(3.3.21) 
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Because of symmetry, the displacement function given by equation (3.3.19) can be used for the second span by using 
the alternative abscissa ( x ) as shown in figure 3.3(b). 


Illustrative Example 4: Continuous Beam With Support Settling 

A two-span, continuous beam that is made of steel with a Young’s modulus E of 30 000 ksi is subjected to a 
load P at the center of its second span as shown in figure 3.4a. The moment of inertia of the uniform beam is 
7 = 100 in. 4 Solve the problem for forces and displacements for the following two load cases: 

Load case 1 : Mechanical load P only 

Load case 2: Settling of the central support, A = 0.25 in. 

Load case 1: Solution for a mechanical load . — 

Step 0— Solution strategy: The coordinates (x, y) with origin at A are shown in figure 3.4(a). The beam is mod- 
eled into three spans. The free-body diagram of the beam is shown in figure 3.4(b). Each beam span has two moment 
unknowns; thus the structure has six unknown moments (M j , M y , . . . , A7^), or n = 6. The problem has five unknown 



2a 2a — 

(a) Support settling in A inches. 




€ e 


(c) Moments used to derive the force deformation relations. 
Figure 3.4. — Two-span beam with setting support. 
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displacements, or m = 5. These are the slope 0j at support A, slope 0 2 at support B, displacement v under load at C, 
rotation 0 3 at C, and rotation 0 4 at support D. Five equilibrium equations can be written along the five displacement 
directions. The problem is one-degree indeterminate because r=n-m = 1. 

Step 1 — Formulate the equilibrium equations: The five EE for the problem are written along with five displace- 
ment directions as follows: 


Along 0j 

= Xj at support A, 

Mj =0 

Along 0 2 

= X 2 at support B, 

-m 2 +m 3 =0 

Along Vj 

= X 3 under load at C, 

a 

Along 0 3 

= X 4 under load at C, 

1 

II 

0 

Along 0 4 

= X5 at support D, 

§ 

ON 

II 

O 


The five EE ([B]{F} = {/-’() can be written in matrix notation as 


(M 6 -M 5 ) 

— — -P = 0 

a 


'-1 

0 

0 

0 

0 

0 



0 

0 

1 

-1 

0 

0 

0 

m 2 


0 

0 

0 

Ha 

-1 la 

—1/a 

Ha 

m 3 

= ■ 

-p 

0 

0 

0 

1 

-1 

0 

m 4 


0 

0 

0 

0 

0 

0 

1 

m 5 

m 6 


0 


(3.4.1) 


(3.4.2) 


The five EE in equation (3.4.2) are not sufficient to determine the six moments; thus, one additional CC is required. 
Step 2 — Derive the deformation displacement relations: The DDR ( {(3 } = \B\ T {X}) yield 


Pl=-*1 


X, 


-Xa 


P 2 =*2 


X 3 

P 5 = X 4 

a 


P3 - ~X 2 + — 

a 


h=—+x 5 

a 


(3.4.3) 
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Step 3— Generate the compatibility condition: The single CC for the problem is obtained by eliminating the five 
displacements from the six DDR: 


202 + 203 +04+05 


or 


[0 2 2 1 1 0 ]] 


0l' 

02 

03 

04 


= 0 


05 

06 


The correctness of the compatibility condition can be verified from its null property ([B] [C] T = [0]). 


(3.4.4) 


-1 

0 

0 

0 

0 

0 ' 

"O' 


o' 

0 

1 

-1 

0 

0 

0 

L 


0 

0 

0 

1/a 

-1/a 

-1/a 

1/a 

L 

1 

1 

0 

= 

0 

0 

0 

0 
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-1 

0 
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0 

0 

0 

0 

0 

1 


0 


Step 4— Formulate the force deformation relations: The FDR for a beam with the end moments derived earlier in 
equation (2.48) are used. For end moments Af ■ and Mj, the FDR shown in figure 3.4(c) can be written as 


0 / = 

Py = 


6 El 


(2Mj + Mj ) 


— — (ilfj + 2M : ) 
6 £/'• * J> 


(3.4.6) 


For element AB, the span is 2a, the end moments are M x and M 2 , and the deformations are 0j and 0 2 . The FDR 
can be written as 


0 1 =^-(2M 1+ M 2 ) 

6EI 


02 =—(M 1 +2Mt) 
6 £7 


(3.4.7) 


Likewise, for element BC with span a, end moments M 3 and Af 4 , and deformations 0 3 and 0 4 , the FDR can be 
written as 


03 

04 


6EI 


6 El 


(2 M 3 +M 4 ) 


(M 3 +2M 4 ) 


(3.4.8) 
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For element CD with span a, end moments M 5 and M 6 , and deformations (3 5 and (3 6 , the FDR become 


P 5 =— (2M 5 +M 6 ) 


+2 M 6 ) (3.4.9 

6 El 

Step 5 — Express the compatibility conditions in terms of forces: The compatibility condition is expressed in terms 
of moments by eliminating deformations between the CC given by equation (2.18) and the DDR to obtain 

{AM X + 8 Mn + 5 + 2M5 + ) = 0 (3.4.10 

6EI 

Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve for 
the forces: The EE and CC are coupled to obtain the following IFM equation: 


"-1 0 0 

0 1 -1 


0 ] \M 1 
0 M 2 


0 0 1/a -1/a -1/a 1/a M, 


0 0 0 
0 0 0 


0 \\M a 


Because of the sparsity of the matrix, the equation is solved easily to obtain the moments: 


Mj = 0 M 2 


= 13Pa M _ o (3 412 ) 

3 32 


The reactions are obtained by writing the EE along the restrained directions at supports A, B, and D (see fig. 3.4(b)): 


At support A, 


M2ZM l = q 

2 a 


or Ra =- — 


(3.4.13a) 


At support B, 


n M x -M 2 M a - A/3 _ 

Ad — u or 

2 a a 


(3.4.13b) 


At support D , 


R D — —————— = 0 

a 


or Rn=- 


(3.4.13c) 


We can verify that the sum of the reactions ( R A + R B + R c = P) is equal to the applied load P. 
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Step 7 — Back-calculate the displacements, if required, from the forces: The deformation displacement relations are 
used to calculate the displacements. The displacement X 3 under load P is related to the deformations (P 4 , P 5 ) as 


*3=-f(P4+P 5 ) 


also, P 4 +P 5 

2 

- 6 £/ ( m 3 + 2 m 4 + 2 m 5 + m 6 ) - 96EI 

(3.4.14a) 

or 

x 3 =— 23 

3 192 El 

(3.4.14b) 

The displacement X 3 is along the negative ^ -direction. Likewise, other displacements are calculated as 


Rotation at support A, 

laM-y Pa 2 
X] ~ ~ 6EI 16£7 

(3.4.14c) 

Rotation at support B, 

x 2 =p, = 4flM2 = Pa ~ 
6EI 8 El 

(3.4.14d) 

Rotation under load at C, 

2 

X 4 =-(P 4 -P S ) = -- — 

4 2 VF4 64EJ 

(3.4. 14e) 

Rotation at support D, 

0 

X 5 =-(p 4 +p 5 +2p 6 ) = ^^ 
5 2 V' 4 Pi po/ 16EI 

(3.4.14f) 


Load case 2: Solution for support settling . — To simplify the calculation, we assume the mechanical load to be 
absent (P = 0) and solve the problem for the settling of supports only. Support settling is accounted for in the right 
side of the compatibility condition in the { 8 R } term. 


{5/?} = -[C]{p}° 


(3.4.15) 


where { (5 }° = ~{B r }X = {B r } A because X = -A (see eq. (2.14d)). 

The column vector {B r } is associated with the EE at support B along the direction of settling A, which for the 
problem is along reaction R B . The reaction R fi can be written in terms of moments (see eq. (3.4.13b)) as 


R B ={B r } T 


' 1/2 a 

T 


-1/2 a 


m 2 

-1/a 

- « 

m 3 

1 la 


M a 

0 


m 5 

0 


m 6 


(3.4.16) 
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{p} 0 =A/B r 7 = A 


(3.4.17a) 


' 1 / 2 a ' 

-l /2a 
-1/a 

!• 

1/a 
0 
0 

/5/?/ = -[C]{p}°=-[0 2 2 1 1 0]{p}°= — 

a 

The IFM governing equation can be rewritten by adding { S/? } and setting the mechanical load to P = 0 as 


-1 

0 

0 

0 

0 

0 




■ 0 ' 

0 

1 

-1 

0 

0 

0 


m 2 


0 

0 

0 

Ha 

-1/a 

-1/a 

1/a 


m 3 


0 

0 

0 

0 

1 

-1 

0 


m a 

. 

0 

0 

0 

0 

0 

0 

1 


M s 


0 

4 

8 

5 

4 

2 

1 


m 6 


12£7A 








2 

a 


The IFM equations are solved to obtain the moments: 

Afj =0 


M 2 


= m 3 = 


3EIA 
4 a 2 


M 4 


= m 5 = 


3EIA 
8 a 2 


M 6=0 


(3.4.17b) 


(3.4.18) 


(3.4.19) 
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Reactions are back-calculated as 


Ra 


m 2 

2 a 


3EIA 

8a 3 


R b = 


M 2 M 4 -M 3 _ 3 El A 

2 a a 4 a 3 


R C = 


M 5 

a 


3 El A 
8 a 3 


(3.4.20) 


In the absence of mechanical load (P = 0), the sum of the reactions (R A +R B + R^= 0) is zero, and symmetry about 
support B is maintained (R A = R c ). The settling A along the negative y-di recti on induces positive reactions at support 
B but negative reactions at supports A and C. 

Calculation of displacement . — To calculate the displacement for the settling of supports, we decompose the total 
deformation {(3} that gives rise to the displacement into an elastic component (|3} e and an initial component {(3|° as 

{(3} = {(3} e +{(3}° (3.4.21) 

The elastic deformation {(3} e is calculated from the moments and initial deformations, and {(3}° is determined 
from the support settlement given by equation (3.4.17a). The total deformation obtained as the sum of the two 
components is as follows: 


e 2 aM 2 A 

Pi — — 

6 El 4 a 


Pl=(Pl+Pl 

11 

0 1 > 

_ 2 aM 2 _ A 
P2 3EI 2 a 


P 2 =0 



P®=-- 

a 



K= ii? ("j +2 " 4) =|; 

P 4 =- 

a 



_ aM 5 _ A 
9 5 3EI 8 a 

p°=0 



_aM 5 _ A 
Pe 6EI 16 a 

p°=° 

is 

(3.4.22) 

Rotation at support A, 

X '^'-Ta 


(3.4.23a) 

Rotation at support B, 

x 2 = p 2 = 0 


(3.4.23b) 
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Rotation at C, 


(3.4.23c) 


Displacement at C, 
Rotation at support D, 


y _ P4 ~ P5 _ 9 A 

4 2 16a 

X 5 = “(P4 + P5 + 2 P 6 ) = 4 ^ 


(3.4.23d) 

(3.4.23e) 


The symmetrical support settlement yields symmetrical displacements. The rotation, or slope, is zero at the 
center support B, 0 2 = X, = 0; and the rotations are equal in magnitude, with opposite signs, at A and D, Q\= X { = 
% = ~Xy 


Illustrative Example 5: Propped Beam for a Mechanical Load, a Thermal Load, 

and Support Settling 

A uniform propped beam of length 2 a with a moment of inertia I of 100 in. 4 is made of steel with a Young’s 
modulus E of 30 000 ksi. It is subjected to a load P at the center span as shown in figure 3.5(a). Solve the problem 
for forces and displacements for the following load cases. 

Load case 1 : Mechanical load P only. 

Load case 2: Thermal load — Temperature is assumed to be uniform along the length of the beam. Along the 
depth, the temperature variation is linear, with values T u and T ( at the upper and lower surfaces, as shown in 
figure 3.5(b). 

Load case 3: Settling of the simple support A by the amount A. 

Load case 1: Solution for a mechanical load . — 

Step 0— Solution strategy: A coordinate system (x, y ) with its origin at A is shown in figure 3.5(a). A second 
coordinate system (x, y) with its origin at B (see fig. 3.5(a)) is also selected to reduce calculations. The beam is 
divided into two elements. Two moments are selected as the unknowns for each element. The beam has four moment 
unknowns (M v M 2 , A/ 3 , M 4 ); thus n = 4. The beam has three displacement unknowns — rotation at support A and 
rotation as well as transverse displacement at B,orm = 3. Three equilibrium equations can be written for the beam. 
The beam is one-degree indeterminate, or r = n-m = l. 

Step 1 — Formulate the equilibrium equations: Three equilibrium equations can be written for the problem (see 
fig. 3.5(c)). 

(1) Rotational EE at support A, Mj=0 (3.5.1a) 

(2) Rotational EE at B, ~M 2 + M 3 = 0 (3.5.1 b) 


(3) Transverse EE at B, 


m,-m 2 Ma-M 3 J „ 

—t L + ^ L + p =0 


We can simplify the three EE to a single equation by using the condition Mj = 0 and replacing Af 3 in favor of My. 



(3.5.2) 


The single EE is expressed in terms of two unknown moments. The problem is one-degree indeterminate, and one 
CC is required for its solution. 
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Enlarged beam cross section 


Temperature variation 


(b) Temperature distribution shown on enlarged beam cross section. 




M 2 —M2 

~a~ ~a~ 


(d) Bending moment for span AB. 



Mq — M2 M2 - M^ 


a a 

(e) Bending moment for span BC. 

Figure 3.5. — Propped beam with settling support. 
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Step 2 — Derive the deformation displacement relations: The two DDR ({ p } = are as follows: 


P 2 =- 


2v 

a 


P 4 =- (3-5.3) 

a 

The deformations (P,, P 4 ) correspond to the moments (M„ M 4 ). The transverse displacement v at B is the dual vari- 
able of the transverse EE at B, and it is considered to be positive along the load P, or negative ^-direction. 

Step 3 — Generate the compatibility condition: The single CC is obtained by eliminating the displacement v from 
the two DDR as 


P 2 + 2p 4 = 0 


„ r [1 2]{j£} = 0 

The null property of the CC and EE matrices ([/J][C] 7 = [0]) can be verified as 


= [ 0 ] 


[-2/a \ta]y 

Step 4— Formulate the force deformation relations: The FDR for the structure are obtained as 

i r 2 « drm 


1 c- c 

p-> = — f 

2 El Jo 


■m- 


dMo 


- dx 


(3.5.4) 


(3.5.5) 


l r 2 

Pa =— \ 

4 El Jo 


l r 2 « dm J 

Ut dx 

3M 4 


For span AB. 


Mo 

M{x)-^~x 

a 


d7X(x) x 


dMt(x) 


= 0 


(3.5.6) 


dMo a dM 4 

For span BC, the origin of the coordinate system is selected at B, and x is measured from B to C (see fig. 3.5(a)). 
M a -Mo _ 


(3.5.7) 


M{x) = - 


-x + Mo 


dfflt(x) _ x 
dMo a 


dtM{x) _ x 
i)M a a 


The deformation p 2 can be obtained by adding contributions from spans AB and BC as follows 


P2 El 1 Jo dMo J a dMo 


2fl _. dw 
2 


(3.5.8) 


(3.5.9) 


The first part of the integral for span AB becomes 


^ [ei 


Mo ) / jci aMo 

— — x — dx = — 

a J[a J 3EI 


(3.5.10) 
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We evaluate the second part of the integral for span BC by using the (x,y) coordinates as shown in figure 3.5(a). 
We change the integration limit from ( a to 2d) to (0 to a) and calculate the moment2£( x) using the origin at B: 


£/) j () ^ + }(' ' ~ 6E/ {2Ml + Ma ] 


or P2 =— (4 M 2 +M 4 ) 

Likewise, we can calculate the deformation p 4 by adding contributions from both spans: 




I j r« , r 2a „ dfft , 

— ( 7K- dx+\ 7K- dx 

I I I JO 3M 4 Ja 3M 4 




-M z - ]\fx\ _ 

-x + M 2 \ - dx 


:(M 2 +2M 4 ) 


Step 5 — Express the compatibility conditions in terms of forces: The CC in deformations (P 2 + 2P 4 = 0) can be 
expressed in terms of moments by eliminating deformations in favor of forces between the CC and the FDR. The CC 
expressed in terms of moments has the following explicit form: 

M 2 + 5 M 4 ) = 0 (3.5.14) 

Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve for 
the forces: The EE and the CC are coupled to obtain the following IFM equation: 

'-2/a 1 / o~| f M 2 I -P 

• = • [ (3.5.15) 

6 5 J[m 4 J 0 

The solution of the IFM equation yields the moments: 


Ma =-- 


Other moments are obtained by inspection: 


M, =0 


M 3 = M 2 = - 


(3.5.17) 
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Reactions ( R A , R c ) at the support are back-calculated as 


At support A, 



5 P 
16 


(3.15.18 a) 


At support C, 


R b = 


m 2 -m 4 

a 


11 P 
16 


(3.15.18 b) 


Step 7 — Back-calculate the displacements, if required, from the forces: Displacement v under load P is 

b= “ P4 = ^7 (2M4+M2)= "S 7 

Displacement v is verified by recalculating it from the DDR (p 2 - —2 via) as 


(3.5.19) 


2$ 2 


IPa 3 

96EI 


(3.5.20) 


The deformations (P t , P 3 ) are not required to determine moments; thus, those calculations were not necessary. 
However, these two deformations are required to determine the rotations at support A and at load application point B. 
The two deformations and the two rotations can be calculated from deformation moment relations and the DDR. The 
deformations (Pj, P 3 ) can be calculated from deformation moment formulas as 


For span AB , 


Pi - -~j { 2M l + m 2 ) ~ 


5 Pa 2 
96EI 


Likewise for span BC, 


P 3 = 


a 

6 El 


(2 M 3 +M 4 ) 


Pa 2 
24 El 


(3.5.21) 


P 4 =-^(2M 4 + M 3 ) = - 


IPa 1 
96 El 


(3.5.22) 


Calculation of rotation requires the deformation displacement relation, which, in turn, requires all three equilib- 
rium equations, which can be rewritten (see eq. (3.5.1)) as 


-1 0 

0 1 

1/a -1/a 


0 

-1 

-1/a 





' 0 


m 2 

* 

0 • 


m 3 


-p 


m 4 \ 




(3.5.23) 
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Rotation 9j, rotation 0 2 , and displacement v represent the dual variables of the EE (1, 2, 3, respectively) in equation 
(3.5.23). The deformation displacement relations ((Pi = lB\ r {X} ) can be written as 


Pi =-9i + — 

a 



a 


P3 --02 

a 

P 4 = V (3-5.24) 

a 


Note that the deformation P 2 calculated earlier is different from the deformation designated p 2 because the 
condition of symmetry is not used here. 


01 - P4 “Pi - - 


Pa- 

&E1 


02 - _ (P3 +P4) : 


Pa 1 

32EI 


(3.5.25) 


Load case 2: Solution for thermal loads . — We obtain the solution for a thermal load by modifying the IFM 
equations that have been developed for mechanical loads. The modification pertains to the term {S/?} in the right 
side of the compatibility condition. To calculate {8/?}, we decompose the total deformations into elastic and thermal 
components. 


P 2 =PI+P2 


P 4 =P1+P4 


(3.5.26) 


Since pf ar| d P 4 have already been calculated for mechanical load case 1, only the thermal deformations 
|p 2 and P4 j need to be determined. These are calculated from equation (2.40) as 


P 2 = 
P 4 = 



m t dm 

El dM 2 

M T dW 
El dM 4 


dx 


dx 


(3.5.27) 
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The thermal moment is calculated as 


rdf . 2 

Mr = 

i-dfl 

EaAT b y dy 

= Eab 

rd/2 

'-d/2 

AT ydy 

= £a& 

rd/2 . 

'T e +T u 

-drv 

, 2 


Eabd 

12 


{T u -T e ) 


El 


Jo 


ra dM(x ) f a dm(x) 

0 3m 2 




£7 Jl J 0 a Jo'' 


|<7x) 


(3.5.28) 


or 


? 2 =- 


uMj 


El 

m t r a 


p i= J?L [-£<** 

El Jo a 


(3.5.29a) 


P 4 : 


qMt 
2 El 


{ 8 /?} = -[C] {(3'} = -[1 2] 


u M 1 


El 
_ aMj 


2EI J 


2A /y’fl 
£/ 


(3.5.29b) 


(3.5.30) 


The IFM equations are obtained by adding { 8 /?} into the mechanical load equations and equating mechanical 
load P = 0 as 


2 

1 

m 2 


0 

a 

a 




6a 

5a 



2 Mja 

_~6EI 

6 E/. 

M 4 


EI 


(3.5.31) 
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Solution of the IFM equation yields 


(3.5.32) 


The moments (M j, Mf) are obtained as 


Mi =0 


i\i 3 — M'y — ■ 


Reactions can be back-calculated from moments as follows: 


M 2 _ 3 M t 

a 4 a 


M 2 -M 4 3 M t 


For thermal loads, the reactions are self-equilibrating ( R A + R c = 0) as expected. 

Calculation of displacements . — Displacements can be calculated from the deformation displacement relations by 
including thermal load contributions. Displacement v at B is given by 


: a(3 4 -^(Pl+Pl)-— — (2M 4 +M 2 )- 


a~ Mf 


For positive M T , or for T U >T ( , the displacement v at B is along the positive y-direction. 

Load case 3: Support settling . — Support settling has to be included in the right side of the compatibility condi- 
tion { Si? } . Since support A settles by A inches in the negative direction of the reaction R A (see fig. 3.5(a)), the EE has 
to be written in that direction as 


• = R} r {F}= 


1/aT \M 2 

oj k 


(3.5.36) 


rif) r 1 fA /a 

{P}°=A{ J B r }=A^ () = 


(3.5.37a) 


{5/?} = -[C]{P}°=-[l 2} 


0 a 


(3.5.37b) 
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The factor (6 Ella) earlier set to unity for mechanical loads has to be retained in {8ft } because CC is a 
nonhomogeneous equation. 



(3.5.38) 


The IFM equation, which includes the effect of support settling, but not mechanical or thermal loads, becomes 


or 


-2 la 

1 Id 

m 2 


0 ' 



m 4 

■ = ■ 

6 El A 

6 

5 _ 


{ a 2 


M 2 


3£7A 

8a 2 


M 4 


3EIA 

4a 2 


M 3 


= m 2 = 


3 El A 
8 a 2 


and M\ - 0 

Reactions can be back-calculated from moments as follows: 

Mj 3EIA 

D — — 


(3.5.39) 


(3.5.40) 


Rr — 


M 2 -M 4 _ 3E1A 
a 8a 3 


(3.5.41) 


In the absence of mechanical load, the reactions are self-equilibrating ( R A + R c - 0). 

We can calculate the transverse displacement v at location B by using the DDR, p 2 or P 4 . Here, the displacement 
is calculated from both p 7 and p 4 to illustrate the participation of initial deformations (P}°. From the transverse EE, 
which is the first equation in equation (3.5.39), the DDR can be written as 


or 


p 2 =p 2 +p 2 - — 

a 

a $ 2 

v = 

2 


Likewise, 


P 4 -Pi + p4 - — 

a 


or 


v = 0P4 
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Elastic deformations are 


Initial deformations are 


Total deformations are 



(3.5.42) 


also, v = a(34=- — A (3.5.43) 

The same displacement value (v = -5A/16) is obtained from deformation p 2 , which has an initial component, and 
deformation P 4 , which has no initial component. 


Illustrative Example 6: Fixed Beam Under a Uniform Load 

A uniform beam of length l is fixed at both ends as shown in figure 3.6(a). It is made of steel with a Young’s 
modulus E and a moment of inertia 7. The beam is subjected to a uniformly distributed load of intensity q per unit 
length. Solve the beam for forces and displacements at its center span. 

Step 0 — Solution strategy: The free-body diagram of the beam is shown in figure 3.6(b). Because of symmetry, 
the moment M and shear force V are equal at ends A and B. Furthermore, the value of shear force can be determined 
from the transverse equilibrium equation as 


2 V = qi or 



The moment EE is automatically satisfied as follows: 


2 


(3.6.1) 


(3.6.2) 
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V V 


(b) Free-body diagram. 



(c) Free-body diagram for half of the beam. 
Figure 3.6. — Fixed beam under uniform load. 


The fixed beam provides no free displacement at either support A or B. For the solution of the problem, any point 
in the span of the beam with two unknowns (a moment and a shear force) can be selected. The center span point C is 
selected because displacement is required at this point. 

A coordinate system (x, y) with its origin at A is selected. The problem has three unknowns (or n = 3), consisting 
of a fixed-end moment M, and M c and V c at C (as shown in fig. 3.6(c)). The problem has two displacements — the 
transverse displacement v and the rotation 0 at C (or m = 2). The problem is one-degree indeterminate {r-n-m- 1). 

Step 1 — Formulate the equilibrium equations: Two EE can be written at point C: 


Transverse EE, 





(3.6.3a) 


Moment EE, Mq -M — t- — — = 0 (3.6.3b ) 

4 8 

Since V c is a known quantity, the moment EE represents the only independent equation, which in matrix notation can 
be written as 



The problem is one-degree indeterminate because the single EE is expressed in terms of two moments 
(M, M c ). Solution of the problem requires one CC. 


(3.6.4) 
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Step 2 — Derive the deformation displacement relations: The DDR for the problem have the following form: 

Pi =G 


P 2 =-0 (3.6.5) 

where the deformations (Pj, P 2 ) correspond to the moments (M, M c ), respectively. Rotation 0 at C is the dual vari- 
able of the EE. 

Step 3 — Generate the compatibility condition: The single compatibility condition is obtained by eliminating the 
rotation 0 between the two DDR: 


Pi +p 2 -0 


(3.6.6) 


In matrix notation, the CC can be written as 


11 i:H o) 

The null property {\B]{C] T = [0]) of the EE and CC matrices can be verified as 


(3.6.7) 


[1 


-1] 


:[ 0 ] 


(3.6.8) 


Step 4 — Formulate the force deformation relations: The FDR for the problem can be obtained from energy con- 
siderations (see fig. 3.6(c)) as 


Pi 


1 2 
— M- — dx = — 

El J o dM El J 0 


7tt 


dm 

dM 


dx 


(3.6.9) 


Because of symmetry, only half the span is integrated when deformation pj is being calculated; then this value is 
doubled to obtain Pj, where 


/ s ,, of qx~ 

W(x) = M + — x — — — 
2 2 


(3.6.10a) 


d7%(x) 

~dM 


(3.6.10b) 


dmx) 

dM c 


(3.6.10c) 


The deformation that corresponds to the moment M c is P 2 = 0 since fM(x)l dM c = 0. The deformation pj that 
corresponds to the moment M can be calculated as 


Pl = 




(1) dx 


(3.6.1 1) 
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„ 2 Ml ql 4 

R, = — + - — 

o 0/1 


El \ 2 24 


Step 5 — Express the compatibility conditions in terms of forces: The CC can be written in terms of moments as 


Step 6— Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve for 
the forces: The IFM equation is 


1 in -T- 


7 °|K 


Solution of the IFM equation yields 


M = — 


M C = ^~ 
L 24 


Step 7 — Back-calculate the displacements, if required, from the forces: The displacement function w(x) is obtained 
by integrating the moment curvature relations: 

(3.6.16) 

dx 2 El 

qlx x 2 ql 2 qlx qx 2 

where MM = M+ 2~ q ~2 = ~n + ~2 T (3A17) 

Integration of the moment curvature relations yields 


1 ql 2 x 2 + qlx 3 qx 4 


El 24 12 24 


+ CjX + c 2 


The displacement boundary conditions, w = 0 at x = 0 and x = (, can be used to obtain the integration constraints 
as Cj = c 0 = 0. The displacement function w(x) can be written as 

w f x ) = lS- qt ~.* l+ gfa3 (3.6.15 

Ely 24 12 24 y 

The displacement at the center span is obtained for (x = It 2) as 


384£7 


The displacement at C is along the negative y-direction, which is also the direction of load q. 
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Illustrative Example 7: Clamped Beam for a Mechanical Load, a Thermal Load, 

and Support Settling 

A uniform beam of depth d, thickness b, and moment of inertia I is made of steel with a Young’s modulus of E 
and a coefficient of expansion of a per °F. It is clamped at both ends (A and B) as shown in figure 3.7(a). The beam 
is subjected to a concentrated load of magnitude P. Analyze the beam for the following load cases: 

Load case 1 : Transverse load P at the center span 

Load case 2: Uniform temperature along the length of the beam (Along the depth, the temperature variation is 

linear with values AT and - AT at the upper and lower surfaces, respectively, as shown in figure 3.7(b).) 

Load case 3: Settling of supports A and B by A a and A g inches, respectively 

Load case 1: Solution for a mechanical load . — 

Step 0 — Solution strategy: The coordinate system ( x , y) with its origin at A is shown in figure 3.7(a). The beam is 
divided into two elements (1, 2); and three nodes (A, B, C). Four moments (M v M 2 , M 3 , M 4 ), as shown in figure 
3.7(c), are considered as the force unknowns of the problem, or n = 4. The problem has two free displacements at 
location C, consisting of the transverse displacement v and the rotation 0, or m = 2. The problem is two-degrees 
indeterminate ( r = n — m = 2). 



P 



€ = 2 a 


(a) Clamped beam under concentrated load. 



(b) Temperature distribution shown on enlarged cross section. 


r 



Figure 3.7. — Clamped beam under a concentrated load and settling support. 
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Step 1 — Formulate the equilibrium equations: Two equilibrium equations can be written at C, consisting of one 
transverse EE and one rotational EE. The transverse EE at C along the displacement v direction (see fig. 3.7(c)) can 
be written as 


Mi - M 2 M 4 - M 3 
a a 

The rotational EE at C along the 0 displacement direction can be written as 

— M9 + M 3 = 0 

The two EE in matrix notation can be written as 




(3.7.1) 


(3.7.2) 


1/a -1/a 
0 1 


-1/a 

-1 


1/a 

0 


M x 

M 2 

M 3 ' 

M 4 



(3.7.3) 


Two EE are expressed in terms of four unknown moments (M v M v M 3 , M 4 ), or the problem is two-degrees indeter- 
minate. Two CC are required for its solution. 

Step 2 — Derive the deformation displacement relations: The DDR ({|3} - \BY {X\) are obtained as 


p 2 =--+0 

a 

P 3 =---e 

a 

P 4 =- (3-7.4) 


The deformations (P l5 P 2 , p 3 , P 4 ) correspond to the moments (M v M 2 , M 3 , M 4 ), respectively. Displacements v 
and 0 are the dual variables of the transverse and rotational EE, respectively. 

Step 3— Generate the compatibility conditions: The two CC for the problem are obtained by eliminating the two 
displacements from the four DDR as 


Pi — P 4 = 0 


Pi +p 2 +p 3 +p 4 -0 


(3.7.5) 


In matrix notation, the CC can be written as 
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(3.7.6) 


10 0 -1 
1111 


Pi 

P 2 

P 3 

P4 


The null property ([B][C] r =[0]) of the EE and CC matrices can be verified as 


1 /a -1/a -1/a 1/a 

0 1-10 


1 1 

0 1 

0 1 

-1 1 


0 0 
0 0 


(3.7.7) 


Step 4— Formulate the force deformation relations: The FDR for a beam element (derived earlier, see eq. (2.48c)) 
are used here. For beam element AC, the FDR are 


^=—(2M l+ M 2 ) 


P 2=— (2^2+Mi ) 
6EI 


(3.7.8) 


Likewise, the FDR for beam element CB are 


P 3 =-^-(2M 3 +M 4 ) 
o hi 


p 4 =— (2 m 4 +m 3 ) 
6 El 


(3.7.9) 


Step 5— Express the compatibility conditions in terms of forces: Elimination of deformations (P } between the CC 
and FDR yields the CC in moments: 


a 2 
6£7 3 


1 -1 
3 3 


-2 

3 


M x 

M 2 

m 3 

M a 



(3.7.10) 


Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the 1FM equations, and solve for 
the forces: The IFM equations are as follows: 


1/a 

-1/a 

-1/a 

1/a 


M x 


-P 

0 

1 

-1 

0 


m 2 


0 

2 

1 

-1 

-2 

« 

M 3 

• = * 

0 

3 

3 

3 

3 


m 4 


0 


(3.7.11) 


Solution of the IFM equation yields the moments as 
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(3.7.12) 


.. Pt ,, _P* 

Ml — Mi — — 

1 8 8 


u P* 

Ma = 


Reactions R A and R r can be back-calculated as (see fig. 3.7(c)). 

M 2 —M\ P 


M3-M4 _ P 


Step 7 — Back-calculate the displacements, if required, from the forces: The displacements can be calculated from 
the deformation displacement relations. The displacement at C under load P is 

a 1 Pi 3 

v = a p 1 =^—(2M l +M,) = — — — (3.7.14) 

F1 6 EI y 1 “ ' \92EI 

The displacement v is along the negative v-direction, which is also the direction of load P. The slope at C is 

6 = P) +P 2 + = +M 2 ) = 0 (3.7.15) 

Slope 0 at the center span of the beam is zero because of symmetry. 

Load case 2: Uniform temperature . — For thermal analysis, the right side of the compatibility conditions, {8/?}, 
which is a nontrivial vector, is calculated as follows: 


{8 R} =-[C]{p}° 




}dx O' = 1, 2, 3, and 4) 

o £/ KJ 


The thermal moment M T is 


rd/2 

1 T = I EaT(y)ybdy 
J—d/2 


(3.7.18) 
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since T(y) - 


2A Ty 


,, EaAT ? 2£ 

Mr = bd = 

1 6 d 


f bd n 
12 


aAr 


or 


Mr 


2EIaAT 


(3.7.19) 


For span AC, the moment function %f(x), see figure 3.7(c), can be written as 


/ \ Mr — A/i 

%(x) = A/ 1 + — -x 


<m f _x N 
dM\ v a , 


<s% x 
dMr a 


Upon integration, the thermal deformations are obtained as 


(3.7.20) 


Pi=-" r 


El J 


1 — — I dx- 


aMf 
2 El 


a AT a 




aMj _ a AT a 
2EI ~~ dT 


(3.7.21) 


The span CB has no contribution to the thermal deformations pj and pi, because 

dm dm _ 

dMi dM 2 

Likewise, the thermal deformations p^and p^can be calculated as 

aMj _ ocA Ta 
2EI _ d 

aMj aATa 


P3=- 


P4=' 


2 El 


The effective initial deformation vector { 5 R ) becomes 


{«} = - 


10 0-1 
111 1 


aATa' 
d , 


4ocA7a fO 
d |l 


(3.7.22) 


(3.7.23) 
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The IFM governing equation for thermal load can be written as 


l/a 

-1/a 

—1/a 

1/a" 


Mi 


0 

0 

1 

-1 

0 


m 2 


0 

2 

1 

-1 

-2 


Mi 


0 

3 

3 

3 

3 _ 


m 4 


24£/aAT 
[ d 


Solution of the IFM equation yields the moments: 


Mi 


1 

M 2 

2EIaAT 

1 

m 3 

d 

1 

m a 


1 


(3.7.24) 


(3.7.25) 


and displacements can be calculated as follows: 





Pt - Pi +P5 - “ (2 M x + M 2 ) = 0 

6 El a 

(3.7.26a) 

Likewise, 

p 2 =o 


thus, 

v = 0 


and 

e = Pi +P 2 =o 

(3.7.26b) 


The temperature distribution, which varies along the depth of the beam but is constant across its length, does not 
induce displacements in the fixed beam. In other words, because the elastic and thermal deformations are equal and 
opposite, they cancel each other. The beam is stressed because of nontrivial elastic deformations. 

Load case 3: Support settling . — Support settling is also accounted for in { bR } , which is the right side of the CC. 
As before, this vector is calculated as 


{5/?} = -[C]{p}° 

{p}°=-[B r ]{x} 

l*}=-{^} <3 ' X27) 

Here, A a and A D are the settling of supports A and B along the negative y-direction, respectively, as shown in figure 
3.7(a). The reactions R A and R B corresponding to the amounts of settling can be expressed in terms of moments (see 
fig. 3.7(c)) as 
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or 


Ra=\ 


— 1/a ' 

T 

M x 


' 0 ] 

r 

M x 

1/a 

► - 

m 2 

• and R B = < 

0 

• * 

m 2 

0 


m 3 

1/a 


m 3 

0 


m a 


—1/a 


m a 



r 

-1/a 0 1 

f 

-A 

a] 


{p}°=-M*}= 




0 -1/a 


l-A, J 


{5/?} = -[C]{p}° 


f" A * 1 

'10 0-1 
111 1 

[~ a b J 



{§*} = 



0 


(3.7.28) 


(3.7.29) 


(3.7.30) 


The IFM equation for support settling can be written as 


1/a 

—\la 

-1/a 

1/a 


M x 


' 0 

0 

1 

-1 

0 

< 

Mo 

6 El 

0 

2 

1 

-1 

-2 


m 3 


a a ~ a b 

3 

3 

3 

3 _ 


M 4 


0 


Solution of the IFM equation yields the moments 


M x 


a a ~ a b 

Mo 

_ 6 El 

0 

m 3 

t 1 

0 

m a 


a b ~ a a. 


(3.7.31) 


(3.7.32) 


For an equal amount of settlement (A a = A B ), M x = M 2 = M 3 = M 4 = 0, and the structure is stress-free. Reactions can 
be back-calculated from the moments as 


R 4 =- 


Mo -M 


Y2EI (Ag - A a ) 

f 3 


Rn =- 


\2El{A B -A A ) 


(3.7.33) 


The reactions are self-equilibrating (R A + R B = 0) when supports settle and there is no other load. 
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Calculation of displacements. — Displacements can be calculated from the deformation displacement relations. 
The midspan displacement v can be written as 


- = <*Pl={(Pl+Pl) 


Pf = — (2 M l +M 2 ) = 
1 6EI v 


&A ~ A B 



{&A +A S ) 

v = 

2 


(3.7.34) 


The displacement v is along the negative y-direction (or along load P). 
The midspan rotation 0 can be written as 


or 


0 - Pi +02 


Pl=Pf+Pl 


(A a +A b ) 

l 


P2 = Pf + P ? 2 

P5=^( A ^-A b ) and P§=^- 


_ 3(A /t -A g ) 

It 


(3.7.35) 


For symmetrical settling of supports (A a = A fi ), rotation, or slope, 0-0. 


Illustrative Example 8: Torsion of a Shaft Fixed at Both Ends 

A circular shaft fixed at both ends and made of two materials with shear moduli Gj and G 2 and polar moments of 
inertia J x and J 2 , is shown in figure 3.8(a). The shaft, which has a total length of t, is subjected to a torque T at a 
distance a from its left support. Analyze the shaft for torque and angle of twist. 

Step 0— Solution strategy: Figure 3.8(a) shows the coordinate system (x, y) with its origin at A. The fixed-end 
reactive torques (T A . T c ) are considered as the unknowns, or n = 2. Only one EE, or m - 1 , representing torque bal- 
ance along the ^-axis, can be written. A dual variable, or angle of twist <p, is associated with the EE. The problem is 
one-degree indeterminate (r-n-m = 1). 
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y 

r ft 

T r 

X 

— ► 

1 A 4s 1 


y Gi ,j a \ 

/ G 2’ J 2 L 

- 1 

1 


(a) Composite shaft. 



t A t A ^ w T w ^ T C T C 

t a t c 

(b) Free-body diagram. 




(c) Solid shaft AB, shown (d) Hollow, tubular 

enlarged. shaft BC. 

Figure 3.8. — Torsion of a circular shaft fixed at both ends. 


Step 1 — Formulate the equilibrium equations: Torque balance yields a single EE for the problem: 

T A + T c ~T = 0 (3.8.1) 

In matrix notation, the EE can be written as 


[-1 - 1 } 


l Tc 


{ ~ T } 


(3.8.2) 


The problem is one-degree indeterminate because two unknown torques are expressed by one EE. 

Step 2— Derive the deformation displacement relations: The DDR ({(3} = [fi] r {X}) are obtained from the EE as 

Pi =-<P 

P 2 =~<f> (3.8.3) 

where cp is the dual variable of the EE and the deformations (j3j, P 2 ) correspond to the torques (T A , T c ), respectively. 
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Step 3— Generate the compatibility conditions: The single CC for the problem is obtained by eliminating (p from 
the two DDR: 


Pi -P 2 =o 

or [1 ~l]{j|;} = {0} 

The correctness of the EE and CC matrices can be verified from their null property ([B][C] r = [0]) as 

“ 1 


(3.8.4) 


[-1 -1] 


-1 


= [ 0 ] 


(3.8.5) 


Step 4 — Formulate the force deformation relations: The FDR for the problem can be determined from the follow- 
ing integral derived earlier in equation (2.45a): 


p = — f -L® 1 

JG Jo JG 

For span AB ( J = J 1 and G = G x ), the torque is 

T = —T a ant 

For span BC ( J = J 2 and G = G 2 ), the torque is 


dx 


o JG dT 


= - 1 
dT A 

= o 
dT c 


(3.8.6) 


(3.8.7) 


T = —T a + T and 


d T 


dT A 


-1 


dT 

dT r 


or 


^=T^\l- TA{ - l)dx+ TF~\ {- T A+T)(-i)dx 

J\(j\ ^2^2 J ° 

T A a { t A -T){l-a) 


Pt = 


J X G X J 2 G 2 
dT 


07 = 0 because = 0 

r z j r T < 


dT, 


(3.8.8) 


(3.8.9a) 


(3.8.9b) 


C 


Step 5 — Express the compatibility conditions in terms of forces: The CC written in terms of torque has the follow- 
ing form: 


a l-a 

J\G\ h G 2 



\ T(e-a) 
1 J 2 G 2 


(3.8.10) 
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Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the IFM equations , and solve for 
the forces: The IFM equation can be written as 


-1 

-l 

T a 


-T 

a l — a 

o 


> = « 

f(e-a) 

J\ G l J 2 G 2 


T C, 


j 2 g 2 


(3.8.11) 


Solution of the IFM equation yields 


J 2 G 2 


a t- a 
- + - 


V. J \ G \ J 2 G 2 J 


■ = T\ 


(3.8.12) 


l-a 

J 2 G 2 


a a 

- + - 


V J\G\ J 2 G 2 ) 


For l = 2a, J x -J 2 = J, and G, = G 9 = G, the torques T A and T c simplify to 



The total deformation tp is obtained by substituting T = 2 T c and (T A = T c ) in equation (3.8.9a) as 


<P = -Pi =~ 


T A a 

T C 

Ug, 

J 2 G 2 



(3.8.13) 


(3.8.14) 


For a uniform shaft, J l =J 1 = J, G, = G, = G, and external torque T applied at the center of the shaft, the deforma- 
tion becomes 


Ta 

2JG 


Ta 

"f 

2JG 


= 0 


The twist angle tp is 0 because of symmetry. The twist angles for the individual shafts ( AB and BC) are 


(3.8.15) 


<?AB ~ - 


T A a 

J\ G \ 


$BC ~ 


T C (t-a) 

J 2 G 2 


(3.8.16) 


Next consider a shaft made of two different materials. Shaft AB is a 4-in.-diameter solid shaft (fig. 3.8(c)) made 
of annealed bronze. It has a shear modulus Gj of 6500 ksi and a length a of 6.5 ft. Shaft BC is a hollow, tubular shaft 
with an outside diameter of 4 in. and an inside diameter of 2 in. (fig. 3.8(d)). It is made of aluminum with a shear 
modulus of 4000 ksi and a length b of l- a = 5 ft. The applied torque T =20 000 in.-lb. Numerical values for the 
composite shaft are 
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a - 78 in. 


7, = = 8it in. 4 


4 . 2 

7jGj =5.271x10 kips -in. 




in. 


J 2 G 2 


= 37txl0 4 kips- in. ^ 


£ - a = 60 in. 


For these values, 



Jl 1.4286' 
{ 8.5714 


10 3 in.-lb 


Twist angle for the first shaft, 
Twist angle for the second shaft, 
Total twist angle, 


q> AB = -54.56 xlO -4 rad 
(p BC = 54.56 x 10 ^ rad 
9 = <Pab +( Pbc =0 


£ = 138in. (3.8.17) 


(3.8.18) 


Illustrative Example 9: Beam Supported by a Tie Rod 

A steel beam of length L, modulus of elasticity E b , and moment of inertia l b is fixed at C and supported by a tie 
rod at B. It is subjected to a uniformly distributed load of intensity q per unit length. The tie rod, which is made of 
aluminum, has a cross-sectional area of A f , a modulus of elasticity of E t , and length l. Determine the force in the tie 
rod, the reactions at C, and the displacement at B. 

Step 0— Solution strategy: A coordinate system (x, y) with its origin at B is shown in figure 3.9(a). Tensile bar 
forces are assumed to be positive. A reaction (R B ) at B and two reactions (Rq 37 at C are considered as the three 
force unknowns of the problem, or n = 3. From the ffee-body diagram of the problem shown in figure 3.9.(b), the 
internal force F in the tie rod is equal to the reaction at B (or F = R b ). Two equilibrium equations— one along the 
transverse y-direction with v as the dual displacement variable and the other along the rotation 0 or the ^-direction 
can be written, or m = 2. The problem is one-degree indeterminate ( r=n — m= 1). 

Step 1— Formulate the equilibrium equations: The transverse and rotational EE of the problem are as follows: 

Transverse EE: Rb + -qL = Q (3.9.1a) 

2 

Rotational EE at B: /V/ c + /? c Z. - = 0 (3.9.1b) 

In matrix notation, the EE can be written as 




r b 



-1 -1 o " 




-qL 



Rc 

■ = ■ 

_^L 

0 -L - 1 


M c, 


2 . 


Since three force unknowns ( R B , R c , M c ) are expressed in terms of two EE, one CC is required for the determination 
of the unknown reactions. 
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(b) Free-body diagram. 

Figure 3.9.- Beam supported by a tie rod. 


Step 2 — Derive the deformation displacement relations: The DDR ({ P } = [B\ r {X}') for the problem are as follows: 

Pi =-v 

P 2 = — v- LB 

P 3 =-0 (3.9.3) 

Deformations (P } , P 7 , P 3 ) correspond to the forces (R B , R c , M c ), respectively. Displacement v and rotation 0 are the 
dual variables of the transverse and rotational EE, respectively. 
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Step 3— Generate the compatibility condition: The single CC is obtained by eliminating v and 0 from the three 
DDR as follows: 


In matrix notation, the CC can be written as 


Pi — P 2 + £$3 - 0 


(3.9.4) 


[1 


-1 


4 


Pi' 

P2 f 


IP 3 


= { 0 } 


The null property ([B][C] r = [0]) of the equilibrium and compatibility matrices can be verified as 


(3.9.5) 





1 



-1 

-1 

0 


-1 


O' 

0 

-L 

-1 




0 





L 




(3.9.6) 


Step 4— Formulate the force deformation relation: The deformation in the CC represents the total deformation, 
which is composed of two components: beam deformations {p} beam and tie-rod deformations {P) uerod , as 

{P} = {P} beam + {P} ne rod (3.9.7) 


Calculation of the beam deformation contributions . — The beam deformation component pj* 3 " 1 due to the 
reaction R B is obtained as the partial derivative of the strain energy U b stored in the beam with respect to R B : 


P beam 
1 


du b 

dR B 


1 r l - 
— 

El Jo 


dm 

dRi> 


dx 


(3.9.8) 


where 


m{x)=R B x 

dm(x) 

= x 

dR B 


(3.9.9) 


or 


Pi 


beam 


f RbL 3 


E b h 



(3.9.10) 


Deformation components (p!f aiT1 = plf 3 " 1 = 0 ) are zero because dm/dR c = dWdM c =0. 

Calculation of the tie-rod deformation contributions . — The deformation component p[ lerod is calculated as 
the partial derivative of the strain energy stored in the tie rod with respect to the internal force (F = R B ) as shown in 
figure 3.9(b). 

pile rod = fl£_ = _L_ [V — dx 

n 3 R B A,E, Jo 3 R b 


since 


o tie rod _ E B l 

Pl “ A t E, 


(3.9.11) 
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Deformations '« tierod ° tierod 
Total deformations are 


dF dF 


= P 3 Crou = 0) because — = — = 0. 


3 Rq 3 Mq 


Pi = 


Ebh 


R B lf qL 


4 'i 


8 


Rr 


A,E, 


P 2 =p 3 =0 (3.9.12) 

Step 5 — Express the compatibility conditions in terms of forces: The CC in forces has the following form: 


^E b I b A t E t 


0 0 


-1 

Rb 

< 

R c 

- 

Me. 


qL* 

*E b I b 


(3.9.13) 


Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve 
for forces: The IFM equation for the problem has the following form: 


-1 


0 


t 


-1 

—L 

0 


jSEf, I b A t E t 

The solution of the IFM equation yields the forces as 

Rr =■ 


0 

-1 

0 




-qL 

r b 

■ R c 

> = < 

qL 2 

2 

M c . 


qL 4 



8 E b l b 


(3.9.14) 


qL * 


*E b I b 


' L 3 


3 E b Ib A,E j 


(3.9.15a) 


t) 



(3.9.15b) 


M c 


qL 1 


*E b l b 


qL 3 
3E b I b 


A,E t ) 


(3.9.15c) 


The forces are calculated for the following numerical values: E b = 200 Gpa, E t = 70 Gpa, L = 3xl0 3 mm, 
i = 7.5xl0 3 mm, q = 12 kN/rn, I b = 20xl0 6 mm 4 , A t = 100 mm 2 . 
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7? s =9.145kN 


R c = 26.855 kN 


M c = - 26.565 kN-m 


Step 7 — Back-calculate the displacement, if required, from the deformation displacement relations: The trans- 
verse displacement (v) at hinge point B can be calculated as 


v = -Pi = 


( R b l ■ 


E bh 


qL‘ 


4 ^ 


Rb* 

A,E t 


(3.9.16) 


The first term in equation (3.9.16) represents the contribution from the beam, whereas the second term represents the 
contribution from the tie rod. 

The axial displacement v tie rod for the tie rod can be determined from its deformation 


tie rod 
1 


R B t 

AE 


or 


tie rod _ _ ^ 

A t E t 


(3.9.17) 


Substitution of the numerical values yields the displacement v tie rod = -9.798 mm, which represents the stretching of 
the tie rod. 

The transverse displacement for the beam can be calculated from its deformation 


pbeam 


1 

E b I b 


3 



or 

For the numerical values of the 


.beam 




qt 


4 't 


E b h 


(3.9.18) 


problem, the transverse beam displacement along the negative y-direction becomes 


v beam = 9 798 ^ (3.9.19) 

The beam and tie rod deform in a consistent manner (v tie rod + v beam = 0), as expected. 


Illustrative Example 10: Three-Bar Truss for Mechanical and Thermal Loads 

A three-bar truss made of steel has a Young’s modulus E of 30 000 ksi and a coefficient of thermal expansion a 
of 6.6X10 -6 per °F (fig. 3.10). The areas of its three bars (A u A 2 , A 3 ) are 1.0, 1.0, and 2.0 in. 2 , respectively. Analyze 
the truss for the following two load conditions: 

Load case 1 : Mechanical loads (P x = 50 kips and P y = 100 kips) shown in figure 3.10 
Load case 2: Two cases of temperature variations 




100.0] 

°F 

ATi' 


100.0' 

A T 2 

- = * 

200.0 

• and 

A T 2 

. = - 

-200.0 

at 3 

easel 

300.0 J 


A7 3. 

case 2 

-300.0 
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Here A 7} represents variation in temperature for bar i. 

Step 0 — Solution strategy: A coordinate system (x, y) with its origin at node 1 is shown in figure 3.10(a). The 
bar forces (Fp F 2 , Ff) are the three force unknowns, or n = 3. Only two equilibrium equations can be written along 
the two displacement directions (X ] , X 2 ) at node 1 , or m = 2. The truss is one-degree indeterminate (r = n-m = 1). 

Load case 1: Solution for mechanical loads . — 

Step 1 — Formulate the equilibrium equations: The two EE ([5]{F}) = {P}) of the problem are obtained from 
the force balance condition at the free node 1 along displacements X x and X 2 , as shown in figure 3.10(b). 


1/V2 0 

-1/V2 -1 


-1/V2 

-1/J2 


Fi 

F 2 ) = 




(3.10.2) 


The truss is one-degree indeterminate because the two EE have three force unknowns. One CC is required for the 
analysis of the three-bar truss. 

Step 2 — Derive the deformation displacement relations: The DDR ({(3} = [B] T {X]) of the truss has the 
following form: 


_ X 

“72 


*2 

V2 


p 2 -~X 2 

„ Xi X, 

(3.10.3, 



(b) Free-body diagram at node 1 . 


Figure 3.10.- Analysis of a three-bar truss. 
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Here, pj, p 2 , and P 3 are the bar deformations corresponding to forces F j, F 2 , and F 3 respectively. The nodal dis- 
placements (Xj, X 2 ) are also the dual variables of the EE. 

Step 3— Generate the compatibility condition: The single CC for the problem is obtained by eliminating two 
displacements from the three DDR: 

p! — V2p 2 +p 3 =0 (3.10.4) 

The CC can be written in matrix notation as 


[l -VI l] 


Pi 

P 2 

P 3 


The null property ([B][C] r = [0]) of the EE and CC matrices can be verified as 


(3.10.5) 


l/VI 0 -1 /VI 

1 

-VI 


o' 

-1/ VI -1 -l/VI 

1 


0 


(3.10.6) 


Step 4— Formulate the force deformation relations: The FDR for the bars of the truss can be obtained as 
(P = IFIAE). The lengths of the three bars are 100 VI, 100, and 100VI in., and their areas are 1.0, 1.0, and 2.0 in. 2 , 
respectively. 


Pi 


A X E X 


100VI 


P 2 



P 3 



(3.10.7) 


Step 5 — Express the compatibility conditions in terms of forces: The CC are obtained in terms of forces by 
eliminating deformations between the CC and FDR: 


100 VI 

E 


[1 -1 1 / 2 ] 


>1 

F 2 

*3 



(3.10.8) 


Step 6— Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve 
for forces: 


l/VI 

0 

-1/VI " 

Fl 


50' 

-l/VI 

-1 

-1/VI 

F2 

• = • 

100 

1 

-1 

1/2 

Fy 


0 
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Solution of the IFM equation yields the forces as 


>1 


5.025 

^2 


42.893 

A 


75.736 


(3.10.10) 


Step 7 — Back-calculate the displacement , if required, from the defamation displacement relations: 


* 2= _ P , = _i^ = 0,43 in. 

E 

1 00 

X l =V2p 1 -p 2 = (2F 1 -F 2 ) = 0.110in. (3.10.11) 

E 

Load case 2: Solution for thermal loads . — Thermal analysis requires the inclusion of nontrivial {SR} in the 
right side of the CC: 


where 


or 


{8f?} = -[C]{f}}° 


{P}° = 


P? 


A7i ei 

-pO 

• = a 

AT 2 £ 2 



at 3 £ 3 

-V2 l]a- 

AT\t\ 
AT 2 £2 ’ 

AT 3 £ 3 


{8tf} = -lOO-Jla (A T y - A T 2 + A T 3 ) 


(3.10.12) 


M easel =-°- 187 (3.10.13a) 

Likewise, the { 8R } calculated for the temperature increase for case 2 becomes 

Mcase2=° (3.10.13b) 

The nontrivial thermal distribution for case 2 represents a compatible temperature distribution that does not induce 
any stress. For thermal distribution case 1, the IFM equation can be rewritten to include the {8/?} term as follows: 
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" i/VI 

o -l /VI 

F\ 


0 

-i /VI 

-l -l /VI 


► — . 

0 

1 

-1 1/2 

VI 


-39.6 


(3.10.14) 


The compatibility condition in equation (3.10.14) has been scaled with respect to 100 VI / E, see equation (3.10.8). 
Solution of the IFM equation yields forces as 


>1 


-13.59' 

Fl 

> = - 

19.22 

F 3 


-13.59 


Calculation of displacement . — Displacement can be calculated from the DDR as 

X 1 =V2|3 1 -p 2 


Pi=P^+Pi and P 2 =P!+P2 


upon substitution, 


X\ =-0.155 in. 


(3.10.15) 


(3.10.16) 


X 2 = -p 2 = -(p? + p' 2 ) = -0.196 in. (3.10.17) 

The compatible temperature distribution thermal load case 2, produces trivial forces (Fj = F 2 = F 3 = 0) but nonzero 
displacements, which are also calculated from the DDR: 


{p}={pr+{pr={pr = 


0.0933 
-0.1320 • 
-0.2800 


X { = VI Pi - P 2 = 0.264 in. 


Y 2 =p^ =0.132 in. 


(3.10.18) 


(3.10.19) 


Illustrative Example 11: Six-Bar Truss for Mechanical and Thermal Loads 

The six-bar truss shown in figure 3.1 1 is made of aluminum with a Young’s modulus E of 10xl0 3 ksi and a 
coefficient of thermal expansion a of 6.0x1 0 -6 per °F. The cross-sectional areas of members 1, 3, 5, and 6 are 1 in. 2 , 
and those of members 2 and 4 are 1 / VI in. 2 Analyze the truss for forces and displacements for the following two 
load cases: 
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Figure 3.11.- Six-bar truss. 


Load case 1 : Mechanical load P = 1 kip at node 1 along the y-direction 

Load case 2: Temperature increase AT = 100 °F (for member 3 only) 

Step 0 — Solution strategy: The coordinate system (x, y) with its origin at node 4 is shown in figure 3.11. The six 
forces in the six bars are the force unknowns, orn = 6. Tensile bar forces are assumed positive. Nodes 2 and 3, which 
are free, with displacements X v X 0 , X 3 , and X 4 , as shown in figure 3.1 1, yield four EE, or m = 4. The truss is two- 
degrees indeterminate (r = n-m = 2). From observation, we can conclude that the force in the sixth bar F 6 = 0, 
which, a priori, will not be assumed. 

Load case 1: Solution for mechanical loads . — 

Step 1 — Formulate the equilibrium equations: Forces in the six truss members (1, 2, . . . , 6) are designated 
F j, F 2 , . . . , F 6 . Four EE at nodes 1 and 2 can be written along the displacement directions (X l ,X 2 , X^, X 4 ). In 
matrix notation, the EE becomes 


1 1/V2 0 

0 l/y/2 1 

0 0 0 

0 0-1 


0 0 0 

0 0 0 

1 / V2 1 0 

-1 / V2 0 0 



F\ 



Fl 


0 ' 


F} 


1000 

• 

Fa 

• — < 

0 


f 5 


0 


F f> 



(3.11.1) 


The sixth column, which is null, in the EE matrix corresponds to the fully restrained force in the sixth bar. Since the 
four EE are expressed in terms of six forces; two CC are required for the solution. 

Step 2 — Derive the deformation displacement relations: The six DDR ({(3) = ) have the following form: 


Pl=*l 


(*3-*4) 

V2 


P2 


{X l+ X 2 ) 

V2 


P 5 - x i 


P 3 =*2-*4 


P 6 =0 


(3.1 1.2) 
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Step 3— Generate the compatibility conditions: The two CC, which are obtained by eliminating the four dis- 
placements from the six DDR, can be written in matrix notation as 

Pi 
P 2 
P 3 
P 4 
P 5 
P 6 

The first CC specifies that fl 6 = 0 is zero because member 6 is restrained at both ends and cannot deform. 

The null property ([B][C] r = [0]) can be verified from the one EE and the CC matrices as 


0 0 0 0 0 r 

1 -yj2 1 -V2 1 0_ 



1 1/V2 0 0 

0 1/V2 1 0 

0 0 0 1/V2 

00-1 -1/V2 


'0 1 

0 0] 0 - 4 i ro o' 

000 1 0 0 

1 0 0 -V2 “ 0 0 

000 1 [00 

|_1 0 


(3.11.4) 


Step 4 — Formulate the force deformation relations: The FDR p = FttAE for the six truss members are as 
follows: 


£ X F j _ 20 Fj 
A X E X ~ E 







(3.11.5) 


Step 5— Express the compatibility conditions in terms of forces: The CC is expressed in forces by eliminating 
deformations between the CC and FDR to obtain 


h 

20 [0 0 0 0 0 ll F 3 _J0) 

T|_i -2V2 1 1 oj ' f 4 [oj 

F 6. 


(3.1 1.6) 
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Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, 
for forces: 


1 

1/V2 

0 

0 

0 

o' 


K 


0 

0 

1/V2 

1 

0 

0 

0 


F 2 


1000 

0 

0 

0 

1/V2 

1 

0 

> 

F3 

> = < 

0 

0 

0 

-1 

-I a/2 

0 

0 


Fa 


0 

0 

0 

0 

0 

0 

1 


f 5 


0 

1 

-2J2 

1 

-2V2 

1 

0 


F 6 . 


0 


Solution of the IFM equation yields the forces: 


Ff 


-545.5 

F 2 


771.4 

f 2 


454.5 

Fa 


-642.8 

f 5 


454.5 

f 6. 


0.0 


Step 7 — Back-calculate the displacement, if required, from the deformation displacement relations: 

20 Fi , 

X } = P] = — = — L090 x 10 - ^ in. 

E 

ryr\ 

Xo =-(3, +V2p 2 = — {F\ - lj2F 2 ) = 5.454 xl0~ 3 in. 

E v ' 

9Q F 

X 3 =p 5 = — = 0.909 xlO -3 in. 

E 

X 4 = -(p! - V2P 2 +P 3 ) = “ ( F 1 “ 2 + F 1 ) = 4545 X 10 ' 3 in - 


and solve 


( 3 . 11 . 7 ) 


( 3 . 11 . 8 ) 


( 3 . 11 . 9 ) 
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Load case 2: Solution for thermal loads . — Only the right side of the CC, {SR} , has to be modified for thermal 
analysis: 


{M = -[c]{P} 0 


W 


0 


' 0 " 
0 

at AT 
0 
0 
0 


{5R} = -[C]{P}° 


f 0 1 


0 

1 


0 0 
-V2 1 


0 0 
-■Jl 1 


0 

at AT 
0 
0 
0 


{5R} = 


0 1 

-at A Tj 


0 

-1.20x10' 


For thermal loads only, the IFM equations can be written as 


'l 

l/y/2 

0 

0 

0 

0 " 


h 


0 

0 

1/V2 

1 

0 

0 

0 


f 2 


0 

0 

0 

0 

1/V2 

1 

0 

1 

h 


0 

0 

0 

-1 

-1/V2 

0 

0 

Fa 


0 

0 

0 

0 

0 

0 

1 


f 5 


0 

1 

- 2 V 2 

1 

-2V2 

1 

0 


f 6 


-6.0 xlO 3 


(3.11.10) 


(3.11.11) 


The {8R} in equation (3 . 1 1 . 1 1 ) is normalized with respect to 20/E (see eq. (3.11.6)). Forces due to AT= 1 00 °F in mem- 
ber 3 are obtained by solving the IFM equation as 
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F\ 


—545.45 

Fl 


771.39 

Fs 


-545.45 

Fa 


771.39 

f 5 


-545.45 

F(i 


0 


Calculation of displacements . — Displacements at the nodes are calculated from the DDR as follows: 


=Pl =Pf +P5 

w = ( 5), = - l090x10 ' 3 

Pi=o 

X x = -1.090 xio -3 in. 

Likewise, other displacements can be calculated: 


*1 


Pi 


-1.091 



-Pi + V 2 P 2 


5.454 

x 3 

► - 

Ps 


-1.091 

*4, 


_ (Pl +P 3 ) 


-5.454 


(3.11.13) 


(3.11.14) 


Illustrative Example 12: A Ring Problem 

A uniform circular ring of radius R, moment of inertia /, and Young’s modulus E is subjected to self-equilibrat- 
ing forces as shown in figure 3.12(a). Determine the bending moment, axial force, and vertical displacement at A. 

Step 0 — Solution strategy: For the ring, the polar coordinates (r = R, 0) shown in figure 3.12(a) are used. Only 
half of the ring needs to be considered because of symmetry. The free-body diagram of the half ring is shown in fig- 
ure 3.12(b). The internal forces acting on the ring depicted in figure 3.12(b) at location B are the normal force N B , 
shear force V B , and moment M B , and at location A they are the normal force N A , shear force V A = PI 2, and moment 
M a . The value for V A = PI 2 is obtained from equilibrium considerations, (see fig. 3.13(c)). There are five unknown 
forces (N a , M a , N b , V b , M b ) but only three EE (LN = 0, XV = 0, 1M = 0) can be written, so the problem is two- 
degrees indeterminate. 

Step 1 — Formulate the equilibrium equations: The three EE are as follows: 
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(c) Transverse EE at A. 


Figure 3.12.- A ring problem. 


X"=° 



^ M = 0 at A 


N b -N a =0 or N a = N b 
P P 

Vd 1 — = 0 or V B — 0 

^22 

PR 

M a -—-M b +N b (2R) = 0 


(3.12.1) 


The first two EE (LN = 0 and IV = 0) are trivial in nature. The moment EE in matrix notation can be written as 
follows: 


[-1 -2R 




(3.12.2) 


Two CC are required for the determination of forces. Rotation cp is considered to be the dual, or displacement, 
variable of the EE. 

Step 2— Derive the deformation displacement relations: The three DDR for the problem are 


NASA/TP— 2004-207430 


89 



Pi =-<p 


P 2 = -2 /?<p 

P 3 =<P (3.12.3) 

The deformations (Pj, p o , p 3 ) correspond to the forces ( M A , N g , M g ), respectively. 

Step 3 — Generate the compatibility conditions: Two CC are obtained by eliminating the rotation tp from the 
three DDR: 


p 2 -27?p 1= 0 


p 2 +2/?P 3 = 0 


In matrix notation, the CC can be written as 


-2 R 
0 


1 0 
1 2 R 




The null property of the |5J and [C] matrices ([S][C] r = [0]) can be verified as 


[-1 -2 R 1] 


-2 R 0 

1 1 
0 2 R 


= [o o] 


Step 4 — Formulate the force deformation relations: The FDR for the ring can be determined as 


(3.12.4) 


(3.12.5) 


(3.12.6) 


Pi 


K 

-f 

El J 


3M 4 


P 2 



0 



p 3 


K 

= ~f 

El J 


dMo 


(3.12.7) 
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The bending moment for the problem can be defined in ranges AC and CB as follows: 


DR 

7?t(Q) = M a + N a R(1-cosQ) sin9 intherange AC, 0<9< 


7%(Q) = M A + N A R(\ - cos9) — 


2 

PR 


dm 

dM A 

dm 

dN A 

dm 


= 1 


= /?(1 - cos 0) 


dM 


in the range CB, — < 0 < 7t 


(for both ranges. 0 < 0 < n) 
(for both ranges. 0 < 0 < k) 

(for both ranges, 0 < 0 < ti) 


(3.12.8) 


B 


jR_ 

El 


k /2 )t 

J \m a +N a R(\-cosQ)-— sin0Kl)rf0+ \\m a +N a R(\ -cos 0)-— 1(1)^0 

0 ic/2 


(3.12.9a) 


or 




Likewise, |3, is obtained as 


(3.12.9b) 


P2 = 1j{ kMa + T RNa “T (3+k) 


P 3 =0 


(3.12.10) 


Step 5— Express the compatibility conditions in terms of forces: Since p 3 = 0, the two CC in deformations 
become uncoupled equations as 


Pj = 0 and P 2 = 0 


The CC in forces can be written as 


(3.12.1 1) 


K KR 


n 


3 kR 
2 


M a 

N a 

M b 


PR f n 
— 1 + — 
2 l 2 


PR 


(3 + x) 


(3.12.12) 
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Step 6— Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve 
for forces: 


-1 

-2R 

1 


m a 


-1 

K 

nR 

0 


Na 

H 

H 3 

( 1 + f) 

K 

3nR 

2 

0 


Mb 


( 3 ri 


Solving the IFM equation yields the forces: 



PR 

M a 


4 

N A 

, — , 

P 



2n 

Mr 




lit 4yJ 


(3.12.13) 


(3.12.14) 


Other forces can be back-calculated as N B = N A = PI 2k and V B = 0. 

Step 7 — Back-calculate the displacement, if required, from the deformation displacement relations: We 
bypassed formulating a shear equilibrium equation at load application point A for the calculation of forces. This EE 
is reinstated for the displacement calculation because the dual variable of this EE represents the displacement along 
load P. The transverse EE at A is 


2V a = P or 



(3.12.15) 


The DDR for this EE can be written as 


p v =2v A (3.12.16) 

where v A is the dual variable of the EE representing the transverse displacement at the load application point. The 
deformation corresponding to the shear force V A can be calculated as 


Pv = 


dU 

dV A 


R_ 

El 


2n 

J* 

0 


9% 


dQ 


(3.12.17) 


Because of symmetry, 
or 





El J 


sAe 

dV A 


Tmt _ idm 

dv^~ dP 


(3.12.18a) 


(3.12.18b) 
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Illustrative Example 13: Three-Span Beam Under a Distributed Load 

A three-span continuous beam made of a single material and under a uniformly distributed load is depicted in 
figure 3.13(a). Each span has a length l, a uniform El, and a load intensity q per unit length. Analyze the beam for 
moment and displacement. 

Step 0— Solution strategy: The coordinate system (x, y) with its origin at A is shown in figure 3.13(a). Four 
reactions (R A , R B , R c , R D ) are considered as the force unknowns of the problem, or n = 4. Two independent EE (one 
transverse EE, YV = 0, and one rotational EE, YM - 0) can be written, or m = 2. The beam is two-degrees indetermi- 
nate (r = n-m = 2). We can, however, use the symmetry condition (R A = R D and R B = R c ) to reduce it to a two- 
variable problem. This two-variable problem has one independent EE and one CC. 

Calculation of shear force: Shear forces at the left and right of support B [v ' B , V B ) and, likewise, at support 
C (y', V< S) are defined. These force shear forces can be calculated from the diagrams shown in figure 3.13(b). From 
figure 3.13(b) we can write 

Ra + = or = q( — Ra 

+ R B + Vfs = V- or = q 1 ' ~R a - R b 

Likewise, Vq and v£ can be calculated as 
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Vc — 2ql—R A -R b 


V r c =2ql-R A -2R B 


The free-body diagram for the beam with the shear forces is depicted in figure 3.13(c). 

Step 1 — Formulate the equilibrium equations: The transverse EE along the y-direction yields 

2R a +2R b = 3q( (3.13.1) 

The reader can verify that the moment equilibrium is automatically satisfied. 

In matrix notation, the EE can be written as 



(3.13.2) 


The single EE is expressed in terms of two unknowns ( R A , R B ). Thus, one CC is required to solve the problem. 

Step 2 — Derive the deformation displacement relations: The DDR ({ P } = [B^IX] ) for the problem has the 
following form: 


Pi=-X 

P 2 =-X (3.13.3) 

Here, Pj and P 0 are the deformations associated with reactions R A and R B , respectively. The dual variable of the EE 
is the displacement X. 

Step 3 — Generate the compatibility condition: The single CC is obtained by eliminating the displacement X 
from the two DDR as 


or in matrix notation. 


Pi -p 2 - 0 


The correctness of the CC is verified from its null property ([B][C] r = [0]) as 


(3.13.4a) 


(3.13.4b) 



= [ 0 ] 


Step 4 — Formulate the force deformation relations: For the calculation of the FDR, the free-body diagram 
shown in figure 3.13(c) is considered. The bending moments required to calculate the deformation are obtained for 
the three spans separately. For simplicity, local coordinates are used for the FDR calculation because these are inde- 
pendent of the coordinates. 

Span AB : The free-body diagram of the span with moment, shear force, and external loads is shown in figure 
13(c). The force equilibrium for the span can be verified by summation of the forces and moments. One can deter- 
mine the bending moment (Mf for the first span AB by considering a local coordinate Xp with the origin at A and the 
axis measuring from A to B as shown in figure 3.13(d). 


NASA/TP— 2004-207430 


94 



y 


q per unit length 



(a) Three-span beam. 



(c) Free-body diagram for the beam. 



ii 

■VCQ 


(d) Free-body diagram for the first span, AB. 



|| n 

(e) Free-body diagram for the second span, BC. 
Figure 3.13.- Three-span beam under distributed load. 
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CD-. 







(3.13.5a) 


- R a xj- ■ 


Moment at B: 


Shear force at B: 


V^=q(-R A 


dVit d7X f 

dR A f dRv 


Span BC : For the central span BC, the local coordinate system (x ( ) is selected with its origin at B as shown in 
figure 3.13(e). 


Moment at C: Because of symmetry, Mq = M g = R A I 


Shear force at C: 


v£=2 qt-R A -R l 


The moment function for span BC becomes 


•7K C - M B -(ql-R A - Rg ) x c y~ 


Wc -R A (l + x c )+R B x c - ^ (f 2 +x~ +2fx c j 


cM c „ , d?K c 

— = {l + x c ) and — = x c 

dR A dR B 


Because of symmetry, span CD (which is identical to span AB) need not be considered separately. 
The deformations (3, and P 2 can be written as 


Pl= i7 


El J J dR A 
Span AB 


d 

' zr a 

Span BC 


l \ r r 

[3 1 — — •] 2 I IK f — dx + 

El 1 J } dR D J 


\AB , o BC 


™ c dR^ dX \ = 2 ^ +P2 


The contribution to the deformation P ( - from span CD is accounted for by doubling the value for span AB 
because the first and last spans are symmetrical. Upon integration, the contribution from the first span, or the 

deformation component pj^ B , is obtained as 


O AB ... 1 fw 3 qt 4 

Pl £7 3 8 
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Likewise, the contribution from the central span, or the deformation component p^ c , is obtained as 


Pf = — 

1 \EI 


TR a ^_ + 5R a 


,3 15 qi A ^ 


Total deformation pj is obtained by adding the contributions from the three spans: 


(3.13.9) 


Pi 


El 


3 Ra 


g3 1 5 R B e 17 qt 
6 8 


Likewise, deformation p, can be calculated as 


(3.13.10) 



r 5R A P 


Rb l 3 

3 


17 q r 

24 


(3.13.11) 


In the calculation of deformation p 2 , there is no contribution p, B = 0 from the first (or third) spans because 
dTUJdRg = 0. In other words, the reaction R B at support B is not present explicitly in the moment function TJtj. 
In matrix notation, the deformation force relations can be written as 


p, 


[3 "I 
6 

Ra 

<li A 

17' 

8 

P2 

El 

5 1 

_6 3. 

r b . 

El 

17 

.24. 


(3.13.12) 


Step 5 — Express the compatibility conditions in terms of reactions: In terms of reactions, the CC (pj - P, = 0) 
can be written as 


'13 

6 


r 

2 


I R* 


R 


1 [ 17 qi 


B) 


(3.13.13) 


Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the 1FM equations, and solve 
for forces: 


-1 

-1 

Ra 


' 3 ' 

13 
_ 6 

1 

2. 

r b. 

• = qC 

2 
17 
12 . 


Solution of the IFM equation yields the reactions: 


Ra 1 


l R B) 


\ = ql 


Other forces and moments can be back-calculated from the reactions as 


(3.13.14) 


(3.13.15) 
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(3.13.16) 


Rc 


-Rb - 


11 ql 
10 


Rd 



The moment in the first span is given by 




2q£xf 

5 


2 

f_ 

2 


(3.13.17) 


The moment function ^ can be used for the third span; however, the coordinate (xy) has to be measured with its 
origin at D and with its abscissa from D to C. The moment in the central span is 

7K r = - + ^ Xc - -- q X }. (here.*:.. is measured from B to C, see fig. 3. 13(e)) (3.13.18) 

c 10 2 2 c V 

Step 7 — Calculate the displacement, if required, by integrating the moment curvature relations: 

For span AB, 


d 2 wj Wf 
dx 2 = ~eF 


1 

£7 


( 

2 

~qtx f -q 

V 



2 J 


Upon integration, 


El 


si. 

15' 


d 

24 


+ c l x f +c 2 


Displacement uy = 0 at xy = 0 and xy= l , 


or 


and 


= 0 

and 


1 ql 3 

c \ ~ 

El 40 



1 

q£x 3 

4 

qx f 

ql 5 x j 

~~El 

15 

24 

40 


(3.13.19) 


(3.13.20) 


(3.13.21a) 


(3.13.21b) 


The displacement function uyfor span AB can be used for the third span CD by measuring Ay from D as mentioned 
earlier. 


For the span BC, 


Upon integration, 


d 2 w c 

dx 2 


= i { -qt 1 | Wc 

El El 10 2 2 

V / 


El 


q£x 3 


4 ■) 
qx c 


20 


12 24 


+ c\x c +c 2 


(3.13.22) 


(3.13.23) 
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The displacement w = 0 at x c = 0 and x c -£ yields the integration constants as 


c 2=0 and 


<?l 3 

120 El 


The displacement function for the central span has the following form: 


El 


i 2 xl_ + qtxl 


20 


12 


4 .3 \ 

qx c | q£ x c 

24 120 


(3.13.24) 


(3.13.25) 


Illustrative Example 14: Portal Frame 

A steel portal frame has the geometrical dimensions and load shown in figure 3.14(a). Analyze the problem for 
moments, displacements, and rotations under load at C. The Young’s modulus of elasticity E is 30 000 ksi, and the 
moments of inertia I are given in figure 3.14(a). 

Step 0 — Solution strategy: The coordinate system (x, y) with its origin at A is shown in figure 3.14(a). For 
analysis, the portal frame is divided into four elements and five nodes (A, B, C, D, E). It has a total of eight moment 
unknowns (Afj, M 2 . . ., M g ), or n = 8. For simplicity, the five displacement unknowns (m = 5) considered are 

( 1 ) Displacement X s , representing the sway of the portal frame (see fig. 3 . 1 4(a)) 

(2) Three rotations (0 fi , Q c , 0 D ) at three locations (B, C, D). respectively 

(3) Transverse displacement X c at load application point C 

The problem is three-degrees indeterminate (r = n-m- 3). For this analysis, only bending deformations are consid- 
ered; axial deformations are neglected. 

Step 1 — Formulate the equilibrium equations: The free-body diagram shown in figure 3.14(b) is sufficient to 
write the EE for the problem. The five EE along the five displacement degrees of freedom can be written as follows: 

(1) Along the sway X s direction — The displacement along the y-direction at the top of the frame or along the 
beam BCD is referred to as the sway displacement^. The force equilibrium along the sway direction for the beam 
BCD yields the following EE: 


M i —M2 M% — Mj 

T\ 

(2) Along the rotation 0 B direction at B: M 2 - M 3 = 0 

(3) Along the rotation 0 C direction at C: M 4 -M 5 = 0 

(4) Along the rotation 0 D direction at D: M 6 - Mj = 0 

(5) Along the transverse displacement X ( . direction at C: 


(3.14.1) 

(3.14.2) 

(3.14.3) 

(3.14.4) 


M 3 — M 4 M 6 


-M* 


- + P 


(3.14.5) 


The EE are simplified by eliminating the length parameters = t A and i 2 = 2^ 3 ). The process makes the EE 
dimensionless, and the load term (P 5 ) becomes an equivalent moment of 2592.0 kip-in. The five EE can be written in 
matrix notation as 
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■ 1 
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0 
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1 

1 " 

Mf 

-\ 

h 

h 


M z 

0 
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-1 

0 

0 

0 

0 

0 

m 3 






0 

0 

0 

Ma 

0 

0 

0 

1 

-1 










M S 

0 

0 

0 

0 

0 

1 

-1 

0 

m 6 



1 

1 

2 

2 

0 

0 

m 7 

0 

0 





m 8 _ 



~2 

~~2 

^2 

^2 




0 

0 

• 0 • 

0 

-18 


(3.14.6) 


Since five EE are expressed in terms of eight unknown forces, three CC are required for their determination. 
Step 2 — Derive the deformation displacement relations: The eight DDR for the problem are as follows: 


Pl=- 


P2 - — 


P5 - ~6c + 
p6 =e D - 


2Xc 
2 

2X C 


p 3 =-0fl — - 

3 h 

X c 

p 4 =0 C+ _C 

h 


P 7 

x, 

Ps = 

l x 


(3.14.7) 


In the DDR, the eight deformations (Pj, (3 2 , . . . , P 8 ) correspond to the eight moments (M v M 2 , . . . , M g ), 
respectively. The dual displacement variables of the five EE are X s , 0 g , 0 C , 0 D , andX" c . 

Step 3 — Generate the compatibility conditions: The five displacements are eliminated from the eight DDR to 
obtain three CC: 


1 1 

-1 0 

1 0 


1 1/3 1/3 0 

0 2/3 2/3 1 

0 0 0 0 


0 0 
1 0 
0 1 


Pi 

P2 

1 P3 
P 4 

P 5 
J P 6 
P7 
Pg 


O' 

■ 0 - 

0 

J 


The null property {[E\[C\ T = 0) of the compatibility and equilibrium matrices can be verified as 


(3.14.8) 


NASA/TP— 2004-207430 


101 



— — — 0 0 0 0 — — — r 1 -! 1 

i o 0 

01 -1 0 0000 100 

„ „ „ 1/3 2/3 0 

0 0 0 1 -1 0 0 0 

1/3 2/3 0 

0000 01 -1 0 010 


00 0 0 — -— 
li l x 

-1 0 0 0 0 0 


0 0 0 1 -1 0 0 0 


0 0 


_ 2 _ _ 2 _ 

e 2 e 2 




0 1 0 
0 0 1 


[0 0 0 
0 0 0 
= 000 
0 0 0 
0 0 0 


Step 4— Formulate the force deformation relation: The FDR for a uniform beam of span t, subjected to end mo- 
ments A/, and Mj, and the corresponding deformations P ( and can be written (see eq. (2.48c)) as 

fp/i ( t Y2 iir^i 


U-Lf 


2 M, 


The FDR can be specialized for the four beam elements of the problem as 


Pi = — (2.4M 1+ 1.2M 2 ) 
6 E 


p 5 = — (0.4M 5 +0.2M 6 ) 

oh, 


p 2 =-2-(1.2M 1+ 2.4M 2 ) 
6 E 


P 6 =2-(0.2M 5+ 0.4M 6 ) 
6 E 


P 3 = — (0.8 Mj +0.4M 4 ) 
6E 


p 7 =— (1.6M 7 +0.8Mg) 
6 E 


p 4 =— (0.4M 3 +0.8M 4 ) p 8 =— (0.8Af 7 +1.6M 8 ) (3 

6E 6 E 

Step 5— Express the compatibility conditions in terms of forces: The CC can be expressed in moments by 
eliminating deformations between the CC and the FDR: 


3.600 

3.600 

0.933 

0.667 

0.133 

0.067 

0.000 

0.000 

-2.400 

-1.200 

0.267 

0.533 

0.467 

0.533 

1.600 

0.800 

2.400 

1.200 

0.000 

0.000 

0.000 

0.000 

0.800 

1.600 


Step 6 — Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and 
solve for forces: 


NASA/TP— 2004-207430 


102 



1 

1 





1 

1 ' 





0 

0 

0 

0 


— 









h 

h 



0 

1 

-1 

0 

0 

0 

0 

0 
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-1 
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0 

0 


M 2 
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m 3 


0 

0 

0 

0 

0 

0 

1 

-1 

0 


m 4 
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0 

0 

1 

1 

2 

2 

0 

0 

< 

m 5 

► = < 

-18 



1 2 

f-2 

^2 





m 6 


0 

3.600 

3.600 

0.933 

0.667 

0.133 

0.067 

0.000 

0.000 


M-j 


0 










Mg 


0 

-2.400 

-1.200 

0.267 

0.533 

0.467 

0.533 

1.600 

0.800 



2.400 

1.200 

0.000 

0.000 

0.000 

0.000 

0.800 

1.600 




Solving the IFM equation yields the moments as 


Mi ' 


136.1 

m 2 


-213.3 

m 3 


-213.3 

m 4 


617.9 

m 5 


617.9 

m 6 


-262.4 

m 7 


-262.4 

M s 


87.0 


Reactions can be back-calculated as follows. 
Horizontal reactions (see fig. 3.14(b)): 


At support A. the reaction H 4 is along the positive x- direction: 

M 2 — Mi 


Hj = - 


= -1.12 kip 


At support B, the reaction H B is along the negative x- direction: 

M-j — Mg 


Hn =- 


= -1.12 kip 
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(3.14.13) 


(3.14.14) 


(3.14.15a) 


(3.14.15b) 



Vertical reactions: 


At support A, the vertical reaction V A can be calculated as 


V A 



5.77 kip 


(3.14.16a) 


At support E. the vertical reaction V E becomes 


V E 


M5 - Mg 

h 


= 12.23 kip 


(3.14.16b) 


Step 7 — Back-calculate the displacement, if required, from the deformation displacement relations: The 
displacement under load P is 

~ ~ = (Pl +P 2 + P3 ) 

*■2 

X c = — (3.6Mj +3.6M 2 +O.8M3 +0.4M 4 ) = -161.1 xlO -3 in. (3.14.17) 

6 E 

Displacement^ is along the negative y-axis, which is the direction of the external load. The sway displacement is 

X = l$ x = — (2.4Mj + 1.2M 2 ) = 122.5 xlO -3 in. (3.14.18) 

6 E 

The portal sways along the positive x-axis. 


Illustrative Example 15: Navier’s Table Problem 

Structural indeterminacy was recognized by Navier (1785-1 836) when he attempted to determine the four reac- 
tions (ftj, R 2 , Ry Rf) along the four legs of a table that was subjected to a concentrated load P with eccentricities e x 
and e as shown in figure 3.15(a). He assumed that the symmetrical table was made of wood and that it was resting 
on a level floor made of a rigid material, such as stone. The distances between the legs along the x- andy-directions 
were 2a and 2b, respectively. Solve the problem for the reactions and displacements. 

Step 0 — Solution strategy: For the table problem, the origin A of the coordinate system (x ,y, z) is selected at the 
center of the table top. The reactions (R v R 2 , Ry R 4 ) are the four force unknowns, or n = 4. The problem has three 
equilibrium equations and three displacement variables (w, 0 X , Q y ) — shown in figures 3.15(b), (c), and (d) — which 
represent the dual variables of the three EE, or m = 3. The problem is one-degree indeterminate {r-n- m =1). For 
the table problem, the reaction along the positive y-direction produces compression in the legs of the table as shown 
in figure 3.15(a). 

Step 1 — Formulate the equilibrium equations: The three EE (the sum of the reactions, or EF = 0, and the sum of 
the moments along the x- andy-directions, or I X M= 0 and 2 y M= 0) can be written for the table problem as follows: 
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or 


(3.15.1) 


X z V = 0 -1 -1 
X x M = 0 a a 


Ij,M = 0 


-b 


b 


-1 

■*f 


-P 

< 

R 2 

► = < 

~e„P 


R’i 


y 




c P 


*4, 


l x J 


[B]{R} = {P] 

where [S] is the (3 x 4) equilibrium matrix, and {/?} and {P} represent the four reactions and the three load compo- 
nents. The three EE are expressed in terms of four unknown reactions, so the table problem cannot be determined 
from the EE alone. Navier was the first to recognize the indeterminate nature of this problem, which is one-degree 
indeterminate. One CC is required for its solution. 

Step 2— Derive the deformation displacement relations: The DDR ({ p} = [B] r {X)) for the table problem can be 
written as 


7V 


-1 

a 

-b 

P 2 


-1 

a 

b 

£3 


-1 

~ a 

b 

Pa. 


-1 

-a 

-b 


1 . . 


w 

- 

0* 




(3.15.2) 


where the four deformations ((3j, P 2 , P 3 , p 4 ) along the four legs of the table correspond to the four reactions (R V R 2 , 
R 3 , /? 4 ), respectively. The three displacements {X}, which are the dual variables of the EE, represent one translation, 
w, along the r-direction and two rotations, Q x and 0 , about the x- and y-axes, respectively. 

Step 3— Generate the compatibility condition:' The DDR contains four equations in terms of three displace- 
ments. In other words, the four deformations are not independent. One relation between the four deformations can 
be obtained by eliminating the three displacements from the four DDR: 


or 


Pi -P2 +P 3 — P4 - 0 


[1 -1 1 -1] 


P 2 

P 3 

IP4J 


= 0 


(3.15.3) 


Equation (3.15.3) represents the deformation compatibility condition ([C]{p} = {0}) for the table problem. Like the 
EE, the CC is independent of the material of the structure. 

The null property of the equilibrium and compatibility matrices ([B][C] r = [0]) can be verified as 


-1 -1 -1 -1 
a a —a -a 
—bb b -b 


1 

' 1 ■ 


o' 


-1 



0 


1 




-1 


0 





(3.15.4) 
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p 






Figure 3.15.- IFM introduced through Navier's table problem. 
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Step 4 — Formulate the force deformation relations: The FDR for the problem can be written as 


Pi 


£R\ 

~AE 


Pi 


£R 2 

AE 


P 3 


m 3 

AE 


P4 


(R 4 

AE 


or 


£ Pi P2 P 3 P4 

= EL = E± = ZJ-~r± = constant 

AE R ] R 2 R 3 R 4 


(3.15.5) 


where l is the table height, A is the cross-sectional area of the legs, and E is the modulus of elasticity of the material 
of the table legs. 

Step 5 — Express the compatibility conditions in terms of forces: The CC in forces is obtained by eliminating 
deformations p in favor of reactions R as 


AE 


(/?! — R 2 + Rn, — R 4 ) — 0 


or 



-1 1 



(3.15.6) 


Step 6— Couple the equilibrium equations and compatibility conditions to obtain the IFM equations, and solve 
for forces: 


'-1 

-1 

-1 

-1 




-p ' 

a 

a 

—a 

—a 


R 2 


-e y P 





< 


> = * 


-b 

b 

b 

-b 


R 3 


e x P 

1 

-1 

1 

-1 


* 4 . 


0 


(3.15.7) 


In the homogeneous CC given by the fourth equation in (3.15.7), {l /AE = 1) is set to unity. Solution of equation 

(3.15.7) yields the forces: 


where r. = e,Jb and r v = e v la. 

a a y y 


Ri 


>■ 

1 

l 


r 

r 2 

P 

l 

+ 

P 

1 

> — < 

\ + r x + r y 

> ~ 


*3 

4 

4 

1 

*4. 


& 

+ 

1 


1 


r =r= 0 


(3.15.8) 
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When the load is placed at the center of the table, r x = r = 0.0 and each leg carries one-quarter of the load. 


*1 


- R 2 - R$ - R 4 


P_ 

4 


(3.15.9) 


Step 7 — Back-calculate the displacement, if required, from the deformation displacement relations: Once inter- 
nal forces are known, displacements can be calculated by back-substitution in the DDR given by equation (3.15.2) 
and the FDR. 


Pi + P2 + P3 +P 4 


4 AE 


(R l +R 2 +R 2 +R 4 )~- 


Pl 
4 AE 


Pi + P2 ~P 3 ~ P4 

4 a 


4 AEa 


{R1+R2-R3-R4)-- 


Pte 


y 


4 AEa 1 


e,=- 


Pi -P2 ~P 3 +P 4 

Ah 


4 AEb 


(/?l ~R 2 ~R 3 +R 4 )- 


PZe x 
4 AEb 7 


(3.15.10) 


When the load is placed at the center of the table (r x = r y = 0), the average displacements become 


f \ 


' e ' 

w 

p 

AE 

Q x 

~4 

0 



0 

v y) 

O 

!! 

II 



The three displacements — uniform displacement (w) along the transverse direction and two rotations (0 r 8 ) are 
depicted in figures 3.15(b), (c), and (d). 

The transverse displacement (w) is along the negative y-direction, which also represents the load direction. Like- 
wise, the rotations (0^, 0 ) are along the negative (clockwise) and positive (counterclockwise) directions, respectively 
(see figs. 3.15(a), (b), and (c)). 

In the solution of the table problem, the tabletop is assumed to be rigid, which can be readily visualized by an 
examination of the displacements given by equation (3.15.10) and figures 3.15(b), (c), and (d). The assumption limits 
the scope of the analysis to a certain extent. 
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Chapter 4 

Integrated Force Method and Dual 
Integrated Force Method for 
Finite Element Analysis 

Internal forces and displacements are the primary unknowns of discrete structures, including frameworks (that is, 
trusses and frames) and continuous structures (such as plates, shells, and solids) idealized by finite elements. For the 
purpose of analysis, such structural models can be designated by two attributes, n and m, such as “structure (n, tri)P 
The number offeree unknowns or force degrees of freedom Cfof) is n. Likewise, the number of displacement 
degrees of freedom ( m dof) is m. Once the n forces {F} are determined, then the m displacements {X} can be back- 
calculated from known forces, and vice versa. The method that treats all n fof as the principal unknowns is known as 
the force method. Likewise, the method that treats all m dof as the principal unknowns is called the displacement 
method. The force method and displacement method are the two fundamental structural mechanics formulations. 

The displacement method, also known as the stiffness method (refs. 33 to 35), has been well researched and 
developed during the past few decades. This method currently dominates the analysis scenario. Its governing equa- 
tion is (|X]{X} = {P}), where \K\ is the stiffness matrix, (X) is the displacement vector, and {P\ is the load vector. 
The method is indirect because forces and stresses are back-calculated from the displacements, and the derived quan- 
tities can be susceptible to inaccuracies. Yet, the method is popular because of its simplicity, versatility, and com- 
puter amenability. The stiffness method parallels the Navier’s displacement formulation in elasticity (see table I in 
the Preface). 

The force method is the direct force determination formulation for structures. One would have anticipated 
that all n internal forces {F} would have been the primary unknowns of the force method with basic equations 
([5] { F} = { P * }), where [S] is the governing matrix, and that the force method in structures would have paralleled the 
Beltrami-Michell Formulation in elasticity (see table I). However, the classical force method, with redundants as the 
unknowns, satisfies neither attribute. Although the redundant method is elegant for manual computation of small 
problems for static loads, the classical formulation is cumbersome for computer automation; and for all practical 
purposes, it has disappeared from current use. The redundant force method, which has considerable historical signifi- 
cance, is described in references 30, 36, and 37. 

A force method that could be considered equivalent to the Beltrami-Michell stress formulation in elasticity was 
not available because the understanding and the development of the compatibility conditions (CC) were immature. 
We have researched and have come to understand the CC of structural mechanics. We express the CC for discrete 
analysis first in terms of deformations { P } as [Cl { p } = 0 (here [C\ is the compatibility matrix), then in terms of all 
( n fof) force variables as [C][G]|F} = {8R} (here [G] is the flexibility matrix and {8/?} is the effective initial defor- 
mation vector). Thus, the classical, redundant-based, ad hoc compatibility has been replaced by the correct 
CC representing the deformation balance conditions. The CC are coupled next with the equilibrium equations 
[B]{F] = {P) (here [/?] is the equilibrium matrix, {F} is the force vector, and \P\ is the load vector), to obtain the 
direct force determination method, which in the literature is referred to as the Integrated Force Method (IFM). The 
IFM with forces {F} as the unknowns parallels the Beltrami-Michell Formulation in elasticity, with the governing 
equation [S]{F] = {P}*. The IFM is as versatile as the stiffness method, and it is amenable to computer automation. 
It produces accurate stress and displacement solutions, even for modest finite element models. 

A dual formulation for the primal IFM, referred to as the Dual Integrated Force Method (IFMD), has been for- 
mulated. The dual method uses the elemental matrices of IFM, but displacements are its primal unknowns with the 
governing matrix ([D]{XJ = {F}). 
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Forces {F} in the dual method are backcalculated as {F} = \G\~' L \R] T {X} . The primal IFM and dual IFMD are 
equivalent formulations and produce identical results for forces and displacements. In the earlier chapters, the IFM 
was developed for elementary, indeterminate structural mechanics problems. In this chapter, the primal and the dual 
methods are introduced for finite element analysis. 


Equations of the Integrated Force Method 

The IFM equations for a continuum discretized by finite elements with n and m force and displacement degrees 
of freedom, respectively, are obtained by coupling the m EE and the ( r — n - m) CC. The m equilibrium equations 
([5]{FJ = {P}) and the r compatibility conditions ([C][G]{F} = {5F}), when combined to obtain the IFM governing 
equations for static analysis, can be written as 


r t b ] ■ 

m 


{F}= — 

[C][G]_ 

{»} _ 


or [5]{Fj = {P *} 


(4.1) 


From forces {F}, displacements {X} are back-calculated with the following formula: 

{X} = [y]{[G]{F} + {p 0 }} (4.2) 

where [2] = m rows of [[S] -i j . 

The definitions of the matrices and vectors in equations (4.1) and (4.2) follow: 

[S] is the ( m x rt) rectangular equilibrium matrix with more columns than rows. It is a very sparse, unsym- 
metrical matrix with full row rank. 

[G] is the ( n x n) symmetrical flexibility matrix. It is a block-diagonal matrix, where each block represents a 
flexibility matrix for an element. 

[C] is the (r x n) compatibility matrix. The generation and properties of the compatibility matrix are explained 
later in this chapter. {8F} = -[C]{(3}° is the effective initial deformation vector. Here, {|3} 0 is the initial deformation 
vector of dimension n. 

[5] is the (n x n) IFM governing unsymmetrical matrix defined in equation (4.1). 

[/] is the ( m x n) deformation coefficient matrix back-calculated from the [5] matrix. 

The IFM has two key equations: equation (4.1) to calculate forces, and equation (4.2) to calculate displacements. 

The familiar process of differentiation used in the popular stiffness method to generate stresses from displacements 
is avoided in the IFM. 

In terms of fundamental operators, an analogy can be made between the IFM and the theory of elasticity 
(ref. 38). The three fundamental operators of elasticity are (1) the equilibrium operator of Cauchy, which relates 
stresses to external loads; (2) the compatibility operator of St. Venant, which controls components of strain; and 
(3) the material constitutive matrix of Hooke, which relate strains to stresses. Likewise, the IFM has three matrices 
that are equivalent to the operators of elasticity theory. These operators, which become matrices in the context of 
finite element analysis, are (1) the equilibrium matrix [5], which links internal forces to external loads; (2) the com- 
patibility matrix [C], which governs the deformations ([CJ{(3} = {0}); and (3) the flexibility matrix [G], which relates 
deformations and forces. Both the equilibrium and the compatibility operators of elasticity and the corresponding 
matrices of the IFM are unsymmetrical, whereas the material constitutive matrix and the flexibility matrix are sym- 
metrical. 

Governing operators of other formations (e.g., Navier’s displacement formulation, Airy’s stress function formu- 
lation, Reissner’s hybrid formulation, or the Hu-Washizu’s mixed formulation) and the matrices of other discrete 
analysis methods (such as the stiffness, redundant force, mixed, and hybrid methods, see table I) are, in principle, 
derivable from the basic unsymmetrical operators of elasticity and the matrices of the IFM (see app. A). Mathemati- 
cally speaking, the derived operators and matrices of other formulations can possess characteristics (i.e., the numeri- 
cal norms, spectral radii, and stability of equation systems) no more superior than the basic unsymmetrical operators 
of elasticity theory or the matrices of the IFM, even when the derived operators and matrices become symmetrical 
(ref. 30). 
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The frequency analysis equation of the IFM without damping is as follows: 


[S]-co 2 


[M][J][G] 

0 


(4.3) 


where [M] is the (m x m) mass matrix, 00 is the circular frequency, and (F} is the force mode shape of the eigenvalue 
problem. 

Forces are the unknowns of the IFM vibration analysis. Displacement modes in IFM, if required, can be back- 
calculated from forces {F} by using equation (4.2). In other words, the IFM provides one set of equations to deter- 
mine forces (i.e., eq. (4.1) for static analysis or eq. (4.3) for vibration analysis) and provides another set for the 
calculation of displacements (eq. (4.2)). 


Equations of the Dual Integrated Force Method 

The Dual Integrated Force Method (IFMD) is obtained by mapping forces into displacements (ref. 5). The basic 
equations of the dual formulation, without initial deformations and damping, are summarized next. 

Static analysis equations of the IFMD are as follows: 

[D]{X} = {P} (4.4) 

Forces can be obtained from displacements by using the following formulas: 

{F} = [G] _1 [F] r {X} (4.5) 

The dynamic analysis equations of the IFMD follow: 

[[D]-co 2 [M]]{X} = 0 (4.6) 

where the (in x m) symmetrical matrix {\D\ = [B][G] _1 [B] r ) is assembled at the element level. 

From the displacement modes, force mode shapes can be back-calculated by using equation (4.5). Like EFM, 
the dual IFMD, which treats displacements as the primary variables, has one equation to calculate displacements (eq. 
(4.4)) for static analysis or eq. (4.6) for dynamic analysis) and has one equation for the determination of forces from 
displacements (eq. (4.5)). IFM and IFMD provide identical solutions for stresses, displacements, and frequencies. 

For design and sensitivity analysis, the primal IFM, however, has some advantage over the dual IFMD (refs. 34 and 35). 


Matrices of the Integrated Force Method 

Three matrices are required for IFM finite element analysis — namely the equilibrium matrix [5], the flexibility 
matrix [G], and the compatibility matrix [C]. The generation of the matrices is illustrated by considering a simple 
rectangular plate-bending element as an example. 

Generation of the equilibrium matrix [B ]. — The EE, written in terms of forces at the grid points of a finite ele- 
ment model, represent the vectorial summation of n internal forces { F} and m external loads {F} . The nodal EE in 
matrix notation give rise to a (m x n)-banded rectangular equilibrium matrix [5], which is independent of the material 
properties and design parameters of the indeterminate structure (n, m). For finite element analysis, this matrix is 
assembled from elemental equilibrium matrices. 

The elemental equilibrium matrices [ B e ] for bar and beam elements can be obtained from the direct force bal- 
ance principle (ref. 37). For continuous structures, such as plates or shells, very few equilibrium matrices have 
been reported in the literature (refs. 36 and 39). Equilibrium matrices for the plate flexure problem are given by 
Przemieniecki (ref. 36) and Robinson (ref. 39). Przemieniecki generates the matrix for a rectangular plate in flexure 
from direct application of the force balance principle at the nodes. Robinson utilizes the concept of virtual work to 
derive the matrix for a rectangular-plate in flexure. The procedures of Przemieniecki and Robinson, which are 
documented in their books (refs. 36 and 39), are not repeated here. 
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Energy-equivalent equilibrium matrices for finite element analysis can be obtained from the IFM variational 
functional (ref. 8). The procedure for generating an elemental equilibrium matrix from the IFM variational functional 
is illustrated next for a rectangular-plate in flexure. The portion of the IFM functional (ref. 8) that yields the equilib- 
rium matrix [B] for a plate flexure has the following explicit form: 

f f d~w r T 

li)is (4 ' 7) 

where \M\ T = (M x , M y , M xy ) are the plate-bending moments, and {e} r = (d 2 w/dx 2 , d 2 w/dy 2 , d 2 w/dxdy ) represent the 
curvatures. The plate domain is D, and the coordinates are x and y. 

By appropriate choice of force and displacement functions, one can discretize the energy scalar U p to obtain the 
elemental equilibrium matrix [B e ]: 

U p ={X} T [B e ]{F} (4.8) 

where the elemental displacement degrees of freedom are symbolized by {X} and the elemental force degrees of 
freedom by {F}. 

The force fields have to satisfy two mandatory requirements: 

(1) The force fields must satisfy the homogeneous EE (here, the plate -bending equations in the element domain). 

(2) The force components F k (k = 1,2 ,...,m) must be independent of one another. This condition ensures the kine- 
matic stability of the element, and the resulting [ B e ] matrix has a full column rank. 

Consider a four-node, rectangular plate bending element of length (a, b) along the x- and y-axes, as shown in 
figures 4.1(a) and (b). For the plate element, 12 equilibrium equations along directions Sj, s 2 , ■ ■ ., s 12 (which corre- 
spond to m dof= 12 for the element, as shown in fig. 4.1(b)) can be written. The internal element forces \ F) (which 
correspond to n fof= 9 for the problem) are different from the 12 nodal directions (5, ,S 2 ,.. ., 5 12 ). 



(a) Nodal displacements. (b) Nodal equilibrium. 


Figure 4.1 .- A rectangular element (PLB4SP) in flexure. 
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For the rectangular plate, the force field is chosen in terms of nine independent forces as 


{F} = {F i ,F 2 ,...,F 9 ) 


as 


M x = F x + F 2 x + F 2 y + F 4 xy 


M y = F 5 +F 6 x + F 7 y + F%xy (4.9) 

M xy = F 9 

Although the variation of the normal moments in the field is linear, the twisting moment is constant. The assumed 
moments satisfy the mandatory requirements. The selection of force fields for advanced finite elements is given in 
references 26 and 27. 

The displacement field that should satisfy the continuity condition (ref. 34) is selected in terms of 12 variables to 
match the three nodal degrees of freedom (the transverse nodal displacement vv ; and two rotations, Q xi and 0 v; , per 
node i for the four nodes). It can be written in terms of Hermite polynomials as 


w(*. y) = 01 (x)H 0 i (y)X J + H 0l {x)H x ! (y)X 2 


+ H x , (x)H 0} (y)X 3 + H 01 (x)// 02 {y)X A 
+ ^02 (y)x 7 + H 02 {x)H 12 (y)X 8 

+ ffl2 ( x )H() 2 (y)X 9 +Hq 2 ( x)H 0l (y)A 10 

+ H 01 (x)H n (y)X n +H n {x)H m {y)X n (4.10a) 

In equation (4.10a), the Hermite polynomials are defined as 


jr 3 -3a 2 x + 2a 3 


H m W - 3 

4 a 


TI , . x 3 — 3 a~x — 2 u 3 

H 02 (x) = 71 

4 a 


Hn(*) = 


x 3 -ax 2 -a 2 x + a i 
4a 2 


H n {x)-- 


x 3 - ax 2 -a 2 x — a 3 
4a 2 


(4.10b) 


where A^, X^,..., X l0 are the 12 degrees of freedom and a and b are the dimensions of the plate element along thex 
and y-directions, respectively (fig. 4. 1 a). The Hermite polynomials for the y-direction can be obtained by changing 
(x and a) to (y and b), respectively, in equation (4.10b). The displacement field equation (eq. (4.10a)) gives rise to a 
linear force distribution for the plate-bending problem. 
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For finite element calculations, the force and displacement fields given by equations (4.9) and (4.10) can be 
written as follows: 


{M} = [Y]{F e } (4.11) 

{w} = [i¥]{l7 e } (4.12) 

Here, {w} = (w), and \M} r = (M x , M „M ) are the displacement and moments at a location within the plate domain; 

{ U } is the vector of element nodal displacements; [N] is the matrix of displacement interpolation functions; { F} e is 
the unknown force vector; and [Y] is the matrix of the moment interpolation functions. The curvature vector { e } , 
which is obtained by differentiation of the displacement field, can be written as 

{s} = [Z]{U e } (4.13) 

where [Z] = [L][A/], and [L] is the matrix of differential operators. 

The expression for the element equilibrium matrix [B e ] is obtained from the strain energy given by equation 
(4.7) by substituting equations (4.11) and (4.13) into equation (4.7). The element equilibrium matrix [ B e ] can be 
symbolized as 


[B e ] = \ D [Z][Y]ds (4.14) 

The equilibrium matrix is obtained by substituting interpolation functions for moments and displacement, then inte- 
grating. For the simple element, the integration can be completed in closed form. Equation (4.15) gives the equilib- 
rium matrix [B e ] thus obtained. The generation of the equilibrium matrix \B e ] illustrated is a general procedure that 
can be applied to any other element type. Henceforth, the matrix obtained from the strain energy function is referred 
to as the consistent equilibrium matrix. 
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(4.15) 



The row and column dimensions of the equilibrium matrix [^correspond to the elemental displacement and force 
degrees of freedom (here m- 12 and n = 9). The equilibrium matrix [B e ] has a full column rank of n = 9. 

Generation of the flexibility matrix [G]. — Generation of the flexibility matrix for the finite element analysis is 
well established in the literature. It is obtained by discretization of the complementary strain energy expression t/ : 

U c =i( {M} T [D]{M}ds (4.16) 

Z J D 

where [D] is the compliance matrix of the material. Substituting the moment interpolation function given by equation 
(4.1 1) into equation (4.16) yields U in terms of the flexibility matrix and force vector: 

U c=\{FeY[G e ]{Fe} (4.17) 

where the element flexibility matrix [G e ] is 

[G e ] = $ D [Y] T [D][Y]ds (4.18) 

The complementary energy for this plate element in terms of moments takes the following form: 

u c = J {m 2 x + M 2 - 2vM x M y + (1 + v)M% } dx dy (4.19) 

where D { = (Eh 3 /]2); E is the Young’s modulus, v is the Poisson’s ratio, and h is the plate thickness. The (9 x 9) 
symmetrical flexibility matrix, which can be integrated in closed form, is obtained as 
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(4.20) 


2(1 + v) 


The symmetrical ( n x n) flexibility matrix [G e ] has a full rank of n = 9. 

Generation of the compatibility matrix [C ]. — The compatibility conditions are controller types of relations. 

For finite element discrete analysis, the balancing of deformations (/(Pj, P 7 , . . . , P n ) = 0) is the central compatibility 
concept. The procedure given earlier to generate the equilibrium and flexibility matrices cannot be applied for the 
generation of the CC. Instead, the CC and the associated coefficient matrix [C] are obtained as an extension of 
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St. Venant’s strain formulation in elasticity. The strain formulation is reviewed first through the plane stress elasticity 
problem and then it is extended to finite element analysis. The strain displacement relations for the elasticity problem 
are 


du 3v 3 m 3v 

E * = a^ Ey_ aJ lxy ~Jty + dx 


(4.21a) 


where £ , £ , and y are the strain components, and u and v are the displacements. The CC in elasticity are obtained 
by eliminating the two displacements from the three strain displacement relations: 


d 2 e x d 2 e y d 2 y xy 

dy 2 dx 2 dxdy 

The two steps that generate the CC are 


(4.21b) 


(1) Establish the strain displacement relations given by equation (4.21a). 

(2) Eliminate displacements from the strain displacement relations to obtain the CC given by equation (4.21b). 

In the mechanics of discrete structures, the deformation displacement relations (DDR) are equivalent to the strain 
displacement relations in elasticity, and the deformations } ji} of discrete analysis are analogous to the strains { £ } 
in elasticity. We obtain the DDR by utilizing the equality relation between internal strain energy (l/2{F} r {p}) and 
external work (1/2 {X} T {P}\ which for a discrete structure (n, m) can be written as 


}W 7 '{P} = !« 7 >} (422a > 

We can rewrite equation (4.22a) by eliminating the load vector {P} in favor of internal forces {P}, by using the EE 
([P]{P} = {P}) to obtain 


\{X} T [B]{F} = ±{F} T m 

or i{F} r {[fi] r {X}-{p}} = 0 (4.22b) 

Since the force { F} is arbitrary and it is not a null vector, we finally obtain the following relation between defor- 
mations and nodal displacements: 

m = [B] T {X} (4-23) 

Equation (4.23) represents the global DDR of a finite element model whose system EE are given as [5]{F} = {P}. 

In equation (4.23), n deformations { P } are expressed in terms of m displacements { X ] ; thus, there are r = n - m 
constraints on deformations, which represent the r CC of the structure (n, m). We can obtain the r CC by eliminating 
the m displacements from the n DDR. In matrix notation, the CC can be written as 

[C]{p} = {0} (4.24a) 

The deformation (P) in equation (4.24a) represents total deformation, consisting of initial deformations (P}° 
and elastic deformations (P} e as { P } = {p) e + (P) 0 . The CC, in terms of the elastic deformations 1PT', are as 
follows: 

[C]{Pr={SF} and {5F} = -[C]{p}° (4.24b) 

Since { P) e = [G]{F}, the compatibility condition in forces can be written as [C][G]{F} = {5 R}. 

The matrix [ C] has the dimensions (r x ri). It is rectangular and banded, with a full row rank of r. The CC are 
kinematic relationships, and these are independent of sizing design parameters (such as the area of the bars, the 
moments of inertia of the beams, and the thickness of the plates), material properties, and external loads. The CC 
depend on the initial deformation in the structure. 
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An important attribute of compatibility conditions .— Because the compatibility condition in total deformations 
is a homogeneous equation ([C]{ |3} = {0}), an alternate matrix [c] obtained by linear combination of the rows of the 
matrix [Q still represents a compatibility matrix that can be used in IFM analysis. The concept is illustrated through 
the 11-bar indeterminate truss (11,9) shown in figure 4.2 with n fof = 11, m dof = 9, and r = 2. Its 11 DDR ((Pi — 
[B] r {X}) have the following explicit form: 


Pi' 
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The DDR for the truss are obtained first by writing the EE ([B]{F} = {P}) and then {p} = [B] r {X}. The nine 
displacements are labeled in figure 4.2. The 1 1 deformations {p} represent elongations due to force in the 
1 1 bars of the truss. Two compatibility conditions are obtained by eliminating the nine displacements from the 
eleven deformation-displacement relations given by equation (4.25) as 

p 1 -V2p 2 +p3-^/2p 4 +P5+P6=0 ( 4 - 26a ) 

p6 - V2P7 +Ps ->/2P 9 +P10 +P11 = 0 ( 4 . 26 b) 


NASA/TP— 2004-207430 


117 




By combining equations (4.26a) and (4.26b), we can obtain the (2 x 11) compatibility matrix [C] as 


[C] = 


1 -V2 1 1 1 

0 0 0 0 0 1 


0 0 0 0 0 
-V2 1 -J2 1 1 


(4.26c) 


For the truss (9, 1 1), the CC is a (2 x 1 1) matrix. Because the CC are homogeneous equations ([C]{|3} = 0), through 
linear combination of the rows of the [C], the following full matrix can be obtained 


\rl- 

1.000 

1.000 

1.000 

-1.417 

-1.417 

0.874 

-0.012 

-0.012 

0.018 

0.018 

-0.012 ' 

L C J - 

-0.012 

-0.012 

-0.012 

0.018 

0.018 

0.874 

1.000 

1.000 

-1.417 

-1.417 

1.000 


(4.27) 


The full matrix jcj can be used in IFM analysis. Theoretically speaking, there is the danger of generating a full 
compatibility matrix, destroying its banded nature. Several different procedures have been attempted to eliminate the 
deficiency associated with the generation of the homogeneous equations. The conclusion, however, is that a direct 
elimination of displacements from the DDR with row pivoting usually produces a well-banded matrix [C]. Elaborate 
schemes do not provide much superior results. The recommendation, therefore, is to follow direct elimination with 
adequate precision to reduce round-off errors, which can increase the bandwidth. With the fast, improved computa- 
tion facility that is currently available, a few additional entries in the [C\ matrix do not pose a problem. It should be 
pointed out that a direct assembly of the CC for a general finite element application, even when possible, would not 
be likely to substantially improve the compatibility conditions. 


Bandwidth of the Compatibility Conditions 

The compatibility conditions of finite element models are banded. On the basis of bandwidth considerations, the 
CC can be divided into three distinct categories: interface, cluster or field, and external CC. For the finite element 
model shown in figure 4.3, the three types of CC are illustrated in figure 4.4. 

Interface compatibility conditions . — Numerous interfaces internal to the structure are created in the finite ele- 
ment discretization process. The interface is the common boundary shared by two or more finite elements. In the 
figure 4.3 finite element model, the common boundary along nodes 1 and 7 is the interface between elements 1 
and 2, the boundary connecting nodes 12 and 17 is the interface between elements 13 and 14, and so on. Consider 
the interface between elements 1 and 2 as shown in figure 4.4(a). The deformations of elements 1 and 2 must be 
compatible along the common boundary defined by nodes 1 and 7, which gives rise to the interface CC. The number 
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(a) Interface CC. 



(b) Cluster or field CC, (c) External CC. 


Figure 4.4.- Bandwidth of compatibility conditions (CC). 


of CC at the interface depends on the element types (membrane, flexure, solid tetrahedron, etc.) and their numbers. 
The maximum bandwidth of the interface compatibility conditions ( MBW^ CC ) can be calculated as 

JT 

MBW icc =^{fof ej ) (4.28a) 

M 

where JT is the total number of elements present at the interface and fof ej represents the force degrees of freedom of 
the element j present at the interface. 

The bandwidth MBW represents the maximum bandwidths of the interface compatibility conditions written 
either in terms of forces {F}, as in [C][G]{F) = 0 (where we are referring to the bandwidth of the product matrix 
([C][G])), or in terms of deformations {(3), as in [CUP) = {0} (where the bandwidth is that of the compatibility 
matrix [C]). The actual bandwidth of the compatibility matrix [C] is smaller than its maximum bandwidth. 

The interface CC of discrete analysis are analogous to the BCC in elasticity (ref. 21). The interface CC are the 
most numerous compatibility conditions in any finite element model. These can be generated by writing the deforma- 
tion displacement relation of the local region (such as shown in fig. 4.4(a) for the interface defined by nodes 1 and 7) 
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and then eliminating the displacements from the local DDR, For the interface shown in figure 4.4(a), there are two 
elements (i.e., JT = 2). Let us assume that both are membrane elements; the /o/ of the triangular element fof t is 3 and 
that of the quadrilateral element fof is 5. Then, the MBW jcc calculated from equation (4.28a) is 8. 

Cluster, or field, compatibility conditions. — Consider any element in the model shown in figure 4.3, for exam- 
ple, element 19. This element, along with its eight neighboring elements, is also shown in figure 4.4(b). The deforma- 
tions of element 19 must be compatible with those of its neighboring elements (14 to 16, 18, and 20 to 23). For such 
a cluster of elements, the CC are referred to as cluster or field CC, which essentially represent St. Venant’s strain 
formulation in the field. The maximum bandwidth of the cluster CC can be calculated as 

JTC 

MBW CCC =Y J {f°fe j ) (4.28b) 

7=1 

where JTC is equal to the total number of elements present in the cluster. If we assume that there are five force 
degrees of freedom for quadrilateral elements and three for triangular elements, the bandwidth calculated from equa- 
tion (4.28b) is MBW CCC = 41. 

External compatibility conditions. — Reactions are induced at restrained nodes. If such restraints on the bound- 
ary exceed the number of rigid-body motions of the structure, then it is externally indeterminate. The degree of 
external indeterminacy R ext can be calculated as 

R ext= N x~ N f (4.28c) 

where N x is the number of displacement components suppressed on the boundary and A^-is the number of boundary 
conditions required only for the kinematic stability of the structure. 

Let us assume that the finite element model shown in figure 4.3 represents a membrane structure. Then, 

R ext = 7-3=4, since the number of actual boundary restraints is N x = 1 and the kinematic stability requirement 
is Nj-= 3. To calculate the bandwidth of the external CC, separate the local region connecting any two restrained 
boundary nodes. Let the number of elements between the two nodes be represented by JTE, then the maximum band- 
width of the external CC ( MBW ecc ) is given by 


JTE 

MBW ecc = ^(fof ej ) (4.28d) 

7=1 

If we assume, as before, that there are five and three force degrees of freedom for the quadrilateral and triangular 
elements, respectively, then the maximum bandwidth of the external CC for the boundary segment shown in figure 
4.4(c) is MBW ecc = 8. 

The interface, cluster, and external CC are the local constraints. All three categories of local conditions are con- 
catenated together to form the system, or global CC, of the structure («, m). The sum of the number of interface (r- ), 
cluster (r ccc ), and external (r ecc ) CC is equal to the r-n-m of the structure (n, m); that is, r = r icc + r ccc + r gcc . 

Our recommendation is to generate the CC by direct elimination of m displacements from the n DDR. Once the 
CC have been generated, their bandwidth can be checked against equations (4.28a), (4.28b), and (4.28d). 

The equilibrium equations and the compatibility conditions are linked through the DDR. We can obtain the EE 
from the DDR through variational concepts of advanced calculus, and a direct elimination of displacements from the 
DDR yields the CC. The null product property (see eq. (2.17)) of the two fundamental matrices ([S] and [C]) implies 
that error in the EE can propagate to the CC and vice versa. 


Illustrative Example 16: Fixed Bar 

Integrated Force Method. — The IFM analysis process can be illustrated through the fixed bar that was solved 
in chapter 3, Illustrative Example 1 . This bar is subjected to thermomechanical loads. Such a bar, along with its end 
conditions and analysis parameters, is depicted in figure 4.5(a). The total length of the bar is 3 L. It is idealized by 
three finite elements consisting of a central element and two boundary elements of equal length L. The cross- 
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(a) Fixed bar. 



(b) Discretized bar (3,2). 



Figure 4.5.- Analysis of a bar by the Integrated Force Method, showing (P-| and P 2 ), 
bar forces (Fj, F 2 , and F 3 ), nodal displacements (X n and X 2 ), and lengths of 
elements (€). 


sectional areas of the boundary elements A are equal, and the area of the central element is 2A. The bar is made of 
steel, its Young’s modulus is E, and its coefficient of thermal expansion is oc. It is subjected to mechanical loads P j 
and P, at one-third and two-thirds of its length. The temperature distribution of the central element is T 2 , and the 
temperatures of the boundary elements are equal to T v The IFM and IFMD solutions are illustrated for this problem. 

IFM solution for the fixed bar . — The IFM solution requires three matrices: [B], [G], and [C] (seeeq. (4.1)). 
The equilibrium matrix [B] and the flexibility matrix [G] are assembled from their elemental matrices. The compat- 
ibility matrix [C] is obtained by eliminating displacements from the deformation displacement relations. The genera- 
tion of the matrices for the problem, following the procedure given earlier in this chapter, is described next. 

Elemental matrices for a bar element: By following the procedure that was illustrated for the plate element in 
equations (4.7) to (4.20), we can obtain the element matrices for the bar. The equilibrium matrix is generated for the 
following displacement and stress fields: 


U — Wj + 


U 2 — U\ 



(4.29) 


<7 = —— (4.30) 

A 

where u x and u 2 are the axial nodal displacements for nodes 1 and 2, respectively, in the local system (x ( , y , ) as 
shown in figure 4.6. The uniform internal force in the bar is F, the associated stress is O, and the area of the bar is A. 
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d Xg 

Figure 4.7.- Nodal displacements for 
a truss element in a global system. 


The displacement in the field u varies linearly along the length of the bar. Equations (4.29) and (4.30) can be rewrit- 
ten using the interpolation functions given by equations (4.1 1), (4.12), and (4.13), as follows: 


{o} = [T]{F} = r|l{F} 

LA J 



[Z] = [£][*] = [“ j] (431) 

The elemental equilibrium matrix becomes 

[B e ] = \[Z] T [Y]dv 
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The equilibrium matrix for a truss element in local coordinates is a (2 x 1) matrix. Its rows correspond to the two 
displacements (u v u 2 ). Its single column corresponds to the internal force F. 

The equilibrium matrix can be written in global coordinates (x^, y g ) by expressing the displacement u in terms 
of the nodal displacement as shown in figure 4.7. In this figure, (Xj, X 2 ) and (X 3 , X 4 ) represent the global nodal dis- 
placements at node 1 and 2, respectively. The bar orientation (local x ( -axis) with respect to the global axes (v,,, y^) 
is 9. Let us define t and m as the direction cosines (for the bar AB) of the angles between line AB and the x^~ and 
y^-axis, respectively. The nodal displacements (u v « 2 ) can be written as 


{ M 1 l= 

" i 

m 

i — 

o 

o 

w 

0 

0 

e m\ 


*1 

X 2 

x 3 

Xa. 




(4.33) 


The displacement in terms of nodal displacements (Xj, X 2 , X 3 , X 4 ) can be written as 


or 


[Z] = [L][N][X] 


(4.34) 


The equilibrium matrix [b| ] in the global nodal displacement {X} can be written as 


KWKH" ¥ I 1 " 1* 


or [«j] = [ ir ][ B 1 <4 ' 35) 

(4x1) (4x2) (2x1) 

The equilibrium matrix for a bar element in global coordination is a (4 x 1 ) matrix that corresponds to the four global 
displacements (Xj, X 2 , X 3 , X 4 ), as follows: 



(4.36) 


The equilibrium matrix [BJ in IFM equation (4.1) is written in global coordinates. Thus, the matrix j given in 
equation (4.36) should be used to assemble the system equilibrium matrix [B], 

The flexibility matrix for the bar element is obtained by substitution in equation (4.18) as 

[G e ] = j[y] r [D][T]dv (4.37) 
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For the bar element, 
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The flexibility matrix is a (1 x 1) matrix because the bar has a single force unknown (F). 

For illustrative example 16, three bar elements are used to discretize the structure, as shown in figure 4.5(b). The 
nodes and forces are labeled to match the quantities in figure 3.1 . The force and displacement degrees of freedom of 
the structure shown in figure 4.5(b) are as follows: 

(1) Force degrees of freedom: Each bar element is idealized by one internal force, or the structure has three 
force degrees of freedom ( n fof= 3; F v F 0 , F 3 ). 

(2) Displacement degrees of freedom: The structure has two displacement degrees of freedom, one at each of its 
two free nodes, 2 and 3, ( m dof= 2; X v X 2 ). For the IFM analysis, the structure is designated as bar (3, 2). It has 

m = 2 equilibrium equations and r = n - m = 1 compatibility condition. The local and global systems for the problem 
are identical, or direction cosines l = 1 and m = 0; thus, the local coordinates can be used to solve the problem. 

Equilibrium equations. — The two-system EE of the bar are assembled from the three elemental EE (see eq. 
(4.32)): 

Bar 1 (refer to fig. 4.5(b)) 


Bar 2 


Bar 3 


dof 

i l^fof 


II 

'-1 ' 

o 

— 1 

— i 

1 


2 

-1 

1 



(4.39a) 


(4.39b) 


3 



'-T 

L j 2 

1 


(4.39c) 


The (2 x 3) system equilibrium matrix [5] is obtained by following the standard finite element assembly 
technique: 


[B] = 


1 

2 


1 

-1 

0 


2 3 <— fof 

1 O' 

-1 1 


t 

dof 


(4.39d) 
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and the two EE can be written as 


-1 

0 



(4.39e) 


The EE given by equation (4.39e) is identical to equation (3.1 .1c) from example 1 in chapter 2. 

Flexibility matrix . — The elemental flexibility matrix for the bar is a (1 x 1) matrix [G e ] = t/AE. The flexibility 
matrix [G] for the structure is obtained by concatenating the elemental matrices along the diagonal as follows: 


[G] = 


A,E 


AnE 


A,£ 


(4.40) 


Compatibility conditions. — The first step in obtaining the CC is to establish the DDR ({(3) = \B]^{X) )■ The 
three-component deformation vector {(3} corresponds to the three elemental expansions due to forces F j, Fy and Fy 
respectively. This deformation vector is related to the displacements (Xj, X 2 ) by the DDR ({|3} = \B\^{X\), which 
has the following form: 


Pi 

r 

P 2 

ii 

,P3. 

L 



(4.41a) 


The two displacements are eliminated from the three DDR to obtain one CC (r = n - m = 3 - 2 = 1). In terms of 
deformations, the CC, [C]{ |3} = {0), has the following explicit form: 


[1 1 1 ] 


Pi 

P 2 

P3 


= { 0 } 


(4.41b) 


The CC constrains the total elongation of the bars to zero (Pj + (3 2 + P 3 = 0). For this simple case, this could have 
been asserted also by observation. 

From equation (4.41b), the compatibility matrix [C] is obtained as 

[C] = [l 1 1] (4.41c) 

The correctness of the matrix [C] can be ascertained from its null property [E][C\ r = 0, here 


1 

-1 1 0 

1 

0 -1 1 


"O' 
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(4.41d) 
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Effective initial deformations . — Prescribed deformations due to other effects (here, thermal expansion) are 
accommodated in the effective initial deformation vector {8/?}. The effective initial deformation vector, which is 
obtained from initial deformations {p}° and the compatibility matrix [C], has the following form: 


p? 1 


V 


P2 

■ = al< 

T 2 

(4.42a) 

P 3 


Ji. 



The vector { 8/?} is obtained from the formula {S/?} = ~[C]{ j3}° as 


{SZ?} = -oc(2 r, + T 2 )L (4.42b) 

Consider the temperature distributions T ] = 7q/2 and T 2 = - T () . For this temperature distribution, the effective 
initial deformation vector {8P} = 0. Since { S/? } is zero, such compatible initial deformations do not induce forces in 
the structure. 

From the definition of matrices [Z?], [C], [G], and (SP), the final governing IFM equation ([S1{F) = {P'\) is 
assembled as 

1 -1 0 

0 1 -1 

1 A 1 /A 2 1 


pi 


-Pi 

V 2 

• = • 

-Pi 

[ f 3. 


E8R/L 


(4.43) 


The CC or the third equation in equation (4.43) is scaled by multiplying it with the factor Ell. The scaling process 
makes the governing IFM matrix [S'] dimensionless. Solution of equation (4.43) yields the forces from which dis- 
placements are calculated by back-substitution (from eq. (4.2)). 

Numerical results . — The parameters of the example are as follows: 

(1) Lengths of the bars: L, = L 9 = L 3 = 10 in. 

(2) Cross-sectional areas: A x = A 3 = 1 in. 2 and A 2 - 2 in. 2 

(3) Modulus of elasticity: E = 30 000 ksi 

(4) Coefficient of thermal expansion: 6xl0~ 6 per °F 

Case 1 — Mechanical loads only: The external loads are P, = 10 kips, P 2 = 20 kips, and 8 R = 0. 

Internal forces are 


Pi 


14 

P 2 

> — 

4 

F 3_ 


-16 


kips 

Nodal displacements can be calculated from equation (4.2) as 
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Case 2— Nonuniform temperature: The thermal load corresponds to an increase of temperature distribution for 
the central span, or Tj = 0 °F, T, = 2000 °F, and = 0 °F. 


Internal forces: 


Nodal displacements: 



Case 3 — Uniform temperature: In this case, Tj = T 2 = 77 = 2000 °F. 



h 


'432' 

Internal forces: 

f 2 


432 • 


A 


432 . 

Nodal displacements: j 

A] 

x 2 \ 

i j 0.024 
[ “ [-0.024 


Dual Integrated Force Method. — A dual formulation of the primal IFM has been developed. This formulation 
is termed the dual Integrated Force Method, or IFMD. The dual method is obtained from the IFM equations by map- 
ping forces into displacements at the element level. Like the IFM, the dual method has two sets of equations. The 
first (or primary set), which represents a symmetrical set of equation (4.5), is used to calculate the displacements 
(see eq. (4.4)). Forces (or stresses) are back-calculated from the secondary set of equations. IFM and IFMD are 
analytically equivalent — producing identical analysis results for stresses, displacements, frequencies, and buckling 
loads. The dual method does not utilize differentiation to calculate stresses. The primary equations of the dual 
method closely resemble the equations of the popular stiffness method. Significant differences and similarities 
between the dual and the stiffness methods are 

(1) The IFMD has two sets of formulas: one for the determination of displacements and another for the calcula- 
tion of forces or stresses. The stiffness method has one set of equations to calculate displacements. Stresses in the 
stiffness method are calculated by differentiating the approximate displacements to calculate strains, which can be a 
source of error, and then transforming strains to stresses. 

(2) The stiffness method utilizes the differentiation of approximate displacements to calculate stresses, which is 
avoided in the IFMD. 

(3) The symmetrical equations of the IFMD can use the solvers that have been developed for the stiffness 
method. In fact, a stiffness code can be changed into an IFMD code without substantial modification. 


The dual method developed and illustrated in this chapter is obtained from the four basic IFM equations: 


Equilibrium equations: 

[B]{F} = {P] 

(4.44a) 

Compatibility conditions: 

[C]{P} = {0} 

(4.44b) 

Flexibility relations: 

{P} e ={P}-{P}°=[G](F} 

(4.44c) 

Deformation displacement relations: 

{PMfifm 

(4.44d) 
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The IFM is obtained from these equations with forces as the primary unknowns: [5]{F} = {P}*. We can obtain 
the IFMD from these equations by considering the displacement {X} as the primary unknown in the following two 
steps. 

Step 1 — Eliminate deformations {[S} between the flexibility relations and the DDR to obtain the following force 
displacement relations: 


{F} = [G]- 1 [P] r {X}-[G]- 1 {P}° (4.45) 

Step 2 — Eliminate forces between the EE and the FDR to obtain the primary equation of the IFMD: 


[B][GT ] [B] T 

W 

= {p}+- 

[e][G]-'[p]°l 

(4.46a) 

_mxn nxn nxm _ 

mx 1 

mx 1 

mxn nxn nxl ] 



m ifmd {X} = {P} iJmd (4.46b) 

mxm mx ' rax 1 


where [ D ] iJme i =[ B ][ G ] l [ B Y 

Wifi*, ={^} + ([«][G]-‘{P}°)} 

The governing equation appears to be similar to the stiffness equation, with some variation in the load term. 

Once displacements are known from the solution to equation (4.46), the force displacement relation given by equa- 
tion (4.45) is used to back-calculate forces. Equation (4.46) is the primary equation, and equation (4.45) is the sec- 
ondary equation of the IFMD. 

Assembly of matrix [D]^ md . — The dual matrix [D]^ is assembled from elemental matrices in a process quite 
similar to the regular stiffness method. The assembly of the dual matrix [D]^ is illustrated here by considering the 
structure with three bars as an example (see fig. 4.5). Let the elemental equilibrium and flexibility matrices of the 
structure be represented by [B le ], \B 2e \ [B 3e ], and [G lc ], [G 2e ], [G 3e ], respectively. The assembled equilibrium ma- 
trix [B] of the structure can be written as 

[B] = [B le j B 2e ! B 3e ] (4.47) 

The flexibility matrix [G] is a concatenation of elemental matrices along the diagonal: 



[G] = 




(4.48) 


By using elemental equilibrium and flexibility matrices, we can define the following elemental pseudostiffness 
matrix: 


["W-'MIgTV] 


(4.49a) 
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Equivalent elemental loads are defined as 


< 449b » 

By substituting equations (4.47) and (4.48) into equation (4.46a) and expanding in terms of elemental matrices, we 
obtain 






[ G '1 

vn 

T 


={P}+- 

[b 2£ 


[g 2£ ' 




[b 3€ 


[g 3£ ] 

CO 

7 


We can write the system matrix of the IFMD as 

\P*\ifmd = \Pifmd-e\ \ ^ ifmd—el i ^ifmd—e3 ] 


(4.50) 


(4.51) 


{P} 


ifrnd 




{p} 

OT 


thermal_l 

thermal_2 

thermal_3 


The pseudostiffness matrix [D\ ifind for the finite element model can be assembled from the equivalent elemental 
matrices by following the standard stiffness assembly technique. The formulation of the mechanical nodal load 
vector {P} also follows the regular stiffness assembly procedure. 

The three-bar example used earlier to illustrate the IFM is used again to illustrate the IFMD solution process. 

Step 1 — Generate the elemental pseudostiffness matrices and load vectors: 

Equilibrium matrix: \b £ 



Inverse of flexibility matrix: 



Pseudostiffness matrix : 


[°w=M[ G fVf 



■ i 

-r 

l i J 

-i 

i 


Equivalent thermal load: 


{P} 


thermal -e 


-1 

1 


AE 


K1 
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Step 2 — Fonnulate the governing equation for the problem by using elemental matrices and vectors: 


Degrees of freedom (fig. 4.5(e)) — > 
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1 1 J 
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^ Element 1 ^ 
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l e J 
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Element 2 
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2 

v y 


Element 2 


2 0 
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'*i 


f-i 



>• 



1 1 j 

-i 
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*2. 



Element 3 

< > 



(4.52a) 


Equation (4.52a) can be simplified for the parameters of the problem (Aj = A 2 = 1-0 in. 2 and = 2.0 in. 2 ) to obtain 
the following equations of the dual method: 


Casel Case 2 Case 3 


' 3 -2 



10 

oc E 

'—21 AT 2 

ojE 

' i A7) —21 AT 2 ' 



* = 




■, — 






~r 


l 


—2 3_ 
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20 


2t AT 2 


-(A T 3 + 21 AT 2 


(4.52b) 


Solution of equation (4.52b) yields the following values for the displacements for the three load cases: 


J X l| _ f 0.46671 j-0.0481 J-0.0241 

Wj inch = l°.5333f { 0.048 J’ | 0.024 J 

We can calculate the forces from equation (4.45) as 


£ l] 


14' 


144' 


'432' 


> = *: 

4 

> f < 

144 

— < 

432 


kip 

— i6 
Casel 


144 

Case 2 


432 
Case 3 


(4.53) 


(4.54) 


The solution generated by the dual IFMD agrees with that of the primal IFM. Both the primal IFM and IFMD 
are based on identical sets of equations. These are the equilibrium equations, the compatibility conditions, the defor- 
mation force relations, and the deformation displacement relations. Only the solution order is changed. That is, in the 
primal IFM, forces are determined first and then displacements are back-calculated. In the IFMD, the reverse order is 
followed; that is, displacements are determined first and then forces are back-calculated. Therefore, the primal IFM 
and its dual IFMD yield identical solutions. The IFMD can be considered to be the true displacement method. The 
governing equations of the IFMD and the popular stiffness method are symmetrical, but the coefficients of the stiff- 
ness matrix and that of dual method pseudostiffness matrix can differ in magnitude. This aspect will be explained 
further in the next section. 
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Illustrative Example 17: Cantilevered Beam 

Integrated force Method. — The IFM solution procedure is further illustrated by considering two simple 
examples. The first is a cantilevered beam idealized by two membrane elements. For the example, all matrices and 
equations required for IFM and IFMD analysis are generated in closed form. Also, results from the stiffness method 
and MSC/NASTRAN codes are obtained and compared. The second cantilevered truss example also is solved in 
closed form. 

The cantilevered beam shown in figure 4.8 is 12-in. long, 2-in. deep, and 0.25-in. thick. It is made of steel with a 
Young’s modulus E of 30 000 ksi and a Poisson’s ratio v of 0.3. For static analysis, it is subjected to a tip load of 
10 lb at the free end. For dynamic analysis, the beam is considered massless with two lumped masses at its free end 
as shown in figure 4.8. Because closed-form analysis is used, a very simple four-node rectangular membrane element 
with five internal force unknowns, designated element IFMRC0405 (IFM rectangular four-node element with dof= 8 
and/o/= 5), is used. The problem is solved via the following methods: 

(1) Integrated Force Method 

(2) Dual Integrated Force Method 

(3) Regular stiffness method 

(4) MSC/NASTRAN stiffness method 


0.25 in. 


2 in> 

I 


Cross 

section 

z' 


Q | 0.00388 slugs (5-lb mass) 


O f 0.00388 slugs (5-lb mass) 


-12 in.- 


Figure 4.8.- Cantilevered beam idealized by two membrane elements (IFMRC0405). 


Generation of the element equilibrium and the flexibility matrices required by IFM/IFMD analysis are as 
follows. 

Elemental equilibrium matrix [B e ]: The elemental equilibrium matrix is obtained by discretizing the strain 
energy for the membrane, which can be written as 



,, du XT dv 

x dl + y ty + ^ 


^ du + dv 
y dy dx / 


dxdy 


where 

m, v membrane displacements 

N r , N... N r ,, membrane forces 

x, y Cartesian coordinates in domain S 

The rectangular element with spans 2 a and 2b and thickness t is shown in figure 4.9. 


(4.55) 
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Figure 4.9.- Rectangular membrane element IFMRC0405. 
Displacement interpolations for the element are as follows: 




(4.56b) 


where X v X,, . . . , are the eight displacement degrees of freedom of the element. The force interpolation selected 
has the following form: 


N X =F, 



Ny = +F4 — 

a 

N xy = F 5 (4.57) 

where Fp F 0 , . . ., F 5 are the five force degrees of freedom. 

The displacement function given by equation (4.56) represents a standard interpolation for a four-node mem- 
brane element. The membrane force interpolation given by equation (4.57) represents a constant shear (N ) and a 
linear variation for normal forces (N x , N y ) along the y- and ^-directions, respectively. The force interpolation will 
produce an acceptable stress distribution at the center of the element. 

Substitution of equations (4.56) and (4.57) into equation (4.14) and integration yields the following (8 x 5) 
nonsymmetrical, equilibrium matrix \B e ] for the element: 

b b/3 0 0 -a 

0 0 —a a 1 3 —b 

b -b/3 0 0 -a 

0 0 —a —a/ 3 b 

b b/3 0 0 a 

0 0 a a/3 b 

-b -b/3 0 0 a 

0 0 a -a / 3 —b 
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Elemental flexibility matrix [G e ]: The elemental flexibility matrix is obtained by discretizing the complementary 
strain energy, which has the following form: 



vAA 


(N v Y N 


V f i 



-2(l + v)^ 



2 


dxdy 


(4.59) 


Substitution of the force polynomials (N x , N y , NJ from equation (4.57) into the complementary energy expres- 
sion and integration yields the symmetrical flexibility matrix [G e ], as follows: 



(4.60) 


where E is the Young’s modulus and v is the Poisson’s ratio. 

System equilibrium equations: The equilibrium equations for the two-element cantilever beam are generated by 
following standard REC04 assembly procedure. The force and displacement degrees of freedom for the two elements 
are depicted in figure 4.10. The system equilibrium matrix [B] of dimension (8 x 10), which is assembled from the 
two elemental matrices is as follows: 


<- Element 1 «- Element 2 


'1 

-1/3 

0 

0 

-3 ! 

-i 

1/3 

0 

0 

-3 ‘ 

0 

0 

-3 

-1 

t 

i | 

0 

0 

-3 

1 

-1 

0 

0 

0 

0 

0 I 

i 

-1/3 

0 

0 

-3 

0 

0 

0 

0 

o ! 

I 

0 

0 

-3 

-1 

1 

0 

0 

0 

0 

0 1 
1 

1 

1/3 

0 

0 

3 

0 

0 

0 

0 

o ! 

0 

0 

3 

1 

1 

1 

1/3 

0 

0 

3 i 

-1 

-1/3 

0 

0 

3 

0 

0 

3 

1 

i * 

0 

0 

3 

-1 

-1 


(4.61) 


The number of entries in any column corresponds to the entries in one elemental equilibrium matrix irrespective of 
the problem size. Thus, the system equilibrium matrix is a very sparse matrix. 





*3 


>*2 


*1 



Figure 4.10.- Elemental equilibrium matrices. 
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Compatibility matrix: The compatibility matrix is obtained from the deformation displacement relations 
({(3} = [6] r {X}). In the DDR, 10 deformations (which correspond to the 10 force variables) are expressed in terms 
of eight displacements (X^, X,, . . Xy. The problem has two CC that are obtained by eliminating the eight displace- 
ments from the 10 DDR as 


[C] = 


0 0 -1/3 
0 0 2/3 


1 0 0 0 0 
0 0 0 0 -1/3 


0 0 
1 0 


(4.62) 


Like the [ B ] matrix, this (2 x 10) compatibility matrix [C] is also sparse. 

Flexibility matrix: The flexibility matrix for the problem is obtained by diagonal concatenation of the two 
elemental flexibility matrices as 


where 


[G] = 


[G], 


[C ] 2 


1 0 -0.3 0 0 

01/300 0 


[G]j =[G] 2 =1.6xl0 -6 


-0.3 


0 


1 


0 0 


0 001/3 0 


(4-63) 


0 0 0 0 2.6 


The governing equations, [5]{F} = {P}, of the IFM (which can be obtained from the equilibrium matrix [B], compat- 
ibility matrix [C], and flexibility matrix [G]) have the following form: 
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(4.64) 
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The first eight rows represent the EE, and the last two rows represent the CC. The last two CC rows are normalized 
by a factor of 10 -7 . Solution of equation (4.64) yields the forces: 


Fl 
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Fi 


135 

f 3 
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Fa 
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f 5 
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f 6 


0 

Fi 
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Fg 
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Fg 


0 
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-5 


From the forces, stresses can be generated from equation (4.57). For the second element, the stress components that 
are accurate at the center of the element are as follows: 


Oj = —^ = 540 psi 

a y = —j— = 0 psi (4.66) 

%Xy p=-20pSl 

The shear stress at t xy = -20 psi agrees with the average shear stress calculated from the strength of materials 
formula (t = P/A = 10/[2(0.25)] = 20 psi). There is no stress in the transverse direction (o y = 0), which agrees with 
the strength of materials results. 

In contrast to the strength of materials results of 720 psi, the normal stress is 540 psi. The discrepancy occurs 
because a (or N x in eq. (4.57)) is uniform along the length of the element with the origin at its centroid (see fig. 

4.1 1). At the elemental centroid, the strength of materials stress value at 540 psi is identical to that of the IFM 
solution. The element with five force unknowns given here to illustrate the IFM analysis procedure is a very simple 
element. This element should not be used for accurate stress analysis. Even for this element, the normal stress distri- 
bution can be improved by using a smaller element at the origin. For example when a smaller first element with a 
reduced length of 1 in. and an ll-in.-long second element are used to discretize the structure (as shown in fig. 4.12), 
the normal stress value improved to o x = 690 psi, which is 96-percent accurate. 

The displacements can be obtained via back-substitution in the formulas ({X} = \f\[G\{F\). The maximum 
value of the displacement, which occurs at node 5 or 6 along the y-direction, is 

8 max = -1-100 xl0~ 3 in. (4.67) 

Sbeamsolution = -1-152 X10- 3 in. (4.68) 

The displacements calculated with the two-element regular model are about 96-percent accurate, whereas the normal 
stress is 75-percent accurate. When the irregular model given in figure 4.12 was used, the stress was 96-percent accu- 
rate, whereas the accuracy achieved for displacement was only 86 percent with the maximum value for displacement 
(‘Vemlar = 0.919xl0' 3 in.). In other words, accurate displacement does not necessarily translate into a corresponding 
level of accuracy in the stresses. 
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Element 1 



Figure 4.11.- Normal stress distribution in 
pounds per square inch for element 2. 
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Figure 4.12.- Irregular model. 


Dual Integrated Force Method . — The pseudostiffness matrix of the dual method is obtained by using the two- 
element model given in figure 4.8 following the procedure for the fixed bar in example 17. The system equation of 
IFMD has the following form: 
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(4.69) 


Solution of equation (4.69) yields the nodal displacements as 




ifrnd 


■ 10 ” 


1.080 
3.344 
-1.440 
-11.008 
1.440 
-1 1.008 
1.080 
-3.344 


(4.70) 


The displacements obtained by IFMD are identical to those obtained by IFM, with the maximum value at 
l.lxlO -3 in., which is identical to that for the IFM solution. 

Stress calculation in IFMD . — In the dual method, the 10-component force vector {F}, which is back-calculated 
from the formula given by equation (4.42), {F} = [G] _1 [F] r {X}, is identical to that obtained for IFM (see eq. (4.65)). 
From forces, stresses that can be computed as indicated for the IFM become identical to IFM results. In brief, both 
IFM and IFMD yield identical results for stresses and displacements, as expected. 
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Regular stiffness method .— The regular stiffness matrix for the rectangular element is obtained by discretizing 
the strain energy (U ) given by equation (4.7 1) for the displacement functions (u, v) given by equation (4.56). 


U 


P 



du' 

3y, 


(i-v) [ r a« 1 3v"| 2 

2 ^ dy dx y 


dx dy dz 


(4.71) 


Using the standard stiffness procedure, one can obtain the equations for the two-element cantilever beam as 
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Solution of equation (4.72) yields the displacement vector as follows: 


(4.72) 


m* 


stiffness 


= 10 


-4 


-0.237 

-0.815 

-0.316 

-2.576 

0.316 

-2.576 

0.237 

-0.815 


(4.73) 


The maximum value of displacement (8 max = -2.576 X 10" 4 ) is much smaller, and it is only 22 percent of that of the 
strength of materials solution. 

For the irregular model shown in figure 4.12, the following displacement solution is obtained: 


—0.1 93 
-0.1 14 
-0.287 


Mi 


irregular 


:10 


-4 


-2.944 
0.287 
-2.944 
0.193 
-0.1 14 


(4.74) 


For the irregular model, maximum displacement (8 max — 2.944) is slightly better than the regular model and it is 
25.6 percent of that of the strength of materials solution. Stresses obtained from the stiffness method for the two 
crude elements are very' inaccurate. 
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MSC/NASTRAN stiffness solution . — The MSC/NASTRAN solution to the cantilevered beam problem is 
obtained with the CQUAD-4 element. This popular element is an advanced element in comparison to the simple 
IFM element RC0405. The final MSC/NASTRAN stiffness equation for the problem has the following form: 


0.0 

1.2 

-0.1 

-2.6 

-1.3 

-3.4 

0.0" 



0 ' 

16.9 

0.1 

3.8 

-1.3 

-4.3 

0.0 

-16.0 

*2 


0 


3.0 

-1.3 

-1.7 

-0.1 

-2.6 

1.3 

*3 


0 



8.4 

0.1 

-8.0 

1.3 

-4.3 

*4 


0 




3.0 

1.3 

1.2 

-0.1 

'*5 

* - ‘ 

0 

sym 




8.4 

0.1 

3.8 

*6 


-5 






6.1 

0.0 

*7 


0 







16.9 

* 8 . 


-5 


MSC/NASTRAN displacement for the problem obtained from solution of its equation is as follows: 


(4.75) 


-0.983 

-3.052 

-1.310 




MSC/NASTRAN 


- 10 " 


—10.036 

1.310 

-10.036 

0.983 

-3.052 


(4.76) 


The maximum stress generated by MSC/NASTRAN. which occurs for the element 1 centroid, is a = 540 psi. The 
maximum value for displacement (8 max = -0.10036xl0 _? in.) is 87 percent of the strength of materials solution. 

The MSC/NASTRAN solutions for the two-element irregular model shown in figure 4.12 are as follows. The 
displacement solution is 


-0.209 

- 0.122 

-1.310 


m, 


irregular 


= 10 " 


-8.671 

1.310 

-8.671 

0.209 

- 0.122 


(4.77) 


The maximum displacement is 72.8 percent of the strength of materials result. Stress for the model is 97.9 percent 
accurate, which is identical to that of the IFM/IFMD solution. 

IFM frequency analysis . — The IFM frequency equation j[S] - to 2 [([M] [7] [G]) / 0]Jf = 0 for the cantilever 
beam depicted in figure 4.8 has the following explicit form: 
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The IFM eigenvalue problem for the beam is an unsymmetrical (10 x 10) set of equations. The four nonzero 
rows in the mass matrix represent the participation of the lumped masses located at nodes 3 and 4. Participation of all 
10 force degrees of freedom is essential for the determination of correct frequency and force mode shapes (refs. 6, 9, 
and 29). A correct dynamic analysis formulation cannot be obtained through any elaborate manipulation of the equa- 
tions of the classical force method (ref. 20), because this approach cannot consider the mass or the inertia for redun- 
dant members. The four frequencies obtained for the problem are given in table II. Table III shows the fundamental 
force mode shape along with the displacement mode shape calculated by back-substitution with equation (4.2). 


TABLE H— FREQUENCIES FOR A CANTILEVERED BEAM WITH LUMPED MASS 


Frequency 

numbers 

Frequencies, Hz 

Analytical 

solution 

IFM 

IFMD 

Stiffness 

method 


Kg||B| 

168.32 

170.629 

170.629 

353.073 

178.707 

■ 

1052.00 

2031.261 

2031.261 

2032.188 

2031.922 


2950.38 

2534.278 

2534.278 

4874.846 

2647.656 


5786.00 

9298.475 

9298.475 

9599.293 

9510.976 


TABLE III. — FORCE-AND DISPLACEMENT-MODE SHAPES FOR FUNDAMENTAL 
FREQUENCY OF CANTILEVERED BEAM 


IFM 

IFMD 

Stiffness 

method 3 

MSC/NASTRAN 
QUAD-4 element 

Force mode 

Displacement 

Force mode 

Displacement 

Displacement 

Displacement 

shape 

mode shape 

shape 

mode shape 

mode shape 

mode shape 

0.0 

9.78xl0' 2 

0.0 

9.78xl0' 2 

9.17xl0' 2 

9.8xl0' 2 

1.0 

.303 

1.0 

0.303 

0.315 

.303 

0.0 

.131 

0.0 

.131 

.123 

.131 

0.0 

1.0 

0.0 

1.0 

1.0 

1.000 

-3.65xl0~ 2 

-.131 

-3.65x10'" 

-.131 

-.123 

-.131 

0.0 

1.0 

0.0 

1.0 

1.0 

1.000 

0.343 

-9.78xl0“ 2 

0.343 

-9.78x1 O' 2 

-9.18xl0' 2 

-9.78x1 0' 2 

0.0 

.303 

0.0 

.303 

.315 

.303 

0.0 


0.0 




-3 .65x1 O' 2 


-3. 65x1 O' 2 









a Force mode shape could not be determined. 


IFMD frequency analysis . — The IFMD frequency equation, neglecting damping, ([D]-co 2 [M])x = 0 for the 
two-element cantilevered beam depicted in figure 4.5, has the following explicit form: 



O 

o 

L3 

-0.1 

-2.6 -1.3 

-3.4 

0.0" 


"0 



16.6 

0.1 

4.0 

-1.3 -4.5 

0.0 

-15.6 


0 




3.0 

-1.3 

-1.7 -0.1 

-2.6 

1.3 


4 


10 6 



8.3 

P 

1 

^4 

oo 

1.3 

-4.5 

-co 2 xl0“ 3 

4 






3.0 1.3 

1.3 

-0.1 


4 



sym 



8.3 

0.1 

4.0 


4 







6.1 

0.0 


0 








16.6 


0 



There are eight equations for the dual method. The four nonzero entries in the diagonal mass matrix in equation 
(4.79) correspond to the lumped masses. The structure of the IFMD eigenvalue equation appears similar to standard 
stiffness equations. (The differences between IFMD and the stiffness equations are examined later.) The frequencies 
obtained by IFMD are given in table II. The fundamental displacement mode shape along with the force mode back- 
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calculated from equation (4.5) are shown in table IE. 

The IFM equations are unsymmetrical, whereas those for IFMD are symmetrical. Either set of equations can be 
used to obtain frequency, stress, and displacement mode shapes. The symmetrical form can use popular eigenvalue 
solution routines (such as DSPGV), which are readily available in the LAPACK public domain library (ref. 36). For 
unsymmetrical eigenvalue analysis, the LAPACK routine DGEGV is used (ref. 36). For static analysis, Harwell 
library routines (MA28AD, MA28CD, MA29BD, MA29CD, MA47AD, MA47BD, MA47CD, MA47ID) are used 
(ref. 40). However, the unsymmetrical IFM version can be more useful than IFMD in design and sensitivity analysis 
(ref. 22). Furthermore, the IFM solution that provides r = n-m number of zero frequencies and associated eigenvec- 
tors corresponding to the r number of compatibility conditions can be used to verify solution accuracy 

Stiffness method frequency analysis . — The stiffness method frequency equations ([AT]- © 2 [A7 ])a =0 for the 
cantilever beam depicted in figure 4.8 have the following explicit form: 



'7.6 0.0 

0.5 

-0.1 

-1.9 —1.3 

-4.8 

0.0"' 


‘0 



17.1 

0.1 

3.8 

-1.4 -4.3 

0.0 

-16.1 


0 




3.8 

-1.3 

O 

1 

N 

i 

-1.9 

1.3 


4 





8.6 

p 

1 

00 

L3 

-4.3 


4 


10 6 




3.8 1.3 

0.5 

-0.1 

-co 2 xl0~ 3 

4 






8.6 

0.1 

3.8 


4 




sym 



7.6 

0.0 


0 








17.1 


0 



W = o 


(4.80) 


Frequency and displacement modes obtained are presented in tables II and III. Stress mode shapes cannot be calcu- 
lated readily in the stiffness method. 

MSC/NASTRAN frequency analysis. — MSC/NASTRAN equations for a two-element beam with a CQUAD4 
element and appropriate condensation have the following form: 



"6.2 0.0 

1.2 

-0.1 

-2.6 

-1.3 

-3.4 

0.0" 


"0 



17.0 

0.1 

3.9 

-1.3 

-4.4 

0.0 

-16.0 


0 




3.1 

-1.3 

-1.7 

-0.1 

-2.6 

1.3 


3.9 





8.5 

0.1 

-8.0 

1.3 

-4.4 


3.9 


10 6 




3.1 

1.3 

1.2 

-0.1 

-co 2 xl0~ 3 

3.9 



sym 




8.5 

0.1 

3.9 


3.9 








6.2 

0.0 


0 









17.0_ 


0 _ 



W = o 


(4.81) 


Tables E and IE give the MSC/NASTRAN results for the problem. 

This discussion is centered around the fundamental frequency because other frequencies are expected to be 
inaccurate and thus observed (see table E) because of the two-element crude model. IFM and IFMD yield identical 
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frequencies at 170.6 Hz, which is about 1 percent higher than the beam solution of 168.3 Hz. The stiffness method 
yields a very high frequency at 353 Hz or about twice that of the beam solution. The MSC/NASTRAN QUAD-4 
element produces a frequency of 178.7 Hz, which is 6.2 percent higher than that produced by the beam solution. 

Only IFM and IFMD provide force mode shapes that are identical, as expected. Displacement modes generated 
by all four methods (IFM, IFMD, stiffness, and MSC/NASTRAN) agree quite well. 


Illustrative Example 18: Single-Bay Truss 

A square aluminum truss with a length of 20 in. and a Young’s modulus E of 10x10^ psi is depicted in figure 
4.13. The areas of both diagonal bars are equal to 4l / 2 in. 2 ; all other bars have an area of 1 in . 2 The solution is 
required for two load cases. Load case 1 consists of a mechanical load of 1000 lb at node 1 along the v-direction as 
shown in figure 4.13. Load case 2 consists of an elevated temperature (AT = 100 °C) for member 3, which has a coef- 
ficient of thermal expansion a of 6X10 -6 per °F. The problem was solved separately for mechanical and thermal 
loads (see Illustrative Example 1 1). A solution for both thermal and mechanical loads can be obtained by 
superpositioning the two solutions. 

Solution for mechanical loads: Both mechanical and thermal loads require the same EE, which are assembled 
from the six elemental matrices to obtain four EE ([fi]{F) = (PI). The equilibrium matrix ([B] ; ) for the ith bar 
element with internal force (F-) and direction cosines ((, m) can be written (see eq. 4.36) as 


[B], = 



(4.82) 


The EE for the six bars are as follows 


dof l 


0 

-1 ' 
0 
1 

© © — 
II 

Cl 

tc 

-1/V2 

-1/V2 

1/V2 

1 

F 2 [B] 3 = 3 

' 0 ' 
1 
0 

2 

0 

t=\ 

2 

1/V2 

l=m= 

4 

-1 

e=o 


m = 0 1/V2 


[B] 4 


— 1 /V 2 
1 /V 2 
1 /V 2 
- 1 /V 2 

l=\ld2 

m=— 1/V2 


0 

-1 

0 

" 0 ' 

0 

0 

r , 0 

1 

Ms = 3 

1 

F 5 [b ) 6 = () 

0 

4 

0 

0 

-1 


e=\ e=o~ 

m = 0 m =- 1 


(4.83) 


The assembly of the elemental EE matrices along the displacements (Xj, X 2 , X 3 , Xf) yields the following EE: 
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|- 20 in. «-| 

Figure 4.1 3.- Single-bay truss. 


lh 

1 l/fi 0 0 0 ol f 2 

0 1 / V2 1 0 0 0 F 3 

0 0 0 1/J2 1 0 F 4 

0 0-1 - 1 /S o oj 

k 

The sixth column is null because the sixth bar (F 6 ) is connected to fully restrained nodes 3 and 4, as shown in figure 
4.13. In IFM, null columns are accepted in the equilibrium equations of indeterminate structures without any adverse 
effects. 

Compatibility conditions: The six forces are expressed in terms of four EE, thus there are two CC. The CC are 
obtained first by writing the six DDR and then by eliminating the four displacements {X} from the DDR. The DDR 
((P) = [£]{X}) for the problem can be written as 


Pl=Xl 

X 3 -X 4 

P4 ~ J2 


O *l+*2 

P2 “ V2 

P 5 =*3 


p 3 = x 2 -x 4 

P 6 =0 

(4.85) 


Elimination of the four displacements (X V X 2 , X 4 ) from the six DDR yields the two compatibility conditions in 

six deformations: 

Pi 
P2 
P 3 

P4 

P 5 
P 6 


0 0 0 0 0 1 
1 - 2 V 2 1 - 2^2 1 0 
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The second CC in equation (4.83), (that is, P 6 = 0) corresponds to the sixth element, with the null sixth column in the 
EE given by equation (4.80). 

In the absence of temperature effect, the six deformations are related to member forces through the flexibility 
coefficients as follows: 


Pt = 



20F] 

E 






P 6 =-^ (4.87a) 
E 


The flexibility matrix [G] from the FDR (P) = [G]{F} can be formed from the diagonal concatenation of the flexibil- 
ity coefficients as 


[20 


40 


[G] = 


20 


40 


(4.87b) 


20 


20 


Thermal effect: For the problem, temperature in member 3 is included in the deformation for that member as 

P3 = P3 +P 3 


P 3 = — — (as before) 


P° =oc7T=12xlO -3 and P? = P^ = p^ = P° = P° = 0 


(4.88) 


The effective initial deformation vector {5 R} for the temperature in member 3 can be calculated as 

{5«} = -[C]{P) 0 4_ 12x ° 1() _, 


(4.89) 


The CC, [C][G]{F) = {8/?}, are obtained in terms of forces by eliminating deformations between the CC in deforma- 
tion and the FDR as 


20 

E 


0 0 0 0 0 1 
1 - 2 pr 1 -ijr 1 0 


Thermal 

0 


-12x10 


-3 


(4.90) 
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IFM governing equations ([S]{F} = IP}*): The EFM governing equation is obtained by coupling the EE and 
the CC: 


Mechanical Thermal 
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1/V2 
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o’ 

Fi 
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0 ' 
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1/V2 

1 

0 

0 

0 

F2 
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0 

0 

0 

0 

1/V2 
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0 
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0 

0 

0 

-1 

-1/V2 

0 

0 

Fa 
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0 

0 

0 

0 

0 

0 

1 

f 5 
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0 

1 

-2 jr 

1 

-2V2 

1 

0 

F (i 


0 


-6xl0 3 


Solving the IFM equations yields the forces as 


Mechanical Thermal 


F x 


-545.5' 


-545.45' 

f 2 


771.4 


771.39 

F 3 


454.5 


-545.45 


* — ' 


and 


>4 


-642.8 


771.39 

F 5 


454.5 


-545.45 

Ff, 



lb 

0 


Displacements back-calculated from the IFM formula ({X} = [T][G]{F} + {(3}°) are 


Mechanical Thermal 


*1 


-1.091 


-1.091 

x 2 


5.454 


5.454 


• — • 


•xlO J and 

-1.090 

X3 


0.909 


,^4 

in. 

4.545 


-5.454 


(4.91) 


(4.92) 


(4.93) 
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Appendix A 

Classification of Variables and Methods 
of Structural Mechanics 


Structural mechanics deal with three types of variables: (1) force F, (2) displacement A, and (3) deformation p. 
These variables are related through four types of relations: (1) equilibrium equations (EE), (2) compatibility condi- 
tions (CC), (3) force deformation relations (FDR), and (4) deformation displacement relations (DDR). The choice of 
primary variables and the requisite relations form the different structural mechanics methods (see table I): (1) the 
Integrated Force Method, (2) the displacement method, (3) hybrid method, and (4) total formulation. This appendix 
provides a brief description of the variables and their relationships. The methods of structural mechanics can be 
developed with the matrices of the equilibrium equations, compatibility conditions, and force deformation relations. 

Consider the equilibrium equations, which relate internal force {F} to load {P} through the equilibrium matrix 
[B], as 


[B]{F} = {P} (A.la) 

Force {F} is considered to be the primal variable of the EE. In this report, the EE is obtained as the vectorial summa- 
tion of forces {F} and loads {P}, see equation (1.2) as an example. Alternatively, it can be obtained as the stationary 
condition of a potential function 7t e that can be defined as 

n e ={X} T [B]{F}-{X} T {P} (A.lb) 

The stationary condition of n e with respect to a certain set of variables yields the EE. The “set of variables” is re- 
ferred to as the dual variables of the equilibrium equations. Displacements {A} are the dual variables of the equilib- 
rium equations, or 


5 {X} (7t ( ’)=^[B]{F}-{P} = 0 (A.lc) 

Forces {F} are the primal variables of the EE, whose dual variables are the displacements {A}. The concept 
illustrated for the equilibrium equation can be extended to other relations. For example, deformations { p } are the 
primal variables of the compatibility conditions ([C](P} = 0), and redundants {R} are its dual variables. Likewise, 
displacements {A} are the primal variables of the deformation displacement relations ([5] 7 '{ A} = { P }), whose dual 
variables are the forces {F}. For the force deformation relations ([G]{FJ = { P }), forces {FJ represent both the pri- 
mal and dual variables. If the primal and dual variables are identical, such as for the FDR, then only its coefficient 
matrix [G] becomes a symmetrical matrix. 

The IFM with force {F} as the unknown is the force method with the governing equation 

m{F} = {P) (A.2a) 

The IFMD with displacement {A} as the unknown becomes the displacement method with the governing equation 

[D]{A] = {P} (A.2b) 
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The governing equations of the hybrid method with force I F } and displacement {X} as the unknowns can be written 
as 


s 

0~ 

F 


p 

-JG 

/ 

X 


0 


(A.2c) 


The coefficient matrix of the hybrid method is essentially uncoupled because forces can be determined first; then, 
the displacements can be recovered by back-substitution as {X} = [/|[G]{F}. 

The governing equations of the total formulation with force {E}, displacement {X}, and deformation {p} as the 
unknowns can be written as 


s 

0 

o' 

F 


p 

-JG 

I 

0 

X 


0 

-G 

0 

I 

Ip, 


0 


(A.2d) 


The coefficient matrix of the total formulation also can be considered to be an uncoupled matrix. All four methods 
given by equation (A.2) should provide the same solution fidelity if the key EE matrix [B] and the FDR matrix [G] 
are generated correctly. From known forces {B}, displacements {X} can be back-calculated. Likewise, from dis- 
placements, forces can also be recovered by back-substitution (see eq. (4.2)). Thus, traditionally, the force method, 
or the IFM, and the displacement method, or the IFMD, are considered to be the two principal structural mechanics 
formulations. Table I summarizes the analysis methods, their primal variables, and their associated variational 
functionals. 
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Appendix B 

Solution Through Superposition 
and Redundant Techniques 

Simple indeterminate problems can be analyzed through superposition and redundant techniques. Both tech- 
niques bypass the explicit use of the compatibility conditions. In these techniques, the problem is separated into 
determinate subproblems that can be solved without the compatibility conditions. The solution to an indeterminate 
problem is generated by adding the subproblem solutions. The redundant and the superposition techniques are illus- 
trated by considering a beam as an example. A beam of length ^ is fixed at both ends and is subjected to a uniformly 
distributed load of intensity q per unit length. The IFM solution for the problem was given earlier under Illustrative 
Example 6. 


Solution Through the Superposition Method 

A uniform beam with constant rigidity El is shown in figure B . 1 (a). The free-body diagram of this beam is 
shown in figure B.l(b). Because of symmetry, the reactions are equal (R A = R B ) and the transverse equilibrium equa- 
tion yields their values: 


Ra=Rb = R = y < b1 ) 

Likewise, because of symmetry, the fixed-end moments at support A and B are equal (M A =M B ), as shown in figure 
B.l(b): 

M a =M b =M (B.2) 

Analysis of the problem requires the determination of the single moment M. The superposition technique 
replaces the real structure with two determinate structures, which are shown in figures B.l(c) and B.l(d). The first 
structure, shown in figure B.l(c), is a determinate structure subjected to the given external load q that produces a 
symmetrical rotation 0j at both supports and B. The second structure, shown in figure B.l(d), is the same determi- 
nate structure, but it is subjected to moment M r producing the symmetrical deformations 0 2 . 

The key idea of the superposition technique is to generate a solution to the original indeterminate problem by 
adding the two determinate solutions. For the problem, the superposition technique constrains the determinate slopes 
(0 ] + 0, = 0) to zero to reinstate the fixed-end boundary conditions (0 = 0). The superposition method has two steps. 
In the first, solutions are obtained for the two determinate subproblems shown in figures B.l(c) and B.l(d). In the 
second step, the two determinate solutions are added to obtain the solution for the indeterminate problem. 

Step 1— Determinate solutions: The two determinate problems depicted in figures B.l(c) and B.l(d) can be 
solved from equilibrium considerations only. 

Solution to first determinate subproblem . — For the problem depicted in figure B. 1(c), the bending moment at a 
location x can be calculated as 


qlx qx* 




(B.3) 
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(b) Free-body diagram. 
q per unit length 



(c) Subproblem 1- simply-supported beam subjected to given load. 



(d) Subproblem 2- beam subjected to end moment M. 



(e) Bending moment diagram. 


Figure B.1 .- Solution by superposition methods. 
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The moment is integrated to obtain the displacement function w by using the moment curvature relation as 

d~W\ 1 f qlx qx 2 'j 


dx 2 El El 2 2 


Upon integration, the displacement function is obtained: 


1 f qlx 3 qx 4 

W i = hCiX+C? 

1 El 12 24 1 7 


The two constants of integration (c ]5 c 2 ) are calculated from the boundary conditions (see fig. B.l(c)) as 

wi = 0 at x — 0 or c 2 - 0 

ql 3 

wi - 0 at x = £ or Ci = 

1 1 24 

and the displacement function has the following form: 


(x) (2 £x 2 -x 3 -( 3 ) 


The slope or rotation at A can be obtained as 


Solution to second determinate subproblem. — Likewise, the second subproblem depicted in figure B.l(d) is 
solved by calculation of the bending moment at a location x: 

m 2 (x) = M r (E 

The integration of the moment curvature relation yields the displacement function w as 

d 2 W 2 ^ „ M r j 

dx 2 El El 


M r x ~ j 

= \-a-\ x + d'y 

2 El 

The two constants of integration (d^, d 2 ) are calculated from the boundary conditions (see fig. B.l(d)): 

w 2 = 0 at x = 0 or d 2 = 0 


W 2=0 at x = ( or d\=— - 


The displacement function has the following form: 


W2= ^S_( x 2_ £x \ 

2EI [ 1 
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The slope or rotation at A can be obtained as 


dw _ M r l 
dx x= q 2 El 


(B.13) 


Step 2 — Superposition of determinate solutions: The bending moment M r is the only unknown of the beam, as 
shown under subproblem 2 in figure B.l(d). The superposition principle adjusts the bending moment M r for subprob- 
lem 2 until the induced rotation 9, matches the rotation 0j in subproblem 1, producing the fixed-end condition with 
no rotation (0 = 0) for the original” problem. The superposition principle can be written mathematically as 


0 1 + 0 2 = 0=0 


<?l 3 M T t = Q 
2AEI 2EI 


M r =-^— = M 

r 12 


(B.14) 


The moment M obtained here is identical to the IFM solution given by eq. (3.6.15). Likewise, the reaction R 
given by eq. (B.l) agrees with the shear force given by eq. (3.6.15). The bending moment diagram for the problem, 
which is obtained by superposing the two subproblem solutions, is shown in figure B.l(e). 

The displacement function for the problem is obtained by adding the two responses for the two subproblems: 


w = Wj + W 2 for M r — — 


El 


x 3 e 2 x 2 x 4 ^ 

12 24 24 


(B.15) 


The displacement functions given by equations (B.15) and (3.6.19) for the IFM solution agree. 


Solution Through the Redundant Technique 

The problem is one-degree indeterminate because of symmetry (see fig. B.l(b) and eqs. (B.l and B.2)), or it has 
a single redundant force as the unknown (r = 1). In the redundant technique, a determinate basis structure is gener- 
ated from the original structure by introducing as many virtual “cuts” as there are redundants. For a structure with 
a single redundant, one cut is introduced at the center span to obtain the basis structure, as shown in figure B.2(a). 
Normally, there should be a moment and a shear force at the cuts of the beam. From symmetry, the shear force is 
zero and the cut has only a moment M, which is considered to be the unknown redundant moment. The redundant 
solution method has two principal steps: (1) determination of the redundants, then (2) solution of the basis structure 
subjected to the external loads and the redundants. 

Step 1— Determine the redundant moment: To determine the redundant moment M, one needs to solve the 
basis determinate structure twice — once for the external load, as shown in figure B.2(a), and again for the redundant 
moment M, shown in figure B.2(b). The two solutions are added to restore the original continuity condition at the 
cut. For this problem, because of symmetry, the continuity condition becomes zero (0 = 0) at the beam center. 
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II 1 ■ 

I ' 

^VtV,h> 

(a) Basis determinant structure with a cut at the center. 



(b) Basis structure subjected to redundant moment M. 



(c) Determinate structure subjected to external load q and redundant moment M. 



(d) Bending moment diagram. 

Figure B.2.— Solution by the redundant force method. 
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(B.l 6) 


This conditioni s enforced as(see figs. B.2(a) and B.2(b)) 

e c +e D = 6 = 0 

where 6 C is the rotation for the external loads (fig. B.2(a)), Q D is the rotation due to the redundant moment M 
(fig. B.2(b)), 0 is the rotation at the center of the real structure (6 = 0). 

Calculation of rotation 9 C : Determination of the rotation 0 C requires the integration of the moment curvature 
relation. The momenta V) from figure B.2(a) can be calculateda s 
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The integration constants are determined from the boundary condition at the origin: 


The rotation at cut C is obtained as 
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(B.l 7) 


(B.18) 
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(B.19) 


Calculation of rotation Q D : Determination of the rotation 0 D requires the integration of the moment curvature 
relation, which from figure B.2(b) can be written as 

d 2 w _m r _M 
dx 2 ~ El ~ El 


Mx 2 

W — H C\ X + Cl 

2EI 1 


(B.20) 
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The integration constants are determined from the boundary condition at the origin: 


w = 0 at x = 0 or Cj = 0 

^ = 0 at x = 0 or Cl =0 (B-21) 

dx 


w = 


Mx 2 
2 El 


and the rotation at cut C is obtained as 


0 


Ml 
D ~ 2EI 


The slope continuity condition (0 C + 0 D = 0) yields 


(B.22) 
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(B.23) 


Step 2— Solution of the basis structure: The solution of the basis structure subjected to the external load q and 
the redundant moment M yields the solution to the original problem. Because of symmetry, only half of the basis 
structure is considered, as shown in figure B.2(c). 

The moment is obtained by adding Dt 1 and "7K r (see fig. B .2(c)) as 


■7K 1 +’M r 


ql 2 qlx 
~12 2 



(B.24) 


The moment diagram for the beam obtained from the determinate structure shown in figure B.2(c) is depicted in 
figure B.2(d). Careful examination will confirm that the moment diagrams given in figure B.l(e), which were 
obtained from the superposition solution, and the diagrams given in figure B.2(d), which were obtained from the 
redundant solution, are equivalent. 

Calculation of displacement: The displacement function can be determined by integrating the moment curvature 
relationship (see fig. B.2(d)) as 
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The boundary conditions (w = 0 and dw/dx = 0) at A (x = 0) yield the value of the constants (c, =c 2 = 0). The 
displacement function becomes 



( 


V 


^ 2 x 2 


24 




(B.26) 


Both the superposition solution and the redundant solution solved this simple problem. Extension of the methods 
for the general solution of problems becomes cumbersome, especially when computers are used. Both methods are 
difficult to use in the dynamic analysis of structures. Even though the methods were favored during the manual calcu- 
lation era, they are currently out of favor. 
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Appendix C 

Standard Strength of Materials Formulas 


Elementary structural analysis typically provides the solution for axial force, shear force, bending moment, and 
torque. The stress and strain induced in a beam because of these forces can be calculated through the use of standard 
strength of materials formulas. These formulas, which are derived in standard strength of materials textbooks, are 
summarized here. 


Axial Force 


The stress o due to a force F in a uniform bar with cross-sectional area A can be determined as 


F 

a- — 
A 


(C.l) 


The strain e due to the stress a can be calculated from Hooke’s law of the material and the Young’s modulus of 
elasticity E as 


o 

e = — 
E 

For an axial member we can write 


(C.2) 


a -Ee = — (C.3) 

A 

The stress and strain distribution for a uniform bar of depth d, thickness b, and area A (A = b d) is illustrated in 
figure C.l. 


d 





(a) Cross section 


(b) Stress distribution. 


Figure C.l .- Stress and strain due to axial force. 
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Shear Force 


The shear stress T in a uniform beam due to a shear force V can be determined as 


x = 


VQ 

It 


^ max 



(C.4) 


where 
A area 

b thickness of the beam 
I moment of inertia 

Q first moment of the area with respect to the neutral axis 

In the shear stress formula given by equation (C.4), the shear stress has a parabolic distribution along the depth 
about the neutral axis (or the x-axis in fig. C.2) because the first moment of the area is a function of the y-axis. The 
shear stress has the maximum value x max =1.5 {VIA) at the neutral axis. 

The shear strain y due to the shear stress x can be calculated from Hooke’s law of the material and the shear 
modulus G, where G = E/[2(l+ v)], as 


Y = § (C.5) 

The shear stress and shear strain distribution for a rectangular uniform bar of depth d, thickness b, and area A, where 
A = bd, is illustrated in figure C.2. 


Bending Moment 

The bending stress o in a uniform beam that is oriented along the x-axis has a linear distribution along the 
y-axis, as shown in figure C.3. The bending stress at the distance y from the neutral axis can be calculated from the 
following formula: 


My 
o = — — 
/ 



(C.6) 


where a is the stress along the beam depth at a location y from the neutral axis and e is the strain associated with the 
stress o. 

The bending stress is zero at the neutral (orx) axis, and it peaks at the extreme fibers of the beam as shown in 
figure C.3. 


Torque 

The shear stress x in a uniform shaft that is oriented along the x-axis has a linear distribution along the y-axis, as 
shown in figure C.4. The shear stress x and strain y at the distance F from the neutral axis can be calculated from the 
following formula: 


x = 


77 

J 



(C.7) 


where J (J = 0.5 tiR 4 ) is the polar moment of inertia of the circular cross section with radius R, and r is the distance 
from the neutral axis, as shown in figure C.4. 
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y 



h— b — H 
(a) Cross section 

Figure C.2.- 



Shear stress and shear strain due to shear force. 


y 




Cross section 

Figure C.3.- Shear distribution due to bending moment. 



Cross section 



T 


Shear stress 



r 


Figure C.4.- Shear stress distribution due to torque. 
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Appendix D 

Sign Conventions 

Sign Conventions for Variables 

Uniform sign conventions are useful in solving structural analysis problems. This appendix summarizes the sign 
conventions that are followed in this report. However, one can analyze these problems without following the sign 
conventions given in this appendix. Consistent sign conventions are sufficient for the solution of structural mechanics 
problems. 

The sign conventions are illustrated for common structural mechanics force and displacement variables — such as 
the normal force F, shear force V, bending moment M, and torque T — and the associated displacements — the axial 
displacement u, transverse displacement v, angle of flexural rotation 0, and angle of twist <p. A Cartesian coordinate 
system with the origin 0 and with x, y, and z as the axes (as depicted in fig. D.l(a)) is used to define the sign conven- 
tions. 

Axial force, also referred to as the normal force F is considered to be positive along the positive direction. For a 
beam oriented along the x-axis, the axial force F is positive when it is directed along the positive x-axis, as shown in 
figure D.l(a). The axial displacement u is considered to be positive along the positive x-direction. In other words, 
both the axial force and the axial displacements are positive along the positive x-axis. 

Likewise, the shear force V and transverse displacement v are considered to be positive along the positive direc- 
tion. For a beam oriented along the x-axis, the shear force and transverse displacements are positive along the posi- 
tive y-axis, as shown in figure D.l(b). 

The bending moment and torque, along with the associated displacements (0, cp) are considered to be positive 
along their vectorial directions. For a beam oriented along the x-axis, the moment M is positive when the moment 
vector is directed along the positive z-axis, as shown in figure D.l(c). The rotation 0, which follows the sign conven- 
tion for the moment M, is positive when the rotational vector is directed along the z-axis. 

For a beam oriented along the x-axis, the torque T is positive when the torque vector is directed along the posi- 
tive x-axis, as shown in figure D.l(d). The angle of twist q>, which follows the sign convention for the torque T, is 
positive when the angle of twist vector is directed along the x-axis. 

In the solution of problems utilizing conditions of symmetry, it can be advantageous to select the beam (or x ) 
axis directed from right to left, as shown in figure D.l(e). The other two axes ( J, z ) are selected according to Carte- 
sian conventions, as depicted in figure D. 1 (e). After the x -, and z -axes are defined, the positive directions for forces 
and displacements follow the earlier conventions. 

For example, for a beam oriented along the positive x -axis, the moment M is considered to be positive when the 
vector is directed along the z -axis. On the beam, the moment takes a clockwise direction (see fig. D.l(e)), which is 
opposite to that shown in figure D.l(c). Likewise, the sign for torque is shown in figure D.l(f). Readers can avoid 
confusion by following consistent conventions — that is, by considering a force (or a displacement) variable to be 
positive when it is directed along the positive axis, and by remembering the vectorial directions for moment and 
torque. 

Typically, structural problems deal with gravity load (P ), which is directed along the negative y-axis. However, 
displacement ( X ) is considered to be positive along the y-axis. The convention results in a negative work term, 

W =~ P y X y 

Traditionally, in the solution of elementary structural mechanics problems, a free-body diagram with specified 
signs for forces, moments, and torque is employed. Upon solution, the positive value for the force variable confirms 
the specified sign, whereas the negative value requires the opposite of the assumed sign. Displacements, rotations, 
and angles of twist follow the sign convention for force, moment, and torque, respectively. 
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z 


(a) Axial force ( F) is positive when directed along the positive x-axis. 



z 


(b) Shear force (1/) is positive when directed along the positive y-axis. 



z 


(c) Bending moment (M\ is positive when directed along the positive z-axis. 



(d) Torsional moment (7) is positive when directed along the positive x-axis. 



(e) Moment (M) is positive when directed along the positive z-axis. 



(f) Torsional moment (7) is positive when directed along the positive x-axis. 
Figure D.1 .- Sign conventions for force, moment, and torque. 


NASA/TP— 2004-207430 


162 




Sign Conventions for Equilibrium Equations 


In the equilibrium equations ([fi]{.F} = \P\), aligning the load vector {P} along the positive axes will produce 
the displacement vector {X} along the positive axes. As an illustration of the sign convention for the equilibrium 
equations, consider a beam with a load P at the center span as shown in figure D.2(a). The transverse EE along the 
positive v-direction can be written as 

-R a -R b = -P 


Likewise, the moment EE at A along the positive z-axis can be written as 


-R B t = ~ 


Pt_ 

2 


In matrix notation, the two EE can be written as 



(D.l) 


In the EE given by equation D. 1 , the load component (~P) is along the positive y-axis and the moment vector (-PC/ 2) 
is directed along the positive z-axis. If there is confusion in the sign convention for the equilibrium equations, the 




Figure D.2.- Sign conventions for equilibrium equations. 
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following technique can be used. Apply a fictitious force or moment along the positive axis. For example, for 
moment EE, a fictitious moment m Q should be applied as shown in figure 3.2(b). The moment EE at A can be 
written as 

PI 

-R B t = m 0 - — 

Then, set the fictitious moment to zero (m (j = 0) to obtain the EE written in the positive direction as 



This technique also can be used to write other equilibrium equations. 
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Appendix E 

Important Symbols 


The following symbols are used in this report. 


A 

IB] 

[**l 

[Cl 

\D\ 

dof 

E 

F 

M 

[G] 

§ 

H o\( x ) 

h 

I,J 

IE 

[J] 

[K] 
l 

M 

7% 

M x ,M y ,M xy 

m 

n 

P 

q 

R 


cross-sectional area 
equilibrium matrix 

equilibrium matrix to back-calculate reactions 

compatibility matrix 

dual matrix of IFMD 

displacement degrees of freedom 

Young’s modulus 

axial, or normal, force 

force degrees of freedom 

symmetrical flexibility matrix 

flexibility coefficient 

Hermite polynomials 

plate thickness 

moment and polar moment of inertia 
internal energy 

deformation coefficient matrix 
stiffness matrix 
length parameter 
moment 

moment function as THix) 
plate-bending moments 

number of displacement variables or equilibrium equations 

number of force variables 

external load 

load intensity 

reaction 
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r 

degree of indeterminacy 

[S] 

IFM governing matrix 

T 

torque; as a superscript represents transpose operation 

AT 

temperature variation 

U 

strain energy 

u c 

complementary energy 

U, V, w 

displacement components 

V 

shear force 

W 

work done by external loads 

X 

displacements 

X 

amount of support settling 

x, y,z 

Cartesian coordinates 

a 

coefficient of thermal expansion 

8 R 

effective initial deformation 

£ 

normal strain 

7 

shear strain 

<P 

angle of twist 

K 

curvature 

0 

rotation 

O 

normal stress 

T 

shear stress 


Acronyms and Initialisms 

cc 

compatibility condition 

DDR 

deformation displacement relation 

EE 

equilibrium equation 

FDR 

force deformation relation 

IFM 

Integrated Force Method 

IFMD 

Dual Integrated force Method 
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